Vol. 7 


meUSTICA 


1957 


SOUND PROPAGATION IN DUCTS 
WITH SIDE BRANCHES 


by W. K. R. Lipperr 


Division of Building Research, Commonwealth Scientific and Industrial Research Organization, 
Melbourne, Australia 


Summary 


The propagation of plane waves in a main conducting duct having a closed side branch 
is represented by formulae and diagrams with the characteristic reflection and transmission 
factors. The practical limits of this representation are discussed in measured examples. 
An effective length of the side branch, which differs from the geometrical length, is derived 
from the measurements, and its usage is shown to extend the range of validity of the for- 
mulae and diagrams. 


Sommaire 


On exprime par des formules et des diagrammes, renfermant les coefficients de réflexion 
et de transmission caractéristiques, la propagation d’ondes planes dans un conduit conduc- 
teur principal ayant un embranchement latéral fermé. On examine 4 partir des résultats 
de mesure les limites pratiques de validité de ce mode d’expression. On déduit de ces ré- 
sultats une valeur efficace de la longueur de l’embranchement latéral, valeur qui est dif- 
férente de la longueur géométrique, et on montre comment on peut l’utiliser pour étendre 
le domaine de validité des formules et diagrammes. 


Zusammenfassung 


Es wird die Ausbreitung ebener Wellen in einem Hauptrohr mit abgeschlossenem Seiten- 
rohr in Formeln und Diagrammen mit Hilfe der charakteristischen Reflexions- und Uber- 
tragungsfaktoren dargestellt. Die praktischen Grenzen dieser Darstellung werden an Hand 
von Mefbeispielen diskutiert. Eine wirksame Lange des Seitenrohres, die sich yon der 
geometrischen Lange unterscheidet, wird von den Mefergebnissen hergeleitet, und es wird 
gezeigt, da} ihre Anwendung den Giiltigkeitsbereich der Formeln und Diagramme erweitert. 
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1. Introduction 


The basic idea of producing interference between 
sound waves propagating in acoustical tubes and 
of suppressing certain frequencies originated with 
Herscuet in 1833 and was first applied by Quincke 
in 1866. Later, when electrical wave filters were 
developed, Srewarr [1], [2], [3], Mason [4], 
Warrzmann and Norrner [5], Linpsay [6] and 
others showed that at low frequencies acoustical 
filters can be devised in analogy to electrical filters. 
As a result there has been an improvement in the 
design of mufflers in automobiles, in air condition- 
ing systems and in other practical applications, and 
textbooks on the theory and practice of acoustical 
filters are now available [7]. 


In recent papers on particular problems of plane 
wave propagation in ducts [8], [9] the present 
author found that there are cases where it is pre- 
ferable to represent the frequency characteristics of 
acoustical filters with specific parameters called 


characteristic reflection and transmission factors. 
It was shown that these parameters are particularly 
useful in certain periodic structures. Apart from this 
the characteristic factors are easily and directly 
measurable, as has been amply demonstrated [10], 
BU Ta ees ie eA 

The purpose of this paper is to apply the new 
method to a fundamental acoustical filter consisting 
of a main conducting tube and a closed side branch. 
The long-wavelength approximation represented by 
the characteristic factors yields clear and informa- 
tive diagrams whose usefulness is illustrated by 
measured examples. The theoretical and experimen- 
tal approach discussed below will help to demon- 
strate the degree of accuracy that can be achieved 
with this kind of first order approximation. It will 
follow from the measurements that the effective 
length of the side branch has to be used in the 
formulae instead of the geometrical length, and the 
method of measuring this effective length will be 
described in detail. 
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2. The characteristic reflection and transmission 


factors for plane waves in a duct with a closed 
side branch 


Fig. 1 shows a main duct of rectangular cross 
section and of unspecified length with height a and 
width 6, and a rectangular side branch with the 
cross sectional dimensions d (depth in direction of 
the main duct) and e (width). The geometrical 
length L of the side branch is the distance between 
the opening plane in the main duct wall and the 
terminating piston at the end of the side branch. 
All walls of the main duct and of the side branch 
are supposed to be hard and rigid. 
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pressure of the transmitted wave at the output re-— 
ference plane and the incident wave at the input 
plane, if the output duct_is non-reflecting. 

An easy way of deriving R and T is by electrical 
analogy [13]. Fig. 2a shows an electrical circuit 
between transmission lines the characteristic im- 
pedance of which is Zp) and the waves propagating 
along those lines have the wavelength /. The circuit 
of Fig. 2a can be made equivalent to the acoustical 
duct of Fig. 1 if electrical voltages and currents at 
the various terminals are taken to represent sound | 
pressures and yolume velocities (particle velocity 
area of cross section) at the corresponding re- 
ference planes. 
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Fig. 2. (a) Electrical analogy to the acoustical duct system of Fig. 1, 


(b) simplified analogy and (c) the same for the special case 


Fig.1. Main conducting duct with a 
closed side branch both of rect- 
angular cross-section specifying 
the various dimensions and the 
positions of input (i, i’) and out- 
put (0,0’) reference planes. 


By restricting the wave propagation to plane 
waves it is assumed that the wavelength / is always 
greater than 2a@ or 2b or 2d or 2e. The position 
of the input reference plane in the main conducting 
duct is specified by i,i’ and that of the output re- 
ference plane by 0,0. Both are perpendicular to 
the axis of the main duct, at a distance / from each 
other and symmetrical with respect to the side 
branch. The output duct being non-reflecting, plane 
waves propagating along the main duct result in a 
standing wave in the input and in a progressive 
wave in the output. 

The symmetrical filter section between the reference 
planes will be specified by the characteristic factors 
R and T [10]. The characteristic reflection factor 
R= Ae!* is the complex ratio of the sound pressure 
of the reflected wave and the incident wave at the 
input reference plane, and the characteristic trans- 
mission factor 7=Be!? is the ratio of the sound 


1=0 


The voltages at the input and output terminals 
are E, and EF, and the incident voltage wave is 
specified by 


Eo ed Qu ylh 


where y is the co-ordinate along the input trans- 
mission line with y=O0 at i,i’, and the outgoing 
wave by 

E, et j2xah 


where x is the co-ordinate along the output trans- 
mission line with x=0 at 0,0. 
At the input terminal 
E, 
E,= 2 +k 
2 
where Ep is the reflected wave. This yields the 
characteristic factors as specified above 


rogfriy “Tio ee 
E E 


0 0 


(1), (2) 


Applying Tueveniy’s theorem to the input transmis- 
sion line at the terminals i,i’ this line can be re- 
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presented by a generator with the generator voltage 
E, (voltage at i,i’ if the line is open at the input) 
and supply resistance Z). The non-reflecting output 
transmission line is equivalent to a load resistance 
Z, at 0,0. This representation of the input and out- 
put transmission lines by lumped circuit elements 
in Fig. 2b is convenient for computing R and T 
from eqs. (1) and (2). For actually measuring R 
and 7, however, it is suggested that the representa- 
tion of Fig. 2a and the method described elsewhere 
[10] be used. 

At this stage it is possible to take advantage of a 
specific feature of the characteristic factors. These 
factors transform easily for different values of I or 
for different positions of the reference planes along 
the input and output transmission lines, a piece of 
which (J/2 in length) is circuited in series to the 
input and output in Fig. 2a and.b. 

Let the characteristic factors for 1=0 be 


Ry=Ajei* and T,=Byei% 
then [10] 
A=4A,, B=B, (3) 
and 
a= Pach ae eae kg (4) 
A A 


This means that R and T can be easily found from 
R, and T, so that it is possible to deal first with 
the simpler duct system with /=0, the equivalent 
circuit of which is given in Fig. 2c. At first sight 
this may appear to introduce inaccuracies, as when 
l=0 the wave pattern at the reference planes is 
greatly disturbed by higher order modes caused by 
the sudden change in the wall impedance at the side 
branch. However, in the acoustical case, by defini- 
tion, the characteristic factors are specified and 
measured from the standing and progressive waves 
in the input and output ducts, which are free of 
disturbances from higher order modes. Therefore, 
the actual disturbance of the sound field at the re- 
ference planes when /=0 is trifling, as long as the 
position of the reference planes is always trans- 
formed (with eqs. (3) and (4) ) to parts of the in- 
put and output ducts which are free of disturbances, 
before any connections to other filter systems or 
sections are considered. 
In the electrical circuit of Fig. 2c the impedance 
_Z, and the shunt impedance Z’ represent the acous- 
tical characteristic impedance of the main duct and 
the input impedance of the side branch duct, which 
are related: 


aaa 5 
25 . (5) 

and i 
U = —jF cot 2a (6) 
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o density of medium (air), 

ec sound velocity, 

S cross section of the main duct, 

S’ cross section of the side branch, 
L’ effective length of the side branch, 
A wavelength. 


In eq. (6) an effective length L’ has been chosen 
instead of the geometrical length L. A justification 
for this and a better specification of L’ will be given 
in connection with measurements discussed later. 

Ry and T, can be found from Fig. 2b and with 
eqs. (1), (2), (5) and (6) 


ig otan 2a 
Ro 7 Topas (7) 
acy eta ye 
25 A 
1 
Y paket ern 8 
: Fee Sian oa 
OES A 
Eqs. (7) and (8) are equivalent to 
; tan 20 
> ae ay ——- (9) 
1 eng) (we 2 
/ at Z ( 4 tan* 22 F 
oh == Sx —arctan (5 pun 2a) (10) 
and 
ee an 
/ ay ne 
eden tan* 22 7 
; 1 S i 
= —arct tates) < 1 
Bo are an (52 an2nx 4 (12) 


These formulae have a simple form and are very 
suitable for experimental checking, as will be shown 
later. 

R, and T, belong to a loss-free filter so that the 
following conditions are fulfilled [10] 


2 pens pie ATE 
Ay? + By =1 a —fyp= +t = 


which means that eqs. (9), (10) and (11), (12) are 
interrelated and that it is sufficient to give a dia- 
gram for one set of equations only. 

In Fig. 3a By and in Fig. 3b fy has been plotted 
against 2 L’/A for various values of S’/S. Only one 
cycle for 0 <2L’/A <1 has been drawn as the 
values of By and of fy are repeated if 2 L’/A is in- 
creased by an integer (V=0, 1, 2, 3,...). 

The usefulness of the diagrams of Fig. 3 is best 
illustrated by examples, and some results of actual 
measurements are reported below. 
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Fig.3. Magnitudes (By) and phase values (f/f) of the 
characteristic transmission factors of a duct sys- 
tem of Fig. 1 with /=0 plotted against 2 L’/A for 
various ratios S’/S. 


3. The measured duct system 


The duct system used for the measurements is 
shown in Fig. 4a. The main conducting duct and 
the side branch have the same square cross section. 
From the specification of Fig. 1 this means that 
a=b=d=e. Fig. 5 shows a photograph of the duct 
system. 

The main duct is made of brass with a square 
cross section of 7.6 cm (inside) and wall thickness 
of 0.95 cm, Other dimensions of the system are 
shown in Fig. 4, Part of the main duct, 25 cm in 
length, consists of hardwood having a square open- 
ing (7.6 cm X 7.6 cm) in the top wall. 


The side branch is made of brass with the same 
cross sectional measurements as the main duct and 
with angle pieces of brass at the ends. The side 
branch duct is put on top of the main duct so that 
the square opening and the cross section of the 
side branch fit accurately. A thin rubber lining 
(0.02 cm thick) has been applied on the outer top 
wall of the main duct to ensure a complete seal of 
the joint with the side branch. Two clamps hold the 
side branch in position. Within the side branch duct, 
the walls of which have been carefully machined, a 
brass piston with a wooden extension can be moved 
up and down. A good seal was obtained by making 
the piston sides slightly (0.01 to 0.02 cm) smaller 
than those of the inner tube walls and by inserting 
long strips of shim brass of appropriate thickness 
between the walls and the piston. 
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The length L of the side branch can be varied 
continuously up to L=62.7 cm. (This value of L is 
obtained by removing the piston completely and 
terminating the side branch end with a metal plate.) 
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Measurements are in cm 
@) Orifice 
=RaSsss 


3 
bo 63 


~<a. 00m === 
to microphone unit 

Fig. 4. Diagram showing a cross section (a) of the 
measured side branch system and (b) of the 
orifice part of the probe tube. 


Fig. 5. Photograph of the measured side branch system. 
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4. Equipment and measuring technique 


The equipment for measuring the characteristic 
factors, including the non-reflecting terminal, micro- 
phone unit and other electrical parts, is the same as 
described in previous papers [8], [12]. By taking 
advantage of the fact that in the present acoustical 
filter systems the input and output sound fields can 
be recorded continuously and in one probing run, 
the probe tube of the microphone unit has been ex- 
tended to about 300 cm. The dimensions of the steel 
probe tube are shown in Fig. 4b. The microphone 
unit is mounted on a trolley which can be moved by 
a motor and pulley system at constant speed along 
a steel ruler. The probe tube which is attached to 
the microphone may thus be moved slowly along 
the bottom of the main duct with the position of the 
orifice varying continuously between distances of 
about 70cm in the output side (negative x co-ordi- 
nates) to about 100 cm in the input side (negative 
y co-ordinates). During each run the frequency was 
kept constant and the sound pressure recorded with 
a high-speed level recorder. The probing procedure 
can also be done by hand for accurately measuring 
the minima positions and the position of equal 
phase in the standing and progressive waves. 


Fig. 6 shows some records at various frequencies 
with a recording time of about 100s for each re- 
cord. In Figs. 6a to d the length of the side branch 
is L=62.7¢cm and the frequencies are 400, 700, 
1200 and 1600 c/s respectively. In Fig.6e L =17 cm 
and the frequency is 2200 c/s. These records were 
all taken with a 50 dB potentiometer of the level re- 
-corder. The position in the main duct is specified 
by a few distance marks at the top of each record, 
corresponding to x= —50cm and x=Ocm in the 
output duct and to y=Ocm, y= —50cm and 
y= —100cm in the input duct. The positions 
x=0 and y=0 for the output and input reference 
planes coincide in these cases, as /=0. 
These examples show that, in general, there is a 
‘standing wave in the input duct and a progressive 
‘wave in the output duct, with a disturbance caused 
by higher order modes in that part of the main duct 
Sihére the side duct branches off. This disturbance 
is restricted, however, to a short distance in the in- 
put and output ducts and the characteristic trans- 
mission and reflection factors can be easily deter- 
mined from the undisturbed wave patterns as has 
been described in earlier papers’ [8], [10], [11]. 
e records (c) to (e) of Fig. 6 are typical examp- 
les to demonstrate how greatly the wave pattern in 
he duct systems measured can vary. In record (c) 
the ratio of minimum to maximum sound pressure 
in the input (Piin/Pmax) is about — 28 dB or 0.04, 
and the ratio of transmitted sound pressure in the 
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Fig.6. Examples of records (with a 50 dB _ potentio- 
meter) showing the sound pressure distribution 
in the main duct at various frequencies (a) 400, 
(b) 700, (c) 1200, (d) 1600 and (e) 2200 c/s; 
L=62.7 cm in (a) to (d) and L=17 cm in (e). 


output and maximum sound pressure in the input 
(Ptr/Pmax) is about — 29 dB or 0.035, which yields 
the values By)=0.07 and A,=0.92 [10]. This 
means that at this particular frequency (1200 c/s) 
the single side branch has an attenuating effect on 
the sound propagation comparable with the trans- 
mission loss of a light wall. 


The opposite can be seen with the record of 
Fig.6d, which was measured at 1600c/s and 
where the attenuating effect of the side branch is 
negligible. In this case the reflection factor Ag is 
practically zero which can be seen from the fact 
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that the standing input wave has degenerated into a gas x AY +Ar o, ; (15) 
progressive one. A 2 


In cases with no minima the recording of the 
sound pressure magnitudes alone does not permit 
direct measurement of the wavelength. In the tech- 
nique applied here, however, the wavelength (and 
not the frequency) is the fundamental quantity 
which must be determined accurately, and actually 
in such cases where there is no standing wave not 
only the wavelength, but also the phase value (/) 
of the transmission factor can be measured directly, 
i. e. without knowing the frequency and phase velo- 
city accurately. This will be explained briefly for 
symmetrical filters and an arbitrary position of the 
reference planes, as it is of general interest in acous- 
tical measuring technique. 

In cases where there are no minima in the input 
duct the orifice of the probe tube can be placed at 
an arbitrary position in the undisturbed part of the 
progressive wave of the input duct, and the distance 
between this position and the input reference plane 
is called dz. The phase value of the sound pressure 
at this position denoted by a (+) sign is measured 
in the usual way by making the phase and magni- 
tude of the microphone voltage equal to those of a 
reference voltage (e.g. of the oscillator voltage) 
with the aid of a phase shifter and an attenuator. 
This can be done with high degree of accuracy by 
checking that the difference of the microphone and 
reference voltages is practically zero [8]. By mov- 
ing the probe tube along the input duct (without 
changing the phase adjustments of the reference 
voltage) other positions with the same phase value 
(+) can then be measured at distances Az+ NV 4 
and those with the phase value (—) (which is x 
different from (+)) at distances Az+ (N+ 5) he 
with Vi = 0,90 9203, 3... 

By measuring various distances which belong to 
different (+) or (—) phase values and numbers 
N, the wavelength 2 can be computed from the cor- 
responding differences, which are always a multiple 
of half a wavelength. 

The first position in the output duct which has the 
phase value (—) is then measured in the same way 
and its distance from the output reference plane is 
called Ax. By keeping in mind that the distance Az 
specified above for A =0 and the usual distance Ay 
between the input reference plane and the first pres- 
sure minimum of the standing wave in the input 
duct for A == 0 are simply related by 


Az = Ay + at (14) 


the phase value / of the characteristic transmission 
factor can then be computed with the earlier for- 
mula for general symmetrical filters [10]: 


That this formula is correct in the present case for 
A=0 follows from the fact that, although Az has 
been chosen arbitrarily in the experimental proce- 
dure above, all positions of equal (+) or (—) 
phase in the progressive input and output waves 
shift with changing Az in such a way that the sum 
Ay + Ax does not change. Thus, from eq. (15) the 
value of f is not changed. 

On the other side, the phase value @ of the charac- 
teristic reflection factor, which, after [10], depends 
essentially on Ay , is not measurable in cases where 
A=0O; it is arbitrary because of R= Aei*=0. 

As demonstrated above 4 and f are still directly 
and accurately measurable when A=O. This rule 
applies, to some extent to cases where A~0 which 
in loss-free cases can be very helpful for determin- 
ing a by measuring f/ first and computing a from 
the general phase relation for loss-free symmetrical 
filters (eq. (13)). 

Another typical record is shown in Fig. 6 e. It can 
be seen from this that the disturbance of the wave 
pattern from higher order modes is greater and ex- 
tends to longer distances in the input and output 
parts of the main duct compared with examples at 
lower frequencies. The frequency in this case was 
2200 c/s or the corresponding value of the dimen- 
sionless quantity O/m=2a/A was O/x=0.95 which 
is rather close to the cut-off frequency (with the 
value @/%=1) below which only plane waves and 
no higher order modes can be propagated. This 
indicates that in a frequency range just below the 
cut-off frequency the simple long-wavelength theory 
of eqs. (9) to (12) has to be used cautiously and 
is likely to fail. It is, therefore, of practical impor- 
tance to determine the limits of applicability of 
these formulae and of the corresponding diagrams 
of Fig. 3. This will be done by discussing some 
measured examples. j 


5. Results of measurements 


The characteristic reflection and transmission fac- 
tors of plane waves travelling along the main duct 
(illustrated in Fig. 4a and 5) have been measured 
at various frequencies and for different lengths L of 
the side branch. In these measurements the method 
of the non-reflecting terminal described above and 
elsewhere [8], [10] has been applied. In Fig. 7, 8 
and 9 the measured values of A), By and of a, fy 
have been plotted against O/x (=2a/A) for three 
examples with L =9.4cem, L = 17cm and L = 32.2cm 
(The subscript zero in the measured characteristic 
factors refers to the position of the reference pla- 
nes specified by 1=0 as described above.) 


ho : 
ACUSTICA 
a, (1957) W. K.R.LIPPERT: SOUND PROPAGATION IN DUCTS 143 
1.0 
08 AE 
bo | 
os 
30 
02 
40 
| -04 + 
p08 ia a k-0 | 
Pe, xt | Oo 
= =F a 4+— ey = 
2 T le Ny rn 
B16 ASS Betas ‘ 
20 ane ene < 08 10 
0 02 04 06 08 10 ; 


Ca 


\Fig. 7. The measured magnitude (Ay and By) and phase 

: values (% and fo) of the characteristic reflec- 
tion and transmission factors plotted against 
@/x=2.a/A for L=9.4 cm and the correspond- 
ing computed curves for L’=10.37 cm. 


ig. 8. The measured magnitude (A, and By) and phase 
values (a) and fio) of the characteristic reflec- 
tion and transmission factors plotted against O/x 
for L=17 cm, and the corresponding computed 
curves for L’=17.91 cm. 


The measured quantities can be computed also 
using eqs. (9) to (12) or represented by the 
agrams of Fig. 3. The parameter is S’/S=1. By 
making L’=L it was found at first that the com- 


~~ 


Fig. 9. The measured magnitude (Ay and By) and phase 
values (a) and fo) of the characteristic reflection 
and transmission factors plotted against O/x for 
L=32.2 cm, and the corresponding computed 
curves for L’=33.62 cm. 


puted curves did not agree with the measured ones, 
especially at the minima of By, (and Ay). Values of 
L’ were then chosen so that there was good agree- 
ment between measured and computed curves at the 
minima and in the lower and medium frequency 
range. This could be achieved when L’ was slightly 
greater than L. In the examples quoted above the 
values for L’ providing best agreement between 
measured and computed curves were: in Fig. 7 
L’ = 10.37 cm, in Fig. 8 L’ = 17.91 cm and in Fig. 9 
L’ = 33.62 cm. With these values of L’ the theore- 
tical curves of By and fy have been plotted in 
Figs. 7 to 9 with the diagrams of Fig. 3 and also 
those curves of Ay and a with eq. (13) or eqs. (9) 
and (10). 

There is good agreement between measurements 
and theory in a wide frequency range between at 
least 0< O/x <0.5. In Fig. 9 even several cycles 
(with V=0, 1 and 2 of Fig. 3) of the magnitude 
and phase curves are repeated very well. In the 
range where 0.5<0/nx<1 the difference between 
the measured and computed curves becomes greater 
with increasing frequency, as can be seen from the 
measured examples, and the approximating eqs. (9) 
to (12) become less accurate the closer 9/x ap- 
proaches the value unity, which is not surprising 


(cf. Fig. 6e). 


In the frequency range, where the approximating 
theory is valid, the accuracy of the computed values 
depends very much on the knowledge of the effec- 
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tive length L’ especially at the frequencies where 


minima of By occur. To get a clearer picture of how 


L’ depends upon L, more measurements were carried 
out with different length Z and the effective lengths 
L’ were determined as above, from the @/z values 
where the first minima of By occur or from the best 
fit between measured and computed curves of By 
(or Ay) plotted against O/z. 
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Fig. 10. The measured ratios (L’—L)/d plotted against 
L/d for a side branch system with the same 
square cross section of both main and side 
branch duct. 


In Fig. 10 the measured ratio (L’ — L)/d has been 
plotied against the ratio L/d where d=7.6 cm 
(3 in) is the width of the square side branch, 
which in this particular case is equal to that of the 
main duct (Fig. 4a). With increasing L/d the curve 
of the measured values in Fig. 10 at first shows a 
maximum and a minimum and finally reaches a 
constant value of (L’—L)/d ~0.2 when L/d>4. 
This curve provides a quick means of finding L’ 
if L and d are given. Fig. 10 refers to the special 
case where d=a and b=e=a. However, similar 
curves for finding L’ can be measured when d= a 


and b=: e (Fig. 1) or when the cross section of the 
side branch is other than rectangular, e. g. is circu- 
lar. This method of determining L’ 
superior to the rough guess of the effective length, 


is obviously 


e.g. L’=L + (a/2), which has been recommended 


iat 


6. Discussion 


The results of the measurements discussed above 
constitute an experimental study which demonstra- 
tes the validity of the approximating equations (9) 
to (12) and of the diagrams of Fig.3. The results 
obtained will vary with different shapes and di- 
mensions of the cross section of the main duct and 
of side branches. However, even in cases of varied 
experimental conditions the particular examples of 
the present study are informative by showing what 
sort of agreement can be expected between theory 
and practice, and the experimental procedure ap- 
plicable. 
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The measurements indicate that accurate deter- 
mination of L’ (e.g. in the experimental way dis- 
cussed above) leads to closer agreement over a much 
greater frequency range between experimental re- 
sults and the long-wavelength theory. 

Eqs. (9) to (12) show clearly how the charac- 
teristic factors depend upon the parameter S’/S. 
The particular shape of S’ or S as well as the length 
of the side branch have an effect on the values of 
L’ —L or of L’ and thus on the positions of minima 
and maxima in the frequency characteristics. 

_In cases of several equal side branches around 
the circumference of the main duct the various 
cross sections of the side branches have to be added 
up to form S’. That a summarizing effect of S’ 
of many side branches (with a single cross sectional 
area small compared with that of the main duct) 
really occurs has been shown by the author in a 
former paper on measurements of sound attentua- 
tion in a duct the walls of which consisted of per- 
forated bricks covered by plywood at the outside 
[14]. 

When a side branch system of Fig. 1 is used as 
an element of an acoustical filter, e.g. by connect- 
ing the output to another discontinuity or by re- 
peating the same side branch at equal distances 
along the main duct, the values of dy, By and 
%, fo have to be transformed to A, B and a, f 
which belong to a position of the input and output 
reference planes symmetrical to the side branch and 
at a distance / apart from each other as shown in 
Fig. 1. This transformation of the characteristic 
factors can be done most conveniently with eqs. (3) 
and (4). The value / has to be made sufficiently 
large, so that the sound fields at the input and out- 
put reference planes are reasonably free from dis- 
turbances of higher order modes (cf. Fig. 6). The 
value of J usually has a great effect on the fre- 
quency characteristics of multi-sectional acoustical 
filters. 

In cases of n-sectional filters each section of which 
consists of a loss-free symmetrical filter as in the 
present case e.g. of a side branch section of Fig. 1, 
the characteristic transmission factor T,,—=B, e/? 
can be found with general diagrams from the values 
B and f of one section, as has been shown else- 
where [9]. There is a practical advantage in the 
use of such diagrams. From eqs. (3) and (4) it 
follows that § alone depends upon / but B does not. 
Thus the effect of the various parameters S/S, L’, 
! and n on the form and positions of maxima and 
minima (or of pass- and stop-bands) of the final 
frequency characteristics can be separated and con- 
sidered step by step. The frequency characteristics 
can be described completely, and specific values of 
various parameters may be selected to meet practical 
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design requirements. This is illustrated in the follow- 
ing example. 

A practical problem is to design an n-sectional 
acoustical side branch filter (with a practically non- 
reflecting termination of the main duct) so that an 
attenuation as large as possible is obtained in a 


| frequency range with a given middle frequency 
fst =c/4 (where A is the wavelength in the duct 


at f,,). At first, it is obvious that the number of 
sections n and the parameter S’/S (cf. Fig. 3a) 
should be chosen as large as practically possible. 
From the general diagrams of n-sectional filters 
(Fig. 3 to 5 in [9]) (where B, is represented by 
B and 3) it follows that B must have a small value 
and / should equal + Na at f, to make the value 
of B, small at the particular frequency and in 


the neighbouring frequency range. From Fig. 3a 


B,=B=0 is obtained when 2 L’/A,, = : or when 


L’ =1,/4 in which case also it follows from Fig. 3b 
that {y= +x2/2. The required value of / can now 


_be determined with eq. (4) by taking e.g. B= —a, 


By = —2/2 and A=1,,=4L’. This yields 1= L’ =1,4/4 
which is the well known condition for large attenua- 
tion of a side branch filter at a particular frequency 
and in the neighbouring frequency range [7]. The 
values of n and S’/S being specified the frequency 
characteristics of B, (and /,) can then be determin- 
ed completely with diagrams in the way discussed 
above. 


7. Conclusions 


Formulae and diagrams for representing the cha- 
racteristic reflection and transmission factors of plane 
sound waves propagating along a main conducting 
duct with a closed side branch have been derived. 
The practical limits of the long-wavelength approxi- 
mation are discussed experimentally for a particular 
side branch system having the same cross-section 


_ area as the main duct. The use of an accurate value 
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for the effective length of the side branch is shown 
to be important and a practical method of measuring 
this length is demonstrated. The technique adopted 
in this paper is considered to be practicable and 
useful in the design of acoustical filters with side 
branch systems, e.g. of n-sectional periodic struc- 
tures. (Received 15th January, 1956.) 
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SCHWINGUNGSUNTERSUCHUNGEN 


von H. Harske 


III. Physikalisches Institut der Universitat Géttingen 


Zusammenfassung 


1. Es wurden die Beschleunigungen ermittelt, die durch die Uberfahrt eines Zuges im 
Schienenkopf erzeugt werden. Im untersuchten Frequenzbereich von 10 bis 100 kHz liegen 
diese in der GréRenordnung yon 10 m/s? je kHz Bandbreite. Die Spektren fallen mit wach- 
sender Frequenz nur allmahlich ab; bei der waagerechten Komponente ist diesem Abfall 
eine Uberhéhung von 8 dB im Bereich von 35 bis etwa 60 kHz iiberlagert. 

Zur Schwingungsanregung der Schiene wurden zwei Modellversuche durchgefiihrt; aus 
ihnen geht hervor, daB der Verlauf des Beschleunigungsspektrums nur bedingt von der Art 
der Anregung abhanet. 

2. Es wurden die Schwingungseigenschaften der Schiene untersucht. Bei der Messung 
der Wellenlingen zeigte sich — in Ubereinstimmung mit friiheren Untersuchungen von 
Naaxe — das Auftreten neuer Dispersionsaste bei etwa 35 und 60 kHz. Ferner wurde bei 
sinusformiger Schwingungsanregung die resultierende Amplitudenverteilung lings der 
Schiene aufgenommen. Die erhaltenen Diagramme zeigen je nach Frequenz sehr verschie- 
dene Gestalt; es treten groBe Welligkeiten auf, die durch stehende Wellen wie auch durch 
Superposition mehrerer Wellentypen verursacht werden. Die resultierende Amplitudenver- 
teilung hangt stark vom Ort der Anregung ab. 

Die Dampfung der Wellenausbreitung im Schienenkopf betragt im untersuchten Fre- 
quenzbereich meist weniger als 0,3 dB/m. Die Untersuchung von SchweiSstellen ergab fiir 
AbbrennstumpfschweiRungen kaum eine Dampfung (<1 dB), Thermit-SchweiSstellen hin- 
gegen dampfen die durchlaufende Welle bis zu 6 dB. 


Summary 


1. The accelerations in the frequency range 10 to 100 ke/s experienced by the rail dur- 
ing the passage of a train are measured. The spectrum falls off slowly with rising fre- 
quency; this depends only on the type of excitation (as shown by model research). 

2. The mode of vibration of the rail is investigated. Dispersion regions appear between 
35 and 60kc/s. A large waviness appears along the rail, ascribed to the superposition of 
standing waves of many types, and a function of the place of excitation. The attenuation 
in this frequency range is less than 0.3dB/m. Electric welds produced an attenuation 
<1 dB, but thermite welds up to 6 dB, for the progressive waves. 


Sommaire 


1. On a déterminé les accélérations qui sont produites dans le champignon d’un rail par 
le passage d’un train. Dans la gamme des fréquences étudiée, soit 10 4 100kHz, ces 
accélérations sont d’un ordre de grandeur de 10 m/s? par kHz de largeur de bande. Les 
spectres ne s’affaiblissent que progressivement quand la fréquence augmente; dans le cas 
des composantes horizontales, il se superpose a cette décroissance une pointe de 8 dB dans 
la gamme allant de 35 a environ 60 kHz. 

Pour exciter la vibration des rails, on a procédé 4 deux essais sur maquette; les résul- 
tats montrent que l’allure du spectre d’accélération n’est déterminée que par le mode 
d’excitation. 

2. On a étudié les caractéristiques vibratoires des rails. La mesure des longueurs d’onde 
a mis en évidence — en bon accord avec les résultats des études antérieures de NAAKE — 
de nouvelles branches de dispersion 4 environ 35 et 60 kHz. On a relevé en outre, dans le 
cas d’une excitation de forme sinusoidale, la distribution des amplitudes résultantes le long 
du rail. Le diagramme obtenu a une forme qui différe beaucoup suivant la fréquence; il 
se produit des ondulations importantes qui sont dues a des ondes stationnaires ainsi qu’a 
la superposition de plusieurs types d’ondes. La distribution résultante de l’amplitude 
dépend fortement du point ou se fait l’excitation. 

L’amortissement des ondes dans le champignon du rail est le plus généralement infé- 
rieur a 0,3 dB/m dans la gamme des fréquences étudiées. Les régions de soudure ne pro- 
duisent qu’un amortissement insignifiant dans le cas de joints soudés électriquement 
(<14dB). Les soudures thermiques amortissent par contre l’onde jusqu’a 6 dB. 
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Einleitung 

t Die seinerzeit von Naaxe [1] durchgefiihrten 
| Schwingungsuntersuchungen ergaben unter anderem 

| den iiberraschenden Sachverhalt, da8 im Tonfre- 
ce die Dampfung einer in der Schiene 

entlanglaufenden Welle mit wachsender Frequenz 
| abnimmt. Es war daher von Interesse, das Verhal- 
ten der Schiene im Ultraschallbereich kennen zu ler- 
' hen. “In neuen Untersuchungen wurden deshalb die 
_ Schwingungseigenschaften der Schiene im Frequenz- 
_bereich von 10 bis 100kHz ermittelt und die Be- 


; schleunigungen gemessen, die durch die Uberfahrt 
B cines Zuges in der Schiene erzeugt werden. 


1. Aufnahme und Spektralanalyse des 
Uberfahrtgeriusches — 


5 1.1. Die MeBmethode 


Es wurden die senkrechten und die waagerechten 
Transversalbeschleunigungen ermittelt. Zu diesem 
_ Zweck wurde am Kopf der Schiene ein magnetischer 
_ K6rperschall-Empfanger befestigt, dessen abgege- 
_ bene Spannung nach Fourier analysiert wurde, 
um die Beschleunigung als Funktion der Frequenz 
zu erhalten, Zur Analyse diente ein Uberlagerungs- 
_ empfinger mit einer Bandbreite von 80 Hz; die Um- 
wandlung der der Schnelle proportionalen Mikro- 
phonspannung in eine der Beschleunigung propor- 
tionale Eingangsspannung fir den Analysator 
wurde durch einen entsprechenden Frequenzgang 
des Vorverstarkers erzielt. 

Die zur Analyse des gesamten Bereiches von 10 
bis 100 kHz erforderliche Zeit betrug etwa 2,5 Mi- 
nuten. Da eine Speicherung des Gerdusches wegen 
der hohen Frequenzen nicht moéglich war, muften 
die an mehreren Ziigen gewonnenen Mefergebnisse 
zusammengesetzt werden; um den ganzen Frequenz- 
bereich zu iiberdecken. Es wurde dazu in Vorversu- 
- chen sichergestellt, da waihrend der Uberfahrt eines 
Zuges eine beliebige Frequenz zu jeder Zeit mit etwa 
gleicher Amplitude angeregt wird und dal gleich- 
artige Ziige mit gleicher Fahrgeschwindigkeit im 
Mittel dieselben Beschleunigungen hervorrufen. Die 
gefundenen Schwankungen betragen 2 bis 5 dB und 
liegen damit im Bereich der bei derartigen Messun- 
gen zu erwartenden Ungenauigkeiten. 


1.2. Die Mebergebnisse 


Die Analyse wurde fiir Giiterziige und fir 
Schnellziige durchgefiihrt. In Abb. 1 sind zwei Bei- 
-spiele fiir die erhaltenen Spektren wiedergegeben. 
Um den Einflu8 von Zufalligkeiten oder verschiede- 
ner Wagentypen zu eliminieren, wurde tiber mehrere 
derartige Spektrogramme gemittelt. Daraus ergaben 
sich die in Abb. 2 gezeigten Diagramme. Es treten 
demnach im Ultraschallbereich noch _betrachtliche 
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Abb. 1. Beispiele fiir die bei der Analyse erhaltenen 
Pegelschreiberstreifen; (a) senkrechte, (b) 
waagerechte Beschleunigungen «des Schienen- 
kopfes. 
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Abb. 2. Die durch die Uberfahrt eines Zuges im Schie- 
nenkopf hervorgerufenen Beschleunigungen 
als Funktion der Frequenz; (a) senkrechte, 
(b) waagerechte Beschleunigungen. Die mit 
»D“ gekennzeichneten Kurven wurden bei 
Schnellziigen, die mit ,,G“ bezeichneten bei 
Giiterziigen gemessen. Der Bezugspegel 0 dB 
entspricht einer Beschleunigung von 0,1 m/s? 
je 80 Hz Bandbreite. 


Beschleunigungen auf, die in der GroSenordnung 
von 10 m/s? je kHz Bandbreite liegen. Alle vier Dia- 
gramme zeigen im Prinzip einen ahnlichen, mit zu- 
nehmender Frequenz leicht abfallenden Verlauf, nur 
liegen die Pegel bei den D-Ziigen infolge gréSerer 
Fahrgeschwindigkeit hoher als bei den Giiterziigen 
(95 bzw. 55 km/h). Ein Vergleich der senkrechten 
mit den waagerechten Beschleunigungen zeigt jedoch 
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bei letzteren auBer einem — leicht verstandlichen — 
geringeren Gesamtpegel eine auffallige Uberhohung 
von 8 bis 9 dB, die bei etwa 35 kHz beginnt und 
sich bis 60 kHz hinzieht. 


1.3. Modellversuche 


Die auf vorstehend beschriebene Weise erhalte- 
nen Spektren unterscheiden sich bei Ziigen mit ver- 
schiedenen Wagentypen nur geringfiigig. Und da 
selbst die Diagramme von Giiterziigen und von 
Schnellziigen (mit ihrer etwa doppelten Fahrge- 
schwindigkeit!) einen sehr ahnlichen Verlauf zei- 
gen, erhebt sich die Frage, inwieweit das erhaltene 
Spektrum tiberhaupt ein Charakteristikum fiir den 
betreflenden Zug darstellt und wie stark es durch 
die Schwingungseigenschaften der Schiene selbst be- 
dingt wird. Als Beitrag zu dieser Uberlegung wur- 
den zwei Modellversuche durchgefthrt: 

Die Schwingungsanregung der Schiene wahrend 
der Uberfahrt eines Zuges erfolgt durch die nur sehr 
geringe Relativbewegung zwischen Rad und Schiene, 
aber die dabei auftretenden Krafte sind sehr groB 
(Achslast bis 20t!). Im Versuch wurde nun das 
rollende Rad durch eine kleine rotierende Eisen- 
scheibe von 8cm Durchmesser und 2cm Dicke er- 
setzt, die an einem festen Ort gegen die Schiene ge- 
driickt wurde. Die Relativgeschwindigkeit dieser 
schleifenden Scheibe gegen die Schiene betrug etwa 
2m/s, die andriickende Kraft etwa 50 kg. Obwohl 
die Anregungsbedingungen bei diesem Modellver- 
such auferst verschieden waren von den Verhaltnis- 
sen bei der Uberfahrt eines Zuges, zeigen die erhal- 
tenen Beschleunigungsspektren wieder die gleichen 
Merkmale. 

In einem zweiten Modellversuch wurde die An- 
regung der Schiene noch starker idealisiert, indem 
auf eine Relativbewegung zwischen anregendem 
System und Schiene ganz verzichtet wurde. Die 
Schiene wurde durch einen mit einer Rauschspan- 
nung gespeisten Korperschall-Sender angeregt, der 
an verschiedenen Stellen der Schiene befestigt 
wurde. Die auf diese Art erhaltenen Spektren zeigen 
zwar grobe Abweichungen voneinander und_ han- 
gen in ihrem Verlauf stark vom Anregungsort ab, 
allen gemeinsam ist jedoch wieder die Uberhéhung 
des Pegels im Frequenzbereich von 35 bis etwa 
60 kHz, ahnlich wie im ersten Modellversuch oder 
wie in praxi beim Uberfahrtgerausch der Ziige. Das 
erhaltene Beschleunigungsspektrum hangt demnach 
nur bedingt von der Art der Anregung ab. 


2. Untersuchung der Schwingungseigenschaften 
der Schiene 


Die weiteren Untersuchungen hatten zum Ziel, 
die Schwingungseigenschaften der Schiene selbst zu 
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ermitteln. Die Messungen erfolgten an geschweifbten 
Langschienen. Die Schwingungen wurden durch 
einen piezoelektrischen Korperschall-Sender erzeugt, 
als Mikrophon diente ein magnetischer Korperschall- 
Empfanger, der sich in einem kleinen Wagen befand 
und in konstantem Abstand langs der Schiene be- 
wegt werden konnte. 


2.1. Die Wellenlingen auf der schwingenden Schiene 


Zunichst wurden die Wellenlangen fiir die Trans- © 
versalkomponenten der im Schienenkopf entlang- 
laufenden Wellen ermittelt. Die Messung erfolgte in 
tblicher Weise durch Aufsuchen der Orte gleicher 
Phase mit Hilfe von Lissajous-Figuren auf dem 
Schirm eines Kathodenstrahl-Oszillographen. Die 
erhaltenen Werte sind in Abb. 3 dargestellt*. Die 
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Abb. 3. Die Wellenlingen der Transversalkomponenten 
der im Schienenkopf entlanglaufenden Wellen 
als Funktion der Anregungsfrequenz; (a) senk- 
rechte, (b) waagerechte Schwingungskompo- 
nente. 


— 


Kurven zeigen den gleichen Verlauf, wie ihn schon 
Naake gefunden hatte: Der Abfall der Wellenlainge 
mit zunehmender Frequenz ist infolge der frequenz- 
abhangigen Phasengeschwindigkeit nicht rein hyper- 
bolisch. Bei etwa 35 kHz (fiir beide Schwingungs- 
richtungen etwas verschieden) tritt fiir beide Kom- 
ponenten ein neuer Dispersionsast auf, d.h. hier 
schwingt der Schienenkopf nicht mehr gleichformig 
iiber den ganzen Querschnitt. Ein weiterer Ast ist 
an den MeSpunkten bei 60 kHz zu erkennen. 


Die Art der MeBmethode bringt es mit sich, dah 
bei Uberlagerung mehrerer Wellentypen stets nur 
die Wellenlange der Schwingungsform gemessen 
wurde, die besonders stark angeregt war. Daher 


Schiene aus Thomas-Stahl, die Kreuze gelten fiir eine 
aus Siemens-Martin-Stahl. Da die Hiarten beider bis 
auf etwa 2% gleich waren, unterscheiden sich die ermit- 
telten Wellenlangen kaum vyoneinander. 
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kann man den Diagrammen zusitzlich entnehmen, 
daB der Schienenkopf oberhalb 35 kHz kaum mehr 
_ ohne Unterteilung schwingt, daB aber die niachst 
_ héhere Schwingungsform (ab 60 kHz) fast gar nicht 
 angeregt wird. Genauso sind die Streuungen der 
MeBwerte bei einigen Frequenzen dahingehend zu 
deuten, das dort zwei verschiedene Wellentypen 
_ etwa gleichstark auftreten. (An verschiedenen Orten 
der Schiene ergaben sich bei diesen Frequenzen 
: durch Uberlagerung verschiedene Wellenlangen.) 

“a 


2.2 Die réumliche Amplitudenverteilung auf der 
| schwingenden Schiene 


Die Schiene wurde an einem festen Ort, z. B. am 
Li Schweifsto$, durch den Korperschall-Sender sinus- 
_ formig angeregt und die resultierende raéumliche 
Be Verteilung der Schwingungsamplituden mit dem 
_ Kérperschall-Mikrophon abgetastet, das mit kon- 
_ stanter Geschwindigkeit auf der Schiene entlang- 
_ gefahren wurde. 

Die gefundenen Amplitudenverteilungen haben 
sehr unterschiedliche Gestalt, wie aus den in Abb. 4 
; De iciceehaieh Ausschnitten hervorgeht. Die Re- 
it _ produzierbarkeit wurde sorgfaltig gepriift; dabei er- 
gab sich, daB Verlauf und Aussehen der Amplitu- 
denverteilung sich sofort anderten, wenn die Anre- 
' gung an einem anderen, wenn auch benachbarten 
- Ort auf der Schiene erfolgte. 


Abb. 4. iele fiir die’ Verteilung der Schwingungs- 
plituden lings der Schiene; senkrechte 


am: 
Schwingungskomponente. 
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2.2.1. Die Abstande der Schwingungs- 


maxima 


Auf den erhaltenen Pegelschreiberstreifen wurden 
die Abstande der Maxima ausgemessen. Waren meh- 
rere Welligkeiten gleichzeitig vorhanden, so wurde 
jede von ihnen einzeln ausgezahlt. Die erhaltenen 
Werte sind in Abb. 5 dargestellt. Durch Vergleich 
mit den aus einer Phasenmessung gewonnenen Wel- 
lenlangen (Abb. 3) zeigt sich, daB bei Frequenzen 
zwischen 10 und 20kHz (bei den waagerechten 
Schwingungen bis 32 kHz) die Maxima einen Ab- 
stand von einer halben Wellenlange besitzen, also 
durch Reflexionen hervorgerufen werden. Der bei 
34 (bzw. 32) kHz beginnende steil ansteigende 
Kurvenzug lat sich durch eine Superposition zweier 
mit verschiedener Phasengeschwindigkeit fortlaufen- 
der Wellen (entsprechend den ersten beiden Disper- 
sionsasten in Abb. 3) erklaren. Die tibrigen noch 
aufgefundenen Abstiinde der Maxima sind keine 
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Abb. 5. Die Abstiinde der Schwingungsmaxima auf der 


Schiene; (a) senkrechte, (b) waagerechte 
Schwingungskomponente. 
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Abb. 6. Die Welligkeiten in der Schwingungsverteilung 
auf der Schiene; (a) senkrechte, (b) waage- 
rechte Schwingungskomponente. Die einander 
entsprechenden Kuryenstiicke in den Abb. 5 
und 6 sind in gleicher Weise bezeichnet. 
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ganzzahligen Vielfachen einer halben Wellenlange. 
Ihre Entstehung mu daher auf eine Uberlagerung 
von mehreren Schwingungsformen zuriickgefihrt 
werden. 


2.2.2. Die Welligkeit auf der schwin- 
genden Schiene 


Aus den gemessenen Amplitudenverteilungen 
lings der Schiene wurde ferner das Verhaltnis der 
Maxima zu den Minima, also die Welligkeit, er- 
mittelt, die quasi als Gewichtsfaktor bei der Beur- 
teilung der Abstande der Schwingungsmaxima ver- 
wendet werden mui. Diese Welligkeiten sind in 
Abb. 6 fiir die einzelnen Kurvenziige der Abb. 5 
dargestellt. 


2.2.3. Die Dampfung der Wellenaus- 
breitung langs der Schiene 


Die Dampfungen der beiden Schwingungskom- 
ponenten wurden ebenfalls den registrierten Am- 
plitudenverteilungen entnommen und sind in 
Abb. 7 wiedergegeben: In Ubereinstimmung mit 
den Untersuchungen von Naaxe, der bereits im 
Tonfrequenzbereich einen Abfall der Dampfung 
mit zunehmender Frequenz gefunden hatte, besitzt 
die Dampfung im Ultraschallbereich auferordent- 
lich niedrige Werte und bleibt im wesentlichen 


unterhalb 0,3 dB/m. 


Es wurde ferner die Dampfung ermittelt, die 
eine Welle beim Passieren einer Stofstelle erfahrt: 

Fir Schienen mit StoSliicken fand Naaxe seiner- 
zeit eine stark mit der Frequenz ansteigende 
Dampfung, die bei 2,35 kHz den Wert von 21 dB 
pro StoB besa. Die hier referierten Untersuchun- 
gen hingegen wurden an_ liickenlosen Gleisen 
durchgefiihrt. Dabei ergab sich fiir Abbrennstumpf- 
schweiBungen eine kaum meBbare Daimpfung; an 


—- - 


++ 


ees" + oe Se 


ie 60 80 kHz 700 
Abb. 7. Die Dampfung einer im Schienenkopf entlang- 
laufenden Welle; (a) senkrechte, (b) waage- 


rechte Schwingungskomponente. 
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StoBen, die durch Thermit-Schwei8stellen. verbun: 
den waren, betrug die Dampfung der durchgehen-— 
den Welle maximal 6 dB: 


2.3. Deutung der Schwingungseigenschaften der 
Schiene 


Zusammenfassend ergeben sich fiir die Schwin- 
gungseigenschaften der Schiene folgende Haupt- 
merkmale: . 

Die Dampfung der in der Schiene entlanglau- 
fenden Wellen ist im-.Ultraschallbereich auferst ge- 
ring und auch die Abbrennstumpfschweifungen 
rufen kaum einen Abfall der Amplitude hervor. 
Dadurch k6nnen sich irgendwelche in der Schiene, 
z. B. durch Schlage, angeregte Wellenziige ungehin- 
dert ausbreiten, zumindest bis zum nichsten Ther- 
mit-SchweiBstoB, an dem sie teilweise reflektiert 
werden. 

Wird die Schiene irgendwie zum Schwingen an- 
geregt, so zeigen sich grofe Welligkeiten in der 
Amplitudenverteilung langs der Schiene: Zu einem 
Teil stammen diese von Reflexionen an den Ther- 
mit-SchweiBstellen her, bei hoheren Frequenzen 
konnen auferdem in der Schiene gleichzeitig meh- 
rere Schwingungsformen auftreten, die sich zu orts- 
festen Schwebungen superponieren. Da die in der 
Schiene vorhandene Gesamtwelligkeit durch Uber- 
lagerung dieser Schwebungen mit den stehenden 
Wellen entsteht, muB die resultierende Amplituden- 
verteilung sehr stark vom Anregungsort der 
Schwingungen abhangen, eine Folgerung, die 
durch den experimentellen Befund voll bestatigt 
wurde (vgl. Abschnitt 2.2). 

Kine Korrelation zwischen der gemessenen Am- 
plitudenverteilung auf der Schiene und den auf- 
genommenen Beschleunigungsspektren bei der 
Uberfahrt eines Zuges tritt nicht auf und war auch 
nicht zu erwarten: Bei der Uberfahrt wird die 
Schiene an mehreren Stellen gleichzeitig, namlich 
durch die Rader, zum Schwingen angeregt, und 
erst die Uberlagerung samtlicher dabei erzeugter 
ortsfester Amplitudenverteilungen ergibt den tat- 
sachlich, durch den Zug verursachten Schwingungs- 
zustand der Schiene. 


Ich danke Herrn Prof. Dr. E. Meyrr fiir sein 
Interesse an diesen Untersuchungen und zahlreiche 
wertvolle Hinweise. 

Die Mittel zur Durchfiihrung der Untersuchun- 
gen wurden vom Bundesbahnzentralamt, Minden 
(Westf.), zur Verfiigung gestellt. 

(Eingegangen am 10. April 1956.) 
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Summary 

A detailed investigation is presented. of the electric input impedance of a piezoelectric 
bar and plate, loaded at one face with the acoustic impedance of an interferometer path, 
and free at the opposite face. The calculations start from the basic piezoelectric equations 
of state and lead rigorously to the strict expression for the total electric input impedance. 
The various properties of this impedance are discussed for variable path length of the 
interferometer, for variable acoustic wavelength and for a variation of the driving electric 
frequency. Under ordinary experimental conditions the input impedance can be repre- 
sented by means of circle diagrams. The theoretical results are illustrated by numerical 
examples for water at 15 Mc/s and air at 1 Mc/s. 

The theory presented provides a rigorous foundation for the well known theory of the 
acoustic interferometer given by Husparp and extends it in various directions. The treat- 
ment includes the basic equations needed for the investigation of highly absorptive media. 


Sommaire 

On présente une étude détaillée de V’impédance électrique d’entrée d’un barreau et 
dune plaque piézoélectrique, chargés d’un cété par Vimpédance acoustique d’un trajet 
interférométrique et libres de l’autre cété. Les calculs sont faits a partir des équations 
piézoélectriques fondamentales et conduisent 4 une expression rigoureuse de l’impédance 
électrique totale d’entrée. On examine les différentes propriétés de cette impédance en 
fonction de la longueur du trajet dans Vinterférométre, de la longueur d’onde acoustique 
et de la fréquence électrique d’excitation. On peut représenter l’impédance d’entrée par 
un diagramme circulaire, dans les conditions expérimentales habituelles. Les résultats 
théoriques sont illustrés par des exemples numériques relatifs 4 l’eau (15 MHz) et a lair 
(1 MHz). 

La théorie présentée constitue une base rigoureuse pour la théorie bien connue de 
Husparp relative 4 l’interférométre acoustique, et peut étre étendue dans différentes direc- 
tions. L’exposé concerne aussi les équations fondamentales nécessaires pour étudier les 
milieux trés absorbants. 


Zusammenfassung 

Die elektrische Eingangsimpedanz eines piezoelektrischen Generators (Stab und Platte) 
wird ausfiihrlich diskutiert; die eine Seite des Generators ist frei, wahrend die andere 
durch die akustische Impedanz einer Interferometerstrecke belastet ist. Ausgehend von den 
piezoelektrischen Grundgleichungen wird ein exakter Ausdruck fiir die gesamte elektrische 
Kingangsimpedanz hergeleitet. Ihr Verhalten wird fiir verainderlichen Reflektorabstand, 
verinderliche akustische Wellenlange und verianderliche Generatorfrequenz untersucht. 
Unter den iiblichen experimentellen Bedingungen la%t sich die Eingangsimpedanz durch 
Kreisdiagramme beschreiben. Die Ergebnisse werden an zwei Beispielen fiir Wasser bei 
15 MHz und Luft bei 1 MHz erlautert. 

Die yorliegende Theorie gibt der bekannten Theorie des akustischen Interferometers 
nach Huszarp eine strenge Grundlage und erweitert diese in verschiedenen Richtungen. Sie 
enthalt auch die Ausgangsgleichungen zur Untersuchung von Medien hoher Absorption. 


PART I 
GENERAL EXPRESSIONS FOR THE ELECTRIC INPUT IMPEDANCE 
OF A PIEZOELECTRIC BAR AND PLATE, LOADED AT ONE SURFACE 


surface is parallel to the radiating face. The width 


z 1. Introductory remarks of the acoustic beam is assumed so large in compa- 


__ The acoustic interferometer is regarded as consist- * This paper is essentially based on work done 
ing of a piezoelectric crystal which is free at one face under a contract between the Office of Naval Research, 
and radiates from the opposite face into the inter- Washington, D.C., and the California Institute of Tech- 


ie : : : nology, Pasadena, Cal. (Techn. Report Nr. 3, Janu- 
ferometer path containing the fluid or solid medium Sees etG23) "Ths publicetion ofthe paper wes radde 


whose properties are to be studied (Fig. 1). The possible by support of the Division of Engineering and 
interferometer path ends at a plane reflector, whose Applied Science, Harvard University. 
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rison with the acoustic wavelength / that the wave 
front can be treated as plane. The piezoelectric cry- 
stal is driven by an alternating voltage V of angu- 
lar frequency wo =2af, the current to the crystal 
being 7. The response of the interferometer to the 
driving voltage is completely characterized by the 
electric input impedance Z. Indeed, if Z is known, 
the response of the device in any driving electric 
circuit can be calculated. 


Recs 


V M 


R 
ir) 


Fig. 1. Schematic outline of the acoustic interferometer. 
C=piezoelectric crystal of length 1, M=me- 
dium filling the path of length d, R=plane re- 
flector, V,Z=voltage and current at crystal 
input. 


l d 


The reaction of the acoustic wave-motion in the 
medium upon the crystal can be expressed in 
terms of the acoustic input impedance Z,, of the 
interferometer path as measured at the boundary 
between crystal and medium. Our first step there- 
fore will be to calculate the electric input impedance 
Z of a piezoelectric crystal loaded at its radiating 
surface by an acoustic impedance Z,,. The opposite 
surface of the crystal is assumed to be free (non- 
radiating). 

The crystal is assumed to have the shape of either 
a thin bar, whose breadth and thickness are small 
in comparison with the length, or a plate, whose 
linear surface dimensions are large in comparison 
with the thickness. The bar is regarded as driven in 
compressional lengthwise vibration, the plate in 
compressional thickness vibration. 

These assumptions are those commonly made in 
dealing theoretically with piezoelectric bars or pla- 
tes. The waves in the interferometer path are treat- 
ed as plane, which implies that the dimensions of 
the radiating cross-section are large in comparison 
with the wavelength. In the case of a plate this con- 
dition is well met. In the case of a bar, however, 
only vibrations of high-order harmonics would meet 
the two conditions that the lateral dimensions are 
small in comparison with the length of the bar, and 
large in comparison with the wavelength. If the bar 


x4) Wy shes 


Fig. 2. Piezoelectric bar of length 1, breadth w and 
thickness t, driven by a uniform electric field E,. 
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is driven in its fundamental mode, the admissibility 
of assuming plane waves will depend upon experi- 
mental conditions and the nature of the desired re- 
sults. 


2. Electric input impedance of a piezoelectric bar 
in lengthwise vibration with an acoustic load 
at one face 


Both the breadth w and the thickness ¢ of the bar 
(Fig. 2) are supposed to be small in comparison 
with the length 1. The ratio w/t is assumed to be 
large enough to ensure a practically uniform electric 
field E,, within the bar. The electrodes are assumed 
to cover the entire lateral faces of area d=wl. All 
surfaces except that at y=I, which is adjacent to 
the interferometer path, are free. Consequently all 
stresses except T,,=T. can be assumed to vanish 
everywhere ([15], p.61); 7, vanishes at y=O. 
Only displacements parallel to y , and the compres- 
sional strain Sy, need be considered. 

The appropriate equations of state for a quartz 
X-cut bar with length parallel to a y-axis are! 


Sos 35a Fo dig Bae (1.1) 
Dp=d oT +-e Eve (1.2) 


For the present purpose we write gq” =1/st ; 
H =d,»/s%, effective piezoelectric stress-constant; 
V =E, tis the driving voltage, assumed to be sinusoi- 
dal in time. F; is uniform throughout the bar, but 
the electric displacement D,, while uniform in the 
x-direction, varies in the y-direction owing to T,(y). 

Equations (1.1) and (1.2) may now be written 
as follows, with subscripts omitted, for any point y 
along the length of the bar: 


T=q®S—HD|t, (1.3) 
D=HT/qe +27 Vt. (1.4) 


Since these equations are linear and V is sinusoi- 
dal, the instantaneous values of 7, S, and D are 
sinusoidal in time as well as being functions of y . 
In what follows, time-dependence will be represent- 
ed by the factor e}* , which itself is suppressed. 

The equation of motion, in Lagrangian coordi- 
nates, for a plane compressional wave with displace- 
ment 7 in the y-direction, is 

3° _ OT (y) 
eae = a (1.5) 

or, since 0/Ot=jo, 
—w? o,.n =9T/dy. (1.6) 


1 See [6], p.67 and [7], p.580; also [2]. For 
other crystals than quartz suitable changes in sub- 
scripts can be made ([6] and [7]). Numerical values 
for quartz are given in Appendix A. 
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From eqs. (1.6) and (1.3), 


c) ( De rELe Vs 
Re a 1 ye es SS iy HG) — — 
ay o Oa, 7105.5 0) 0c dé a: ‘) 


OT. wo, wo, HV 

pee eb "Og LY 7 

or By? ie g ey (1.7) 
Introducing y.? = w? o,/g¢® , we get 
oT “ HV 
rt Ta —y2th 

aye ts ieee (1.8) 


In general 7, is complex. The losses in the crystal 
are included in the complex stiffness g”. They are 
treated in the same way as the losses due to visco- 
sity in fluids. 


Setting gq =qotjoF, (1.9) 
where q¢o is real, we have 
2 
2 WO” Oc 
er aa he 1.10 
ED (1.10) 


For smal] losses (w F < q.9), we may write 


pence] / 2 (1- ee) 
V0 2 deo 


or, introducing the crystal phase velocity c,.= 


V4e0/ 0c and, as usual, a space-attenuation constant 
a, we have 
w ; 21 
Vor jo=- 
Co he 


jo. (1.11) 


Although the crystal losses are treated as small, 
still eq. (1.11) can be regarded as correct for large 
losses also; but in this case the phase velocity c, 
and the acoustic wavelength 4, in the crystal are 
functions of and F/q,.9, as is found by taking the 
square root of eq. (1.10). 

Further losses can arise if imaginary components 
of H and ¢? are taken into account. Within the 
range of frequencies observed so far, no such ima- 
ginary component of H — due for instance to re- 
laxation effects — seems to be indicated and we 
therefore regard H as real. A complexe «’ would 
imply dielectric losses in the crystal. For quartz, at 
least, these losses are extremely small. We shall 
treat &? as real; the imaginary component could be 
introduced if necessary. 

The solution of eq. (1.8) for 7 can be written 


T=asin(y,y+%) —HV/t. (1.12) 


The constants a and # are determined by the boun- 
dary conditions. At the free surface y =0 the stress 
T vanishes. From (1.12) we have therefore 


asind=HV/t. (1.13) 


At the opposite surface y=/ the crystal is in con- 
tact with the medium in the interferometer path. 
The reaction of the medium upon the crystal can be 
characterized by the acoustic input impedance co 


of the interferometer path, which is defined in the 
usual way by 


Zm = p(l)/u(l) r) (1.14) 


p(l) representing the excess pressure and w(l) the 
particle-velocity at the boundary between crystal 
and medium. Both pressure and _particle-velocity 
must be continuous at this boundary. Since a stress 
is usually called positive when extensional, whereas 
a pressure is called positive when compressional, we 
write p(/) = —T (1). Furthermore, u = 07/dt=j 7. 
Thus the boundary condition at y =/ can be written 


T (1) /n (0) =—joZy. (1.15) 


The expression for Z,, will be computed in Part III, 
Section 1. From eqs. (1.6) and (1.12) 7 is found 
to be 


n= = BSE cos(yey +), 
ved 
and, introducing ~=7.1, we obtain by eqs. (1.12) 
and (1.15) 
gein(e +0) Vit = jo 


Tks = 
acos(y +9) 764 


E 
It is convenient to introduce here the dimensionless 
quantity 

Zm = OL p/¥e Gg? . (1.16) 


Thus we obtain 


a) ee 7 97 cos (GEO). (1.17) 


By developing sin(g+”) and cos(~+¥#) in terms 
of ™ » and 8" # and replacing asin ® by means 


of eq. (1.13), we find by a short calculation 


acos ia LY 1= 008 = j 2m sin P 


: : (1.18) 
t sIn P —J Zp Cos /P 
In order to establish a relation between the cur- 
rent J and the voltage V, which will lead to the elec- 
tric input admittance Y, we integrate eq. (1.4) 
with respect to y : 
l l 


[oa 4 [rayee™ (1.19) 
q’ 
0 


The displacement-current 7 through the crystal is 
given by 
l 


l 
taw{ Pay=jow [ Day. (1.20) 
0 6 


From (1.12) we find 


I 
[Tay- 4 [cos #(1 — cos) + (1.21) 
“ Ve 
0 


+ sin # sin py] — 


FAVE 
ae 


+ % 
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Replacing asin? and acos® by egs. (1.13) and 


(1.18) and inserting (1.20) and (1.21) into (1. ye 


we find finally, with y,=9/l, 
2 
Tajov el 7) 4 
ENG epee 
a 2 ee ca 
@ git. <-Ljzpeotp 


The input admittance Y of the crystal is defined by 
I=YV and is found immediately from this equa- 
tion to be 

2 A H? tan p/2 — j 2n/2 
pgt 1—jzncoty 


Meri: 


Y=jwC,+jo° (1.22) 


(1.23) 


her Sai Mad Ne 
where apes 
The admittance Y can be represented by the electric 
circuit shown in Fig. 3. C, is the parallel capaci- 
tance of the crystal, with an effective permittivity 


ceellee 
\ é qd 


(1.23 a) 


Fig. 3. Equivalent electric circuit for the loaded piezo- 
electric crystal. 


The impedance of the branch in parallel with C, 
is given by 


. pgit 1-jzncotge 

32m AH? tan y/2—j Zm/2° 
Kq. (1.24) is the rigorous expression for the electric 
input impedance (or admittance) of a piezoelectric 
bar, loaded at one end with an acoustic impedance 
Z., (expressed by z,,). By including imaginary parts 
of the quantities p, g”, «7, and H in eqs. (1.22), 
(1.23) and (1.24), all conceivable losses in the 
crystal can be taken into account. 


(1.24) 


3. Electric input impedance of a piezoelectric 
plate in thickness vibration with an acoustic load 
at one face 


The piezoelectric plate (Fig. 4) is assumed to be 
infinitely extended in the y- and z-directions. The 
acoustic wave motion is treated as a plane problem 
(9/Oy=9/3z=0), the only displacement £ being 
in the x-direction and the only strain being S,. In 
a finite plate other strain-components,. together 
with all six stress-components, may be present, de- 


Fi: > er ee ee ey » 2 Se ie 7 Py A en be a? ee 
ve MA Diet to > nt Vor VCR Lat Sarees re 
/ = 4 ‘ > . 
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oy 
z 


pending on the piezoelectric and elastic coefficients — 


of the particular crystal; they are disregarded in the 
present treatment. It will therefore be assumed that 
associated with the strain S,, only the single stress 
T, need be considered. 


Z4(3) 


(1) 


Fig. 4. Piezoelectric plate of thickness /, driven by an 
electric field Ez. 


The equations of state for a quartz X-cut plate 
then become ([7], p. 580; [2]) 


Ty =ctyS, E11 E,, (1.25) 
Dy =e), 53+ er ky. (1.26) 


Owing to the polarization produced by the strain 
S,(x), the electric field E, is a function of x over 
the thickness J of the plate. The electric displace- 
ment D, , however, is independent of 2, as no free 
space charge exists within the plate. We write eqs. 
(1.25) and (1.26) in the simplified form 


T =q® S—HE, 
D=HS+éE. 


(1.27) 
(1.28) 


These equations for the plate are similar to eqs. 
(1.3) and (1.4) for the bar. The numerical values 
of the constants, however, are in general not the 
same as for the bar (see Appendix A). 

Proceeding in the same way as in the case of the 


bar, we shall establish the differential equation for — 


T . It is desirable to use the electrical parameter D , 
which is independent of 2, instead of EF. By elimi- 
nating E between eqs. (1.27) and (1.28) one finds 


H? HH 
B= Ge Te Sea 
q ( 5 a) 7s (1.29) 
Introducing 
Ke 
=a (14 = (1.30) 
we have 
Tay SD, (1.31) 
é 


By a procedure analogous to that used in deriving 
eq. (1.8), we obtain 


eT 
O22 


+72T= yo 


3D, (1.32) 
é 


where, in the case of the plate, y.2=w?9,/q?.— 


et ys a + eo > Pars io 
abe le is 
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Applying the same boundary conditions? as in eq. 
~ (1.15) at x=/ and the condition T=0 at x=0, we 
obtain the solution for 7 corresponding to eq. (1.32) : 


T = asin(y.2+0) —" p. (1.33) 
é 


The two equations determining the constants of inte- 
_ gration a and @ are found from the boundary con- 
ditions; they are 


asin =" p (1.34) 
é 
and 
: H : 
asin(p +0) = —,D+jaz%mcos(pt+), (1.35) 
€ 


where 
Ly, 
Ve ad 


and g=y,l. (1.36) 


zn = 


From eq. (1.35) it follows by means of eq. (1.34) 
that 


acos0= 


HD 1 = cos Gaia sin Pp (1.37) 
é> sing —j Zn cos P 

On eliminating S between eqs. (1.27) and (1.28) 

we find 


% e gh q 
T= —(H+ 7 B+ i (1.38) 
and integrating eq. (1.38) with respect to x, we 


have 


I 
S 7 EDI 
Phen fr eles : 
[ra ( tae) + (1.39) 


i) 
On the other side we find from eq. (1.33) 


U 1 
[rac [acos #(1—cos@) + 
“ Ve 


(1.40) 
+ asin? sin yp] — ae 
@ 


Inserting asin ® and acos® from eqs. (1.34) and 
(1.37) we obtain after a short calculation 
U 
[raz- EM tan p/2 — j 2m/2 Bey (1.41) 
Vee 1—j 2m cot Pp € 


0 


Introducing this result into eq. (1.39), we get 
S gE S gE 
(w+ a Pr DE (H+ -f}\. 


_ 24D tan y/2~j2n/2 
j,6 /1—jz,cotp 


2 The metallic coating on the plate is assumed to be 
so thin that its acoustic impedance is negligible. 
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Now, by eq. (1.30) 
& ee pS EP 
nit oars 2 (eM) 88, 
and we find, introducing the current /=jwD A, 
were a! 2H? tan g/2—j %m/2 
joAe jpwA(e)2Q? 1—jzp cot 
ShZ, 
where Z, the input impedance of the crystal is 
eae 2H? tan y/2—j zm/2 


Yoo es een 
D oo eras 


with Cy=A él, (1.43) 


According to eq. (1.42) the impedance Z can be 
represented by the equivalent circuit of Fig. 5, 
where 
ye j 2H? tangy/2—j2,/2 
oCr peSq? 1—-jzecoty * 


(1.44) 


Eq. (1.42) gives rigorously the electric input im- 
pedance of a plate in thickness vibration; as in the 
case of the bar, all conceivable losses are included. 
if p , q?, H, and &* are regarded as complex quan- 
tities. 

A 


areal ea 
q Z 


Ow 


Fig. 5. Equivalent electric series circuit for the loaded 
piezoelectric plate. 


Although Fig. 5 corresponds to the simplest ana- 
lytical expression for the impedance of the crystal 
plate, as given by eq. (1.42), it is of advantage in 
many cases to transform the circuit of Fig. 5 into 
the more usual form shown in Fig. 3. It can easily 
be shown that the two circuits are equivalent in that 
they represent at all frequencies the same input im- 
pedance or admittance with respect to the points A 
and B, if the following relations exist: 

2,1 iCieC hs... a 
Inserting Z; from eq. (1.44) and C;=C, from eq. 
(1.43) into eq. (1.45) we obtain after a short cal- 
culation 


PL 
Teeth Boa 
1 20 AP Sent 
H? tang/2_. (: 2 He.) 
se n/[2) (2 cote — 2 — 
Br [2 30!) PoP 9 sap) 


jtan 9/2+2z,/2 


As already stated, the expression for Z, is more 
complicated than that for Z;. For the crystals com- 
monly used, however, the quantity |H?/eSq?| is 
small in comparison with unity; this fact permits 
a simplification of eq. (1.46). 
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PART I 
ELECTRIC INPUT IMPEDANCE OF THE FREE PIEZOELECTRIC 
BAR AND PLATE 


1. The free bar 


The formulas for a free crystal, that is in ab- 
sence of any external load, follow immediately from 
the previous general formulas. If the surface at 
y =l is free, the stress 7(J) vanishes. Therefore in 
eq. (1.15), the condition that T (1) =0 is equivalent 
to Z,=0 and z,,=0 according to eq. (1.16) for a 
bar and (1.36) for a plate. Thus we obtain the im- 
pedance Z, for the free bar, corresponding to the 
upper branch of the equivalent circuit in Fig. 3, by 
setting z,=0 in eq. (1.24): 

Zig= — C2a1) 
Usually this impedance is represented in the neigh- 
bourhood of the series resonance @, by the circuit 
shown in Fig. 6. Series resonance occurs when the 
imaginary part of Z, vanishes, or, in terms of C 
and L, when w,2CL=1. 

With 


ol 
P= yel= a 


c 


—jal 


and regarding the losses in the crystal as small so 
that al <1, we have 


we cot P 24 al apes al cot (2 Ele al 
2 Deh? Ce aD s \2 Ce i5) 
or, retaining only first-order terms in al/2 and de- 
veloping the cot-term, 
DE a Oe OR LS ON Wie O21 
9 Cot 5 ee wr: Bas —j@ | co Ts 
w 1/2 e 
eae eer S. 
sin? (« 1/2 ¢,) 

The imaginary part of Z, therefore vanishes when 
cot(wl/2 ec.) =0 
or 0,/2¢.= 


(2.2) 


(2n+41) x/2, 


where n=O, 1, 2, 3,:°-. Thus, as is well known, 
series resonance for the bar occurs at 


(2.3) 


@,= (2n+1) xe,/l. (2.4) 


The real part of Z,;, representing the resistance R), 
of the bar at resonance, follows from egs. (2.2), 
(2.3), and (2.1) 3 


iy Re a 


"4A BP cy ae: 


3 The consideration of small imaginary parts of q” 
or H would only lead to corrections in second order 
terms and can therefore be omitted. 


As seen from eq. (2.2) the real part of Z, is a func- 
tion of w; this means that R, varies with Aw = 
@—q@,. But as al is assumed to be small, terms of 
order al: Aw can be ignored near resonance, and 
R, can be regarded as a constant as is commonly 
the case. 

The quantities Zj,and C), for the bar in Fig. 6 are 
found from eqs. (2.1) and (2.3) by using the rela- 


tions * 
| dZ, | 1 
Ly= . aes he a and Cy= oe ie (2.6) 
After short calculations one finds easily 
EyR 
peer ee 
SANTA c= CDR) 
8 AH? 1 
and = x? g® t (2n41)”’ 


in agreement with the well-known values of these 
quantities ([5], p. 297). 

Antiresonance is here defined as occurring 
when the imaginary part of the total admittance Y 
becomes zero, that is, when Im(j @ C,+1/Z,) =0 
(see Fig. 3). Disregarding terms of order higher 


than a/, the angular frequency @, for antireso- 
nance follows from eqs. (2.1), (2.2), and (1.23): 
W,l Wal H 
ac oie = lg a Fr 4 (2.8) 


Assuming that H?/e? g® <1, and moreover that 
w, l/2 c, is very close to (2n+1) x/2, we can write 
w,l 


a5 = (2n+1) 5 +6 


: (H?/e7 g® €1). 
DiGa 


Then after a short calculation it is found from eq. 
(2.8) that 
2H? 
~ (2n4+1) net g? 
From eqs. (2.3), (2.9), and (2.10) it is seen that 


w D) Op sist 4 os 
a— Ws aCe 1 vas Bg’ (2n+1)? 7? 


<1 (H2/e? g® <1). (2.9) 


(2.10) 


This is the difference in angular frequency between 
antiresonance and series resonance for a free bar. 
For the usual crystals it is a very small quantity 
(see Appendix A.1). 


4 Since Z, is represented in the neighbourhood of re- 
sonance by Z,=jw Ly) +1/j o Ch + Rp, dZ,/dw =j (Ly + 
+1/w? Cy), and since 2? Ly Cyp=1, | dZ,/d@ |n—ox 
=2 Ty. 
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In the foregoing treatment, as well as in what fol- 
lows, the effect of the electric circuit (matching net- 
work, generator, etc.) is disregarded. So far as the 
use of a crystal in the interferometer is concerned, 
the theory shows that only crystal-resonance is of 
interest. Therefore no heed will be given to the cha- 
racteristics of the external electric circuit. 


2. The free plate 


The impedance Z, corresponding to Fig. 3 for 


this case is obtained from eq. (1.46) if we set 


4 

4 EN ee 
y > (2n+1) 0 &q? 
a 


2m =0. We find 
D 2 
aigaty (- G8 TP 
Zig = Ta dig: Fy? (Z cot Pots ap : (2.11) 


The term H?/e5 g? may be attributed to the effect of 
polarization, which makes the electric field E, vary 
along the z-direction (Fig. 4). Such an effect does 
not occur in the bar, where the electric field is per- 


' pendicular to the direction of mechanical vibration. 
_ As already mentioned, H#/e*q? is a small quantity 


for the crystals commonly in use. When this term is 
neglected, the expression for Z, in eq. (2.11) has 
the same form as the expression for the free bar in 


meq. (2.1). 


Series resonance for the free plate occurs, 
according to eq. (2.11) when® 


D FH? 
Re E (g mre >)|=9. (2.12) 


For small losses (al <1) we make use of the same 
development of (~/2) cot(q/2) that led to eq. (2.2). 


_ Assuming that the imaginary parts of g?/H? and of 


eS are small quantities, we find the condition for 
series resonance to be 


wgl wgl ( I? 
= —}|, De Aes 
20. “ane 2 Co fie eg? ( ) 


_ A simple approximate solution of this equation fol- 


lows under the assumption that Re(H?/s5 q?)<1, 


"by setting w, 1/2 c. = (2n+1) (x/2) —6,, where 0, 

is a small quantity. Retaining only first order terms 

in 6, and ignoring an imaginary part of H?/e*q”, 
we find from eq. (2.13): 


. Li (2.14) 


5 Eq. (2.12) is identical with the result obtained by 
a different approach by R. Becumayn in his theory of 


the vibrating plate ([1], eq. 41). Becumann’s corrective 
term in g? (c’ —c in his notation, eq. (32 a)), however, 
is too great by a factor 2, due to the omission of a fac- 
tor 1/2 in his expression for the piezoelectric inter- 


action-potential y (in Brcumann’s notation); see also 


R. Becrmann [2 a]. 
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Thus series resonance for the plate occurs at® 
Tl Co Peabo) SH% 
l (2n41)-al = 9? © (2:19) 
(HP eto? <1); 


Ms = (2n+1) 


The real part of Z,;, representing the loss-resistance 
R, (Fig. 6) for the plate, is found in the same way 
as for the bar. From eqs. (2.11), (2.12), (2.14), 
and using the development eq. (2.2) for small los- 
ses, we obtain by a short calculation, assuming 


al<1 and 6,<1, 


Di 
ji els 2 oy re Able oe i 
(Rp) o= 0, aE A Fach (2.16) 
° |} 8 
a L R 


Fig. 6. Representation of Z, in Fig. 3 near resonance by 
C, L and R in series. 


The quantities L, and C, for the plate, corres- 
ponding to the equivalent network of Fig. 6, are 
obtained by using eq. (2.6). Neglecting terms in 
dp, we have 


D 73 
qg?l 
jamie Mea) 
8AM ce2 
: (2.17) 
8 A H? 1 
and — ae 


cheat pa rt WAT Oe et ae 
in agreement with the well known values of these 
quantities ([5], p. 323). 

Antiresonance of the free plate occurs when 
Im(1/Z,+j@C,)=0. From egs. (2.11), (1.43), 
and (1.45) we obtain 


Re ( cot 2) =0 


and for small losses, replacing t cot ? by eq. (2.2), 


9» 
we find 


a“ 


w, W2e,= (2n+1) x/2. (2.18) 


From eqs. (2.15) and (2.18) it is seen that for a 
free plate the difference in angular frequency be- 
tween series resonance and antiresonance is 
(2c,/l) dp, where 6, is given by eq. (2.14) when 
H?/e5 g? <1. Thus? 

I? 4 Oa 


= = CO Je5 gq? € 1). 
S yD 28 
eg (2n+1)*ax (2.19) 


W,—-Ws= 


6 It is easily shown that eq. (2.15) agrees with 
Brcumann’s expression for series resonance, as given in 
[1], also in [5], p. 314. 

7 Eq. (2.19) agrees with Brcumann’s results, [1], 
eqs. (61) and (61a). 
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PART IJ 
MECHANICAL AND ELECTRICAL IMPEDANCE 
OF THE INTERFEROMETER PATH 


1. Mechanical impedance of the interferometer 
path 


The mechanical input impedance of the interfero- 
meter path at the radiating boundary of the crystal 
is by definition, eq. (1.14), 

Td) pig), 


(Zn) =a — aay =— A 


if we choose the coordinate system of the interfero- 
meter path as indicated in Fig. 7. T, p, and u denote 
stress, excess pressure, and particle velocity of the 
medium in the path. R is a plane reflector parallel 
to the crystal boundary. The differential equation 


{- ——|.., 


pont S10) 


(3.1) 


Fig. 7. Simplified sketch of interferometer with crystal 
C at x= —d and reflector R atx=0. 


for a plane compressional wave motion, in Langran- 
gian coordinates (displacement £), is 
O2é Sy ets Op o2é 
SS CS . = ot ee = 
Cn 92 Cuore Oz 2x2 
where @,, is the density and g,, the stiffness of the 
medium. To include the losses we assume gy, to be 
complex: 


(322) 


In = Im0 + J OF: (3.3) 
The solution of eq. (3.2) is 
F(z) =E(e"™ + Term), (3.4) 
where Ym? = (2 —j an) = 2 . (3.5) 
Cn m 


The phase velocity of the plane waves according 
to the solution eq. (3.4) is, by definition, ¢, = 
«w/Re(ym). If the losses in the medium are assumed 
to be small, c,, can be expressed by 


Se Gmo/Om (w Fy < qmo) : (3.6) 


The quantity J’ in eq. (3.4) is the reflection coeffi- 
cient of the reflector R. The particle velocity 
u—Oé/3t=j w E at the reflector (x =0) is given by 
eq. (3.4): 


(u)2-9=joé (1+). (3.7) 


If we construe the two terms in eq. (3.4) as belong- 
ing to two wave components progressing in the posi- 
tive and negative x-directions, the ratio of the par- 
ticle velocity u, of the positive wave to u_ of the 


negative wave at the reflector is given, according to 
eq. (3.7), by 


(2=) =| leas (3.8). 


u, 
In general I” is complex; the modulus | J’ | indicates | 
the diminution of the reflected wave-amplitude with 
respect to the incident amplitude, and is associated 
with the energy absorption in the reflector. The 
phase angle © indicates the “phase-jump” which 
the reflected wave undergoes with respect to the in- 
cident wave. For a perfectly rigid reflector (w);_9 
is zero by definition, and therefore J/’= —1 (or 
|I’|=1 and O=x) 8. 
The stress T = qm 0/9 is, according to eq. (3.4), 


T (x) = =j Ym Ym & (en Te eae (3.9) 


Thus the mechanical input impedance Z,,, as defin- 
ed in eq. (3.1), is found by using eqs. (3.4) and 
(3.9) atz=—d: 


TP Ym 9m 1—-T emt 


Z i eae on (3.10) 
"jor © 14 Pe ume 
Introducing the dimensionless quantity } 
ra Tere (3.11) 
we have, at x= —d, 
ie tat AG 
pa ae (3.12) 


The dimensionless quantity z,, in eq. (1.16) for the 

bar and eq. (1.36) for the plate becomes 

Pes a ee Ym Qm 1—r 
Ye Ve ¥eQc Ltr 


where g.=q" is the stiffness of the bar, and g,=q? 
is that of the plate. We introduce 


(3.13) 


=m 


Ym Yin 
Ye de 


(3.14) 


m=My+jm,= 


8 The suggestion is sometimes found in the litera- 
ture that some kind of “reflection-losses’” may occur 
also at the boundary between crystal and medium and 
that such losses might affect the value of J’. There is 
no room, however, for such considerations within the 
scope of the theory given here. The boundary condi- 
tions at the interface between crystal and medium re- 
quire continuity of pressure and displacement in the 
two adjacent media. Nothing else has to be introduced. 
The situation would be different if dissipation by heat 
transfer or surface layers due to impurities were taken 
into account; but such effects are beyond the scope of 
the theory treated here. 
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and have 2rq = 7. 


(3.15) 


From eq. (3.5) and the expression preceding eq. 
(1.8), the quantity m can be expressed by 


i tata -|/2adm m= |/2etee +j © Fin/qmo) 
Ve Ve Gc Ve Qc Geo (1 + j © Fe/qeo) 
@m Cm es 
Oc Ce 1+j@ F,/qeo 
(3.16) 


or m=my+jm= 


If the losses in both crystal and medium are small, 


we can write 
Om Cm ji+ 22 (2 re 


m—=mg +jm,= 
¢ Co 2 m0 =I 


where 

Om Em _ 

Oc Ce (@ Fn Gmo) 

om (7= ), (w F. < qeo) 
2 Gm0 = Ye 


The quantity m, can also be expressed by introduc- 


Mo = 


(3.17) 


ee 


ing the space-attenuation parameter a, as defined 
by eqs. (1.11) and (3.5): 
it — ve (CnC Cn a). (3.17 a) 
a) 


2. Electrical impedance associated with the 
interferometer path 


The impedance Z,;, as given in the general form 
of eqs. (1.24) and (1.46) for the bar and the plate, 
can be thought of as containing two parts: one be- 
longing to the driving crystal, the other to the load, 
that is to the interferometer path, as represented in 
our formulas by the quantity z,,. It is possible, 
indeed, to split the expression for Z, into the sum 
of two such parts. It can be verified easily that Z,, 
according to eqs. (1.24) and (1.46), can be written 
as follows: 


a) Bar: 
gt’ 19 P 
Oi: a ae 78) 
E 
Gt 9/2 
rr a ae , 
7. cotn 
2m 
b) Plate 
D 2 
Bs Ce 5 
Zi o ATP (Z cot 5 Sg? ~ (3.19) 
: tie 1B 
«w AH? 2 Pp 
aS Tasty 2 
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The first terms in eqs. (3.18) and (3.19) corres- 
pond exactly to the expressions for Z, found for the 
free bar in eq. (2.1) and for the free plate in eq. 
(2.11). The second terms may be construed as re- 
presenting that part of the total electric impedance 
Z, which is linked with the interferometer path. 

Writing Z, as the sum of the input impedance of 
the free crystal Z,, and the impedance of the inter- 
ferometer path Z;;, that is Z;=Z,.+Z4;, and intro- 
ducing y=y.l we find 


Lii=% 5 Be (3.20) 
pee 
where 
DO ag arr ae 
Cees te Oh aaa (8-21) 
and 
Be ig ieee ree Ee 
(Zo) plate = Supp Ae = Ve Oc 2A aie (ora) 


The equivalent circuit of Fig. 3 for the loaded 
erystal can be replaced by that shown in Fig. 8; 
Z,, is the impedance of the free crystal, which may 
be represented in the neighbourhood of resonance by 


the RLC-branch in Fig. 6, with the values of R, L 
and C as given in Part II. 


C 


Fig. 8. Equivalent electric circuit for the loaded piezo- 
electric crystal. Zi, impedance of the free crys- 
tal, Z1; impedance of the interferometer path. 


Our main interest lies in the impedance Z1;, as 
represented by eq. (3.20). This impedance is given 
by the same expression for the bar as for the plate, 
but z) might have different values in the two cases. 
Inserting z,, from eq. (3.15), we obtain 


(3.23) 


3. Circle diagram for the electric impedance Z;; 
of the interferometer path 


The usual purpose of the acoustic interferometer 
is the determination of the velocity c,, of sound and 
the absorption coefficient a,, of the medium filling 
the interferometer path. This is accomplished by 
changing the distance d or the acoustic wavelength 
A in the medium and observing, for example, the 
corresponding changes of the current / to the crys- 
tal (Fig. 1) at constant voltage V. All changes of 
parameters of the matching electric circuit connected 
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with the crystal are associated with the changes of 
the electric input impedance Z;=Z;,.+2Z;; of the 
loaded crystal (Fig.8). The behavior of Z;, be- 
longing to the free crystal has been treated in 
Part II. In all cases except when the frequency is 
varied, Zj, is constant, because it is independent of 
the variables of the interferometer path. Our main 
interest lies in the variation of Z;;, associated with 
the interferometer path. 

Eq. (3.23) gives Z,; as a function of r and 
gy =y_,l. First we note that, according to eq. (3.5), 
Ym can be written for any medium as 


jms Nar —j Am « (3.24) 
Cm 
Then, according to eqs. (3.11) and (3.8), 
es Pa aoe eee (3.25) 
or 
ae | ef (O-40dlhgy) : eo ztmt = Re-i, (3.26) 
where 
R=|I'|-0 =? 
ane ol 884 oe 2 ene 


ie Cm 


The variation of r is effected either by a) changing 
the pathlength d at constant acoustic wavelength 2, 
(variable path interferometer), or by b) changing 
Am at constant d (fixed path interferometer). The 
reflection coefficient |J'|ei® is treated here as a 
constant, independent of the frequency. 

Case a): 

The variation of d usually is of the order of magni- 
tude of a few half-wavelengths /,,/2 . Each change in 
d by 1/2 changes the term exp(j°4ad/A,) in eq. 
(3.26) by 2x and causes a periodic change in r by 
a full period. The term exp(—2a,,d) in eq. (3.26) 
decreases continuously with increasing d. If a, is 
not unusually large, this term can be regarded as 
“slowly” varying in comparison with the “rapidly” 
varying preceding factor in eq. (3.26). 

Case b): 

The term exp(—22a,,d) in eq. (3.26) is cons- 
tant; the periodic changes in r are effected only by 
an alteration of the number of half-wavelengths 
within the interferometer path, either by a change 
in the sound velocity c,, (due to changes in pres- 
sure, temperature or composition of the medium) 
or by a change of the driving frequency f. 

The first step in the graphical representation is 
to consider the locus of r. In the complex r-plane 
r=Re® describes exactly a circle with centre 
at the origin, when /,, is changed steadily (Fig. 9a). 
A change of d, however, causes r to describe a 
spiral, due to the factor exp(—2a,d) in eq. 


F.E.BORGNIS: ACOUSTIC INTERFEROMETER 


ACUSTICA 
Vol. 7 (1957) 


Fig.9. Locus in the r-plane of the quantity r= Re)? 
in eq. (3.26) with variation of (a) the acoustic 

» wavelength Am and (b) the path length d. In 

(b) the ratio between two successive radii is 


given by Rn+i=Rne “mm. 


(3.26). In what follows we assume a 
to be a relatively small quantity; we 
may then replace the continuous spiral by a se- 
quence of circles whose diameters decrease by a 
factor exp(—a,4) each time when r has undergone 
a change of one full period (Fig.9b). Only for 
media with unusually large absorption will the spi- 
ral contract so quickly that this approximation has 
to be replaced by a more exact treatment. This can 
be done by using the general expressions of eqs. 
(3.23) and (3.26). 

We next find the locus of Z;; in eq. (3.23) when 
r describes a circle around the origin, as indicated 
in Fig. 9, Z,; in eq. (3.23) can be written in the 
following form: 


Zhi Le I m/2 
<1 = = SS eee 
Zo (1+p)? l=py cee 
re 
l+p 
where P=i5 cot i (3.28 a) 


The term q/2 will be regarded as a constant, as is 
the case when the driving frequency is kept cons- 
tant. Even if the interferometer is used at varying 
frequency, d/A, is usually a large number and 
cot y/2 is a very slowly changing function near crys- 
tal resonance (#2) in comparison with the rapid 
change of r so that, as a first approximation, cot p/2 
can be regarded as constant in eq. (3.28) (see also 
Part V). The kind of function connecting the va- 
riables Z,; and r in eq. (3.28) is called a homo- 
graphic or linear substitution. It is known 
that a circle in the r-plane is always transformed 
into a circle (or a straight line in special cases) in 
the Z,;-plane (see, for example, [8], p. 605, Sec. 
258) 9. 

To obtain the circle in the Z;-plane correspond- 
ing to the r-circles as indicated in Fig. 9, the follow- 
ing steps are to be performed in conformity with 
eq. (3.28) [8]: 


® Circle diagrams in connection with the determina- 


tion of sound absorption have been introduced also by 
C. W. Kosten [14]. 
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1. Transformation of the origin O, in the r-plane 
to a point P= — (1—p)/(1+p:). 

2. Rotation of the plane through 180° around a 
line through P whose angle with the real axis 
is given by one-half of the angular coordinate 
of m/(1+p)?. : 

3. Inversion with respect to P of the r-circle in 
its displaced position with a constant of inver- 
sion equal to |m/(1+p)?|. 

4. Displacement of the original origin O, to the 
point O= (m/2)/(1+>p), which is the origin 
for the Z;;/zo-circle. 

As these steps include only displacements of ori- 
gins and an inversion, it is obvious that every 
r-circle is transformed into a circle (or straight line) 
in the Z;;-plane. The position of the centre C as well 
‘as the radius @ of the corresponding Z;;-circle can 
be calculated by the exact expression eq. (3.28). 
This calculation will be given in the following sec- 
tion under the simplifying assumption that cot p/2 
is constant and real; this assumption is justified 
under the usual experimental conditions. Fig. 10 
shows schematically the positions of the r-circle and 
the Z;;-circle. Since the real part of Z;; is always to 
be positive, the Z,;-circle must be situated at the 
right side of the imaginary axis of the Z,;-plane. It 
is easy to see that the inner part of the r-circle is 
transformed into the inner part of the Z,;-circle. 


Fig.10. Corresponding circular loci in the r-plane 
(left) for r=Rei® and in the Z4j\-plane 
(right) for Z4;/z3=2x-+jy. The circles are de- 
scribed by r and Z4;/Z) as indicated by the sense 
of the arrows with increasing ”, that is with 
increasing frequency. 


The greatest radius R occurring in the r-plane 
corresponds to the value d=0, as seen from eq. 
(3.27). R equals |I"| in this case; Ryax=1 for a 
perfect reflector. Owing to the absorption in the 
medium, R decreases steadily when d increases. If 
d tends toward infinity, R tends toward zero; d= 
corresponds to a crystal radiating into an infinitely 
extended medium, where no reflected wave exists 
(r=0). The corresponding value of Z;; follows from 


eq. (3,23): 
ea. == - (3.29) 


This value of Z;; represents the electric impedance 
equivalent to the radiation-load of a crystal radiat- 
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ing from one face into an infinitely extended me- 
dium. The sequence of r-circles with decreasing 
radii R is represented in the Z;;-plane by a sequence 
of circles with decreasing radii @. But whereas the 
centre O, in the r-plane remains fixed, the centre C 
in the Z;-plane moves towards the point Cy given 
by eq. (3.29), since the coordinates (xq, yo) of C 
are functions of R. Fig. 11 gives a qualitative pic- 
ture of the correspondence between circles in the 
two planes. Actually, as we saw above, a change in 
d causes r to describe a spiral around O,, whose 
radius decreases slowly if the absorption in the 
medium is not extremely large. The corresponding 
locus of Z;; is also strictly a kind of spiral whose 
“centre” is slowly moving as indicated by the se- 
quence C,— C, —C, —Cy in Fig. 11. 


Fig. 11. Location of corresponding circles in the r-plane 
(left) and the Z1;/zo-plane (right). Cy denotes 
the centre with zero radius. 


In practical interferometry, using a vari- 
able path-length d, only a few consecutive 
cycles in the r-plane are described, whose radii 
change slowly under ordinary conditions. Thus Z;; 
describes the same number of cycles in the Z;;-plane 
with slowly decreasing radii and a small displace- 
ment of C, when a, is not extremely large. A full 
cycle in the r-plane is performed when the change 
in (4xd/2,,) equals 2x (eq. (3.26)), or when 


A(2 didn 1, (3.30) 


—2a,,d 


At constant J’ the initial radius R=|J'|e “™ 


changes into R’=|I'|e “"m@*4 
using (3.30) at constant 2,,, 


R'/R=e “m4 —¢ 4n'm, (3.31) 


in each cycle, or, 


If ay dm <1, R'XR or 24; ~Z,; according to eq. 
(3.23). In other words, when we start at a certain 
distance d corresponding to a certain value of Z;; 
and change d into d+2,,/2, we obtain very appro- 
ximately the same value of Z;;, if %, 2m can be as- 
sumed to be small. Or, a change of d by /,,/2 causes 
Z,;; to describe very approximately a complete cycle 
in the Z;;-plane. This is not the case, however, if the 
product 2,,4,, is relatively large. 

In the case of the fixed-path interferometer, 
where /,, is varied at constant d, the term 


—22.d. ; ° 
e ™ in eq. (3.27) remains unaltered. Here we 
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have to deal with only one circle of constant radius 
in the r-plane, corresponding to one circle with cons- 
tant radius in the Z,;-plane. A change in A, corres- 
ponding to eq. (3.30) causes Z;; to describe ex - 
actly one full cycle, provided the term cot y/2 can 
be regarded as constant in eq. (3.23). This means 
that an increase by one of the number of half-wave- 
lengths filling the path-length d leads to exactly the 
same value of Z;;. 


4. Radius and coordinates of the centre of the 
Z,;-cirele 


In the preceding section we dealt qualitatively 
with the conformal representation of the r-circle 
and the corresponding Z;;-circle. Now we establish 
the analytical expressions for the coordinates x, Yo 
of the centre C and for the radius @ of the Z;;-circle, 
which corresponds to an r-circle with radius R and 
centre O, in the r-plane. In eq. (3.23) we put 
Zy/% =w=2+jy;r=Re; cotgy/2=b (assum- 
ing 6 to be a real constant); m=my+jm,, and 
obtain 


Z+jy— (3.32) 
ae! — (my +jm) ie! me Sam ¥} 

2(1+ Re) — (jm, 2m,) bb Real) 
Setting the real and imaginary parts of the left side 
of eq. (3.32) separately equal to zero and eliminat- 
ing cos®# and sin®? from these two equations, for 
example, by squaring and adding them, we obtain 
the equation for the locus of w=Zj;/z) which cor- 
responds to an r-circle with radius R and the origin 
as centre. The result is 


0. 
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(cot p/2 <1 and real; m, < mg). 
When ‘cot y/2=b=0, these formulas are valid for 
any value of m,/my, as is seen by eq. (3.33). 


Fig. 12. Location of circles in the plane of w=2+j y 
=Z;;/z) corresponding to concentric circles in 
the r-plane with different values of radius R. 


Eq. (3.34) gives immediately 


_ my | m g 
Yor wae ea cot > (3.35) 
2 me ee] 
and = ne c 4 (3.36) 


Thus the centres of the circles for different values 
of R are located on a straight line, whose angle y 


{%—[2m,(1+ R®)]/D?}?+ {y=[2m,(1 +R?) +|mPb(1 —R?) 7 / D2} 


_ 4 rmg? (1+ RY)? + [2 my (1 +R?) + |m PbO —R)P _ [mPO-R) 


D4 


where 
D? = 4(1 — R?) + 4m, 6(14 R2) +m? 6?2(1 — R?), 

In most practical cases the quantities cot p/2=b 
and m,/my can be regarded as small (see also 
Part IV). Ignoring all terms of order b?, m,b, and 
higher, we find the following expressions for the 
coordinates xo, yo of the centre C and for the ra- 
dius @ of the quantity w = Z,;/z9 : 


i , (3.33) 


with the x-axis is given by tany=my,/my. The ra- 
dius @ increases with increasing %¢ towards @= co 
at = co or R=1. The smallest possible value for 
to is given by x =m)/2 with e=0 or R=0. 

Fig. 12 shows the positions of different circles 
for different values of R ; it confirms the more qua- 
litative picture in Fig. 11. 


ACUSTICA 


mt ee F.E. BORGNIS: ACOUSTIC INTERFEROMETER 163 


PART IV 
GENERAL PROPERTIES OF THE ELECTRIC INPUT IMPEDANCE 
OF THE ACOUSTIC INTERFEROMETER 
UNDER VARIOUS EXPERIMENTAL CONDITIONS 


1. Quantitative considerations 


Following the more qualitative considerations of 
Part If] we now give attention to the values of the 
more important parameters in the main equation 
(3.23), as encountered under practical conditions. 


a) The quantity m=my)+jm,: 

The imaginary part m, is very small in compari- 
son with mg under ordinary conditions. Neglecting 
losses in the crystal in eq. (3.17 a), that is assuming 
Oe Co <p Cm, and introducing cm = Wy 4m/2%, we 
obtain from eq. (3.17): 

Bia Gnas (4.1) 
™Mg 2% 
The order of magnitude of a,/fo? (in cm~'s”) is 
known to be 1078 in gases and 1019 in liquids 1°. 
The velocity of sound may be assumed to have 
the order of magnitude of cy, =3- 104 cm/s in gases 
and c,=10°cm/s in liquids. The quantity op Am 
may be assumed therefore to be in the order of mag- 
nitude of 
Om Cm =3°10~%: fy in gases, 

fo 10714: fy in liquids. 
Even up to the megacycle range the quantity a, Am 
can be regarded as small in comparison with unity 
(especially in liquids). Still one must be careful in 
special cases, where the assumptions made above 
may not be justified; for some substances 0, and/or 
Cm may be of larger order of magnitude; then fy 
must be smaller to justify the assumption that 
a As Sk: 

Therefore, according to eq. (4.1), m, can be 
neglected ordinarily in comparison with my and m 
may be regarded as a real quantity | m |= my. 


b) The quantity r=Re 3}; 
The foregoing considerations, which resulted in 
the estimates for a4, according to eq. (4.2), jus- 
tify also the assumption mentioned already in 
Part III concerning R, namely that R, according to 
_ eq. (3.27), is a slowly varying function, when d is 
{ steadily changed by an amount of a few half-wave- 
lengths. Thus the spiral in the Z;;-plane mentioned 
in Part III, Section 3, due to a variation of the path- 
length, contracts very slowly with the change of d 

and may reasonably be replaced by a sequence of 
4 circles, as outlined. 


Onn ‘Amy = 


(4.2) 


10 See [3], for example pp. 464 and 538. 


c) The quantity cotg/2: 

According to eq. (1.11) p=y.J=2n(l/A.) —jacl; 
the crystal is usually driven at or near its fundamen- 
tal mode, that is l~/,/2 and y/2~ (x/2) —j a, A,/2. 
The attenuation in crystals is known to be extremely 
small, so that we may regard ~/2 as a real quantity, 
neglecting the term j a, 4,/2. 

In the neighbourhood of v/2=2/2 the cotangent 
function can be replaced by (x/2) Ad/d,, where 
Aad is defined by 4=1)+ AA =21+4 Al in a piezo- 


electric plate. Usually 4 equals 21 very approxi- 


mately and the cotangent term in the denominator 
of eq. (3.23) is a small quantity which can be ne- 
glected in comparison with the other term. The or- 
der of magnitude of m~mzy is 10~° for gases and 
10-1 for liquids. Thus in gases the weight of the 
cotangent term is considerably smaller when com- 
pared with the weight of the first term, which con- 
tains the factor 1/m, than it is in liquids. The co- 
tangent term may become important when r in eq. 
(3.23) comes very close to —1, that is when 
R=|I'| e ~"m is close to unity. 

It is therefore advisable to choose the driving fre- 
quency so that cot ~/2 equals zero, that is at reso- 
nance for the free bar and antiresonance for the 
free plate (see Part II). If the interferometer is used 
at variable frequency, it will depend upon the 
experimental conditions (\I"| and a, °d), whether 
the cotangent term in the denominator of eq. (3.23) 
has to be taken into consideration or not. 


2. Connection between Z;; and the formulas of 
Hubbard for the acoustic interferometer 


J.C. Hussarp [12] was the first to develop a 
theory for the acoustic interferometer. Upon his 
theory are based all subsequent investigations with 
this device. The difference between his theory and 
that given here is that Hupparp’s method, although 
leading by an ingenious procedure quite quickly to 
his well-known results, nevertheless makes it diffi- 
cult to perceive the significance of the various as- 
sumptions involved. On the other hand, the theory 
offered here starts from the basic piezoelectric 
equations of state and leads in a rigorous and 
straightforward way to more general results. 

Hussarp showed that the reaction of the inter- 
ferometer load upon the electric input impedance of 
the driving crystal can be expressed by an impe- 
dance, which for convenience we call 
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Zy = const(P—jQ) , 


in series with the impedance of the free crystal. 

This is, as it should be, quite in agreement with 

the foregoing considerations, which led to the equi- 

valent circuit of Fig.8. The impedance Z,; in 
our notation must consequently be identical with 

Husparn’s Zy, if the assumptions implicitly involv- 

ed in his calculations are introduced in our eq. 

(3.23). These assumptions turn out to be the fol- 

lowing: 

1. cot y/2=0: The crystal is driven exactly at 

its resonant frequency in the case of a bar 

(cf. eq. (2.4)) or at antiresonant frequency in 

the case of a plate (cf. eq. (2.18) ). The losses 

in the crystal are negligible. 

. m=my; m,=O0: The losses in the medium 
are not unusually large, so that m, can be ne- 
glected in comparison with my (cf. Part IV, 
Section la). 

3. QO =x: The phase-jump of the reflected wave 
in the interferometer path equals 180° ex- 
actly (eqs. (3.8) and (3.27)). 

Introducing these assumptions in the general ex- 

pression (3.23), we obtain 


iw) 


—j0 
Z,, = mom lee ; 
1 5. Liha: (4.3) 


After separating the right side into real and imagi- 
nary parts, we find 
z 1— R? 
Zyi= = 9 ( eS Rie (4.4) 
2 \1+R?4+2Rcost) 
‘ 2 Rsin d 
+j = 5 


L4+R242Rcosd? 
Comparing this expression with Hupsarp’s impe- 
dance 


Zy = const (P—j Q) (4.5) 


and introducing 


P= (2wd/e,)—nxn and R=|T| emt 
according to eq. (3.27) with O=z, we get 


p 1-|rPe te" 
ri Bis ea oO ay ees ; 
2.— 2 2md 
1+|f'e 42m4 Aw aay. at agg 
Cm (4.6) 
Sy eee me Leia 
2|T'le PE sa ig ooo 
(9 Pee Se Wes Cm 
> —4a,,d +) 22,4 2wd 
14+|Pe “™—2|Tle “™*cos — 


Cm 
in conformity with Hupparp’s expressions P,_» and 
Q,-o (in his notation), as given, for example, by 
F. E. Fox ({9], formulas (6) and (7)). The cons- 
tant in eq. (4.5) is found by using eqs. (3.17), 
(3.21) and (3.22), and introducing y=g/l=x/21 
= OV 0c/Ge=HCcV Oc/ Ge /L: 
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m Sy) 
; 4 ie for the bar, 
const = ete Ree 
2 Qm Cm for the plat 
LAH) ae 


We may repeat that Hupparn’s formulas do not 
include the cases where a) the frequency differs 
from c,/21, b) the reflector induces a phase jump 
of the reflected wave different from a, and c) the 
absorption coefficient a, is unusually large. Eq. 
(3.23) is rigorous, covering all possible cases. Still, 
under ordinary experimental conditions, one would 
find no appreciable numerical difference between 
the two formulas. 


3. Maxima and minima of the impedance Z;; 
of the interferometer path and their relations 
to a,,|I"|, and O 


Assuming ordinary experimental conditions we 
shall describe some characteristic features of the 
periodic changes of the impedance Z,;, as given in 
eq. (3.23), due to changes of the path length d or 
the acoustic wavelength /,, . 

By ordinary experimental conditions, the follow- 
ing is meant: 

1. fy = const. 

2. p/2=n/2 ; cot y/2=0. The losses in the erys- 
tal are negligible, so that ~ is real. The cor- 
responding frequency fy is given by fy = ¢./2 1; 
thus v~/2=a2l/A=x/2 and the cotangent term 
in eq. (3.23) vanishes. This means that the 
crystal is driven at its resonant frequency in 
the case of a bar, or at its antiresonant fre- 
quency in the case of a plate (see Part II, 
eqs. (2.4) and (2.18)). 

3. The coefficient of absorption a, in the me- 
dium is not unusually large, so that the quan- 
tity ,/, is small in comparison with unity 
(see eq. (4.2)). Then the locus of Z;;, due to 
a change in d of amount 4,,/2 , may be regard- 
ed as practically a circle. When d remains 
fixed and A,, is changed, the locus of Z,; is an 
exact circle, regardless of the value of a 4m - 

If the assumptions mentioned under 1. 2. and 3. 
are valid, Z;; describes a circle as d or A, is varied. 
The radius @ and the coordinates xc, yg of the cen- 
tre C of the circle described by the dimensionless 
quantity w= Z,;/z) has been given in eq. (3.34). 

For experimental reasons it may be desirable to 
know also the locus of the admittance Y4;= 
1/Z,;. The corresponding radius @ and coordinates 
%c, Yo of the centre C of the circle described by the 
dimensionless quantity w =z Y1; are easily found 
by inverting the w-circle as illustrated in Fig. 13. 
The inversion is effected by w=1/w or in coordi- 
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nates, according to the notation in Fig. 13, by 


#0 ee y 
ee 4, 
pune y eye (4.8) 
Where *+jy=w, and z+jy=w. A circle 


(x—2c)? + (y—Yyo)?=0” in the w-plane corres- 
ponds to the circle (~—2%@)? + (y— yo)? =@" in the 
_w-plane; the values of 2, yq and @ of the inverted 


circle are found by means of eq. (3.34) with 
cot p/2=0: 
a _ 2m 1+ R? See page foe 
~ |m-P 1—R? ign y Ag les bee RR? 
O= ia] — ; (cot~/2=0) (4.9) 
y 
z} 2 5 
0 wee 0 


Fig. 13. Loci of the quantities w—=Zj;/z) (left) and 
W=2) Yi; (right) corresponding to the impe- 
dance and admittance of the interferometer 
path. The sense of the arrows corresponds to 
an increase in frequency. 


Fig. 13 shows that the maximal and minimal va- 
tues of | w| and |w|, occurring at the points A and 
B or Band A, are eee at any m=my)+jm,, by 


0 |max = ye + ye sie 0, 


[20 |aax = Vx? + +4 NG iat O’s 


| w min = = a2 4+ y + i > |w lait — x +9 — QO. 
From eqs. (3.34) and ie 9) we obtain 
zo|m| 1+R 
| Z4i max = =o 1-R’ 
2 1+R 
LaF Bate ape pay zy 23 
[Yiilnex= 2 Tm| IR (4.10) 
1}. _Zolm|1—R 
[Ziyi [ee >. 2 1 aR: > 
Od — 
| V4 |min od im | ime: 
The following relations result: 
Zz |m |? 
|Z; = ; Zi | min = ‘ue | ass 
4, 
| Yi lmax * | Yi nin = (4.11) 
zq° | m| 
[21% loin a7 (Fy | Yuiln e a = (4) 
| Zi ee 1 Bt R ‘ V4 et ax ul + R ; 


The products of the maximal and minimal 
values of |Z,;| or | Y;;| are independent of R, that 
is of |J”| and a,,d. On the other hand their ratio 
is independent of the properties of the crystal and 
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of m, but depends upon the absorption in the me- 
dium. From this ratio the quantity R can be deter- 
mined, which involves the absorption of the medium 
(%m@) as well as the losses (|J'|) at the reflector. 
Measurements at two different values of d allow one 
to separate the losses in the medium from those at 
the reflector, as is well known [10]. Eqs. (4.11) are 
obtained under the assumption that the change in 


R due to the factor e “™” is small; therefore they 
should be applied to adjacent maxima and mi- 
nima at each of the two path-lengths. 

The coefficient of absorption a, and the coeffi- 
cient of reflection | Z’| can be calculated as follows, 
if the ratio | Z4i |min/| Z\i faa or | Vij |min/| Lae hoe is 
known for two different path lengths d, and d, 
Introducing 


a ae Vee oe 
and Kean dan Fils (4.13) 
we have, according to eq. (3.27), 
R=|F le, *m=e-*, (4.14) 
and from eq. (4.12) 
a es tanh 3 or X= its (4.15) 


Denoting by o, and o, the values of o measured at 
the two distances d, and dy, we have the relations 


X,=2 a4 d,—In|F|=In[(1+ 0,)/(1—9,)], 


Xy = 2 am dg—In| I’ | =In[ (1+ 02) /(1—99)], 

whence we find (4.16) 
Se (1 — 94) (1+ 99) 

eae) (ae a) 2 (1 + 04) (1 — 9g) set 

and. | P|=etm@+4)] /(1— 94) (= 99) ; (4,18) 
(1+ 0,) (1+ 6,) 

also |r |= —e%nts} meet (4.19) 

1+ 012 


When d can be made so large that 2 a, 
the following simplified formula may be used, which 
is obtained from eq. (4.16), neglecting In| J’| : 

1 l+o 


ae ae 


An 


(2a,d>In|I"|). (4.20) 


On the other hand if d can be made so small that 


e —2a,,d mls 


\I'| can be measured directly. Eq. 
(4.19) ee becomes 


[F|~ = (2a,nd<1). (4.20a) 


Returning to assumption 3. at the beginning of 
this section, m, is very small in comparison with my 
under ordinary conditions. This means that yc in 
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eq. (3.34) or (4.9) is very small in comparison 
with 2. Thus the centres C and C of the circles in 
the w- and w-planes lie practically on the real axis; 
the maximal and minimal values of Zj; or Y4; are 
practically pure resistances and conductances. Under 
the condition that m, is ignored in comparison with 
my, C and C are located exactly on the real axis 
and the difference in # between the points A and B 
or B and A in Fig. 13 is exactly to 40 =x; from 
eq. (3.27) it is seen that this corresponds to a va- 
riation Ad =x c,/2 @=4,/4. Indeed, according to 
eq. (4.4). Z1; is real when sin? =0 or P= tna; 
in this case | Z4; | reaches its maximum (n=l, 3, 
5,...) or its minimum (n=0, 2, 4,...). Therefore 
the difference in # in proceeding from a maximum 
to the next minimum or from a minimum to the 
next maximum in |Z; | (or | Y,i| of course) is 
equal to a. Or, if my < mg, a change in path-length 
of A,/4 leads from the real maximum of | Zi: | or 
| Yai | to the adjacent real minimum or vice versa. 
This is not so, however, if the absorption in the 
medium is unusually large; in this case the corres- 
ponding change in path-length differs from 2/4 
and the maximal or minimal values of the moduli of 
Z,; or Yy; occur at complex values of Z4; or Yy;. 
Since for m, < my the maxima or minima of | Zi 
or | Yy;| occur at 0= tnx, we find from eq. (3.27) 
that the corresponding values of d=d,, are given by 


2Hdy 


0=itnn= C) 
Cm 
O \Am 
or dn = [n+ = a (m,<mp). (4.21) 
n=1,3,5,... correspond to a maximum and 
n=0,2,4,... to a minimum of | Z,;|. The phase 


jump @ of the reflected wave can be determined by 
means of eq. (4.21), when the distance d,, and the 
number of quarter-wavelengths filling the interfero- 
meter path are known. The difference in dy, 
corresponding to a maximum and an adjacent mini- 
mum in | 24; ' or | ¥4;|, is always 2»/4, whatever O 
may be, if © can be regarded as a constant. 

Eq. (4.21) can be applied also to the fixed- 
path interferometer, where d is kept constant and 
the acoustic wavelength /,, is changed. 

With either type of interferometer, for a perfect 
reflector (O =x), eq. (4.21) states that a maxi- 
mum of |Z;;! occurs always when an even num- 
ber of quarter-wavelengths, that is a whole number 
of half-wavelengths, fills the interferometer path; a 
minimum of | Z;;| occurs at an odd number of 
quarter-wavelengths in the path. If the phase jump 
© differs from a, the number of wavelengths filling 
the path at a maximal or minimal value of | Z,;! is 
no longer a whole multiple of A,,/4. 
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The difference, however, in the number of 
quarter-wavelengths filling the interferometer path, 
and corresponding to a maximum and an adjacent 
minimum of |Zj;| or |¥y|, is always unity, 
whatever the value of the constant O (or |I’| of 
course, which quantity does not enter eq. (4.22) at 
all) may be. 


4. Determination of the velocity of sound c,, 
from the maxima and minima 
of the impedance |Z1i| 11 


As is well known, the velocity of sound c, in the 
medium can be found from the distance between 
maxima and minima of |Z1;| (or | Vy; |). Under the 
assumptions 2. and 3, of the preceding section, 
Z,;/z9 describes a full circle in the w-plane (Fig. 13), 
when the change in # amounts to 2a (see eqs. 
(3.26) and (4.3)). Thus the distance between two 
adjacent maxima or minima of | Z1i | (or | ¥4;|) 
corresponds to a change 


s(i28) (88) a 


= rt (4.22) 
for constant ©, as seen from eq. (3.27). 
If the variation of Z;; (or Y,;) is caused by a 
change of path-length Ad , we obtain from eq. (4.22) 
Cn= — Ad =2 f Ad 


as (m,<my), (4.23) 
where Ad is the change in path-length necessary to 
proceed from one maximum of | Z;;| to the adjacent 
maximum or from one minimum of | Zi | to the ad- 
jacent minimum. 

If the variation of Z,; (or Y4;) is produced at 
constant frequency fy by a small change 4/,, in the 
acoustic wavelength /,,, owing, for example, to a 
change of pressure, temperature or composition of 
the medium, the corresponding change Ac, is found 
from eq. (4.22) to be 

Sis Aen __ Ulm dom 
Cu Onde 2am 
This is the relative change in c,, corresponding to a 
change in 2, by an amount Alp =Acm/fy great 
enough to make | Z;;| change from one maximum 
to the adjacent maximum, or from one minimum 
to the adjacent minimum. 

If the variation of Z;; (or Y,;) is produced by a 
change of the driving frequency f= w/22, the ne- 
cessary change 4/ to proceed from one maximum of 
'Z;:| (or | ¥,;!) to the next or from one minimum 
to the next is, according to eq. (4.22), 

Af =Aw/2x=c,/2d, 


whence the velocity cy is 


Cn =2d Af. 


(4.24) 


(4.25) 


11 See also [4]. 


ee . = i 7 ‘ . al 
~ * : : ; 
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A change of 2 x in # always corresponds to a full 
“period” of Z;;, that is to a closed circle in the 
w-plane; it is obvious that the foregoing formulas 
(4.22) to (4.25) hold not only for adjacent maxima 
or minima of | Z1; | but also for the more general 
case in which one proceeds from an arbitrary 
value of Z;; (or Y4;) to exactly the same adjacent 
corresponding value. 

Still, as pointed out in the preceding section, 
under ordinary conditions (m,<my) 
the change in # necessary to proceed from a ma xi- 
mum of |Z;;| to the adjacent minimum equals 
a, that is just half the amount found necessary to 
proceed from a maximum to the adjacent maxi- 
mum in the general case. The foregoing formulas 
‘can therefore be modified for the case in which a 

maximum of | Z1i | and an adjacent mini- 
mum are chosen for the determination of c,,. The 
formulas read in this case as follows: 


Cm =4 f+ Ad’, (4.23 a) 
ees (m,<m))  (4:24a) 
Ga AN i fe \ ; 

Cn =4d: Af . (4.25 a) 


5. On the form of the periodic curves 
for |Z,;| or | ¥4;| 


Assuming again m, < mp, that is circles in the 
w- and w-planes (Fig. 14) with centres on the real 
axis, we now see what happens if @ changes, for 
example, from %,=0 at B (or B) to Jp= +2/2. 
This change in @ corresponds to a change in d of 
+4/8, as is seen from eq. (3.27). The correspond- 
ing values of w=Z,;/z) and w=2z) Yy; are found 
from eq. (4.3): 


: Ee 
teed p= IE (406) 
2 17jR my, 1+jR 
whence 


“eee: Greces: 25] | wp |? = ee 
0S ee rr 
By means of eqs. (3.34) and (4.9) the coordinates 

of wp =%p+jyp and Wp = fp +jyp can be written 


as 
Res 2 yk 
erty ht’ 14+R?°? (4,20) 
1-R 1—R 
pe See Se : 
yo= = 7 Ree YD 1a Re 


Fig. 14 shows the location of wp and wp, denot- 
ed by D and D, for different values of R (cf. the 
notation in Fig. 13). The circles vary from large 
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radii for R=1 to small radii for R<1; at R=O 
the radii @ and @ are zero (eqs. (3.34) and (4.9)) 
and the “centres” are at x= mp/2 and xq =2/m). 
According to eqs. (4.27) and (4.28) the moduli 
| wp |=mp/2, and |wp |=2/my, are independent of 
R. The circles with O and O as centres and radii 
m,/2 and 2/my therefore intersect the w- and w-circ- 
les for different values of R at the points D and D 
whose coordinates are given in eq. (4.29). 


gy 


Fig. 14. Location of the points D (impedance) and D 
(admittance) for ®9=2na+zx/2, and of D’ 
and D’ for #=2nn—n/2 in the w- or 
w-plane. The points A or A correspond to 
0 = (2n+1) x, the points B or B to P=2 nz. 


Fig. 14 can be used to represent the w-plane as 
well as the w-plane with the notation given in 
Fig. 13; only the scale-values have to be chosen pro- 
perly. Incidentally, if mp =2 the scale values are the 
same for both planes, as seen from eqs. (4.27) and 
(4.28). 

From Fig. 14 qualitative curves for (2/my z9)| Z,;| 
or (mp 29/2)| ¥4;| plotted against ® can be found, 
as shown in Fig. 15. If R is close to unity, the w- or 
w-circle is large, and D or D (or, correspondingly, 
D’) are much nearer to the minimum B or A than 
to the maximum A or B of (2/mpz»)| Z| or 
(mo Zp/2)| Yi;|. At any R the quantity (2/m zy) | Z1;| 
(and also (mp z,/2)| Yi; |) has the same value at 
= + 2/2. Thus at R~1 the curves drop down 
rather rapidly from their maximal values to the 
value 1 at #=x/2 and have comparatively flat mi- 
nima. 

When R <1 the circles in the w- or w-plane are 
small and the points D or D have nearly the same 
distance from B as from A (or from A as from B). 
Thus, as R decreases, maxima and minima assume 
more and more closely the same shape. At R=0, 
corresponding to a reflector located at infinity or to 
a perfect absorber, the periodicity vanishes altoge- 
ther, and a constant value of | Z;;| or | ¥y;! is reach- 
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ed, namely the value for a crystal radiating into an 
infinite medium (eq. (3.29)). 

It may be noted that the curves for (mp 2/2)| Y1i| 
are obtained by shifting the curves for (2/mg Z9) | Z1i| 
to the right or to the left by an amount of ® =< in 
Fig. 15, since 


[ (2/rmy 29) | Z1i pce = [(megjzo/ 2) | Yas lle 
according to eq. (4.3). 


2 

Fig. 15. Representation of (a) sera |Z4:| and (b) 
0 0 

mou Yii| versus ®+ 2nzx for different va- 


lues of R. 


The curves in Fig. 15 permit the following con- 
clusions. At the same path-length d and the same 
reflection coefficient |J"| the curves for |Z;| or 
| Yi; | show sharp maxima and flat minima as long 
as the absorption ,, is small. With large absorp- 
tion this difference between maxima and minima 
tends to disappear. Or, for the same medium and 
the same reflector, the curves for | Z| or |Y,;| in 
the neighbourhood of two different distances d, and 
dy such that a, dy <1 and a, d; > 1, show a sharp 
maximum and a flat minimum at d,, whereas at d, 
this difference disappears. Different curves for | Z;; | 
or | Y;;| always intersect at whole multiples of 
&=+x/2; the multiplicity can be found by the 
number of quarter-wavelengths along the path- 


length d. 
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We may add the following remark at this point: 
the total input-impedance Z, of the interferometer 
is given by the sum Z;,-++Z4;, as shown in Section 2 
of Part III. In all cases where the crystal is driven 
at constant frequency, the constant term Z;, has to 
be added to Z;;. When y/2 =x/2, as we have been 
assuming, Z;, is the real “loss-resistance’” due to 
losses in the crystal and mounting. Z,; contains the 
factor mp2 (eq. (4.3)), and is therefore propor- 
tional to m//?, as seen from eqs. (3.21) and (3.22). 
Since my is about 4000 times smaller in gases (at 
atmospheric pressure) than in liquids, at the 
same frequency the crystal is considerably 
more loaded by a liquid than by a gas. 

Under ordinary circumstances, that is at frequen- 
cies not exceeding the order of magnitude of 
1 Mc/s, a gaseous medium contributes essentially 
in the neighbourhood of | Z;;| to the total impedance 
Z,, whereas in the neighbourhood of | Zi; ce the 
value of Z,, is predominant. Fig. 16 illustrates this 
case. The distance OB, corresponding to | Z4; Nn as 
only a small fraction of 0,0, corresponding to Z;,. 


0; 


Fig.16. Total input impedance Z;=Zj.+Z1; for a 
gaseous medium. 


The minimum of Z, is chiefly determined by Z;,, 
and is therefore rather flat. The sharp maximum of 
Z, is due to the sharp maximum of | Z;;|, as shown 
in Fig. 15a, occurring at O=n, 3n, 50,... or 
d,,=n,,/2 according to eq. (4.21). 

In liquids, owing to the much greater value of 
mg, the crystal is considerably loaded, up to very 
high frequencies, even when the impedance | Zi | 
has its minimal value, because the ratio | Z,;|/| Z;¢ | 
is so much larger than in gases. The position of O, 
in Fig. 16 is then rather near to O. The behavior of 
|Z,| is therefore approximately the same as that of 
| 24: | alone (Fig. 15) ¥. 


2 See also [9], (Fig.2c). The characteristic difference 
in the shapes of the curves for current versus d in 
gases and in liquids, as experimentally observed and 
described in that paper (see also [3], Fig. 186, p. 245) 
is a consequence of the characteristic features of | Zi; |, 
as described above and also mentioned by Fox. The 
particular curves for the current as given by Fox are 
essentially due to the particular circuit connected with 
his interferometer and to the frequencies applied. 
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6. Locus of minima and maxima of | Z;;| and | Y,;| 


If the path-length d is steadily increased, the mi- 
nima of | Z;;| or | Y4;| rise and the maxima fall con- 
tinuously, both finally tending towards the same 
value as d approaches infinity. Assuming the same 
conditions as mentioned in Section 3, and making 
use of eq. (4.10), we see that the locations of the 
minima and maxima of both |Z; | and | Y4i| are 
associated with the function 


{(d) = 1—-R _ 1—exp(—2am,d+In|J"|) 


14+R  l+exp(—2a,d+In|I')) 


or j(d) =tanb| 39 d— 3 In|). (4.30) 


(a) 1/f(d) 
(b) f(d) 


— 


_ Ee higir par etens 


Fig.17. Locus containing (a) the maxima and (b) the 


Indeed, f(d) contains all the minima of 
| Z4;|+2/(zo|m|) or of | Yi; | * Z| m|/2, whereas the 


maxima of these quantities are on the reciprocal . 


curve 1/f(d) , as seen from eq. (4.10). 

Fig. 17 shows the curves for f(d) and 1/f(d) 
according to eq. (4.30). The origin O is at a di- 
stance of (1/2) In|I’| from the origin O’ of the 
curve for f’ =tanha,d. If |r| is close to unity 
and a,, d small, we can set | J"|=1—«, where ¢ is a 
small quantity, and obtain from eq. (4.30) as a first 
approximation 


f(d) <a, d + 5 (4.31) 
(dential mele 15) 


It will be noted that the minima are located approxi- 
mately on a straight line, which intersects the ordi- 
nate-axis at {(d) =¢/2 = (1/2) (1—|I"|). The maxi- 
ma lie on a hyperbola 1/f(d). This corresponds to 
remarks, found in the literature, that the “peak en- 
velope” is rather a hyperbola than an exponential. 
The envelope is never a simple exponential, as seen 
from eq. (4.30). 


It may be noted that our statements concerning 
maxima and minima are related to |Z; | or Yui 
and not to the total impedance | Z,|=|Z:.+ Zi; 
which is necessarily the observed impedance. Whe- 
ther or not |Z;,| may be neglected depends on the 
actual ratio | Z;,| to | Z;;|. In many cases, especially 


a m | Zo 
minima of |m| 2. Z| and erat Yii| ver- at small values of m and when R comes close to 
sus Omd. unity, proper consideration must be given to | Zte he 

PART V 


THE VARIABLE-FREQUENCY INTERFEROMETER 


For those cases in which movable parts should be 
avoided, the fixed path interferometer offers con- 
siderable experimental advantages [13]. The change 
of |Z,;| at constant path-length can be effected by 
a variation of the acoustic wavelength /,, in the 
medium, as seen from eq. (3.26). This variation 
might take place at constant frequency fy, as, for 
example, when the pressure, the temperature, or the 
composition of the medium is varied. From this 
variation the corresponding change Ac,, of the velo- 
city ¢» can be found, as shown by eq. (4.24). Still, 
the absolute initial value of c, must be known in 
order to plot the changes Ac,, against the variable 
parameter that causes /,, to change. 

A determination of c,, itself can be accomplished 
by varying the driving frequency at constant path- 
length. If one could observe the corresponding 
change in | Z;;|, cm could easily be found by means 
of eqs. (4.25) or (4.25a). Also the absorption- 
coefficient a,, could be calculated by formula (4.17) 


or (4.20). The difficulty here is that the quantity 


determined experimentally is Z,)=|Zre+Zil, 
where now Z;, also varies due to the change in fre- 
quency. 


At constant d the locus of | Z,;| is practically a 
circle, whose centre, under ordinary conditions 
(m,<my), lies practically on the real axis. The 
reason why this locus is not strictly a circle is 
seen from eq. (3.23): by changing the frequency 
the term cot g/2, which is zero for the reference 
frequency fy=¢,,/21, deviates from zero. Accord- 
ing to eq. (3.34) the centre C of the circle for Z;; 
moves away from the x-axis, as the frequency, and 
consequently @ , changes. Still, under most practical 
conditions the “coupling” term cot /2 in the deno- 
minator of eq. (3.23) may be kept so close to zero 
that it can be neglected in comparison with the term 
2(1+r)/m(1—r) (see also Part IV, Section 1c). 
Quantitatively we may say that the cotangent term 
may be ignored, if 


170 
x | Af| < ae 1—R 
2 fo Soe 

ng A 2 EN eae) (Gh) 
fo | m | 


A more serious complication, however, is involv- 
ed by the fact that the crystal gets off resonance 
when f is varied. In the neighbourhood of the refe- 
rence-frequency fy we can make use of the equiva- 
lent circuit of Fig. 8 and write 


1 : 
Zye=Ro+j(wL- | = Rd 42h Aw G2) 


C 
This situation is illustrated by Fig. 18. As f changes, 
the origin O, for the impedance Z; = Z4. -+ Z;; moves 
perpendicularly to the x-axis. Let 0,’ be the location 
of O, at dw = Aw, . It is obvious that the maximum 
(or the minimum) of Z, now differs by 0,’ E from 
the maximum O,A (or the minimum O,B) belong- 
ing to dw =0. (E is the intersection of 0,’C with 
the circle of radius CO, around the centre C). Fur- 
thermore the difference p’ — #4’ between the values 
of # belonging to B’ and A’ deviates from x if A’ 
and B’ move away from A and B. Thus both the de- 
termination of a, (associated with |Z; |max and 
|Z1i|min) and of cy, (associated with 40 = Jy — Oy’) 
are affected. 

Adequate analytical expressions might be calcu- 
lated from eqs. (3.18) or (3.19), giving the neces- 
sary corrections for the determination of a, and 
C,. This requires, however, more detailed calcula- 
tions, which will not be dealt with here. Corrections 
may be ignored if the angle ¢ in Fig. 18 is suffi- 
ciently small. This condition can be expressed by 


2L Aw 


R, + 2 Xo 
aaL? (5,8) 


where, according to eqs. (3.21), (3.22), (3.27) 
and (3.34) 
Vo? Qelt 14 R? 


Mo Ta oa 1_R for the bar, 


Ee oveie, Flee bbe 
Mp 4A 1_R2 for the plate, 
and R= | i: i 


For example, the necessary change Af to describe a 
full circle of Z;;, corresponding to 40=22, is 
given from eq. (3.27) by 
Af _ An/2 
fo d 
To keep 4/ as small as possible, the path-length d 
should contain a large number of half-wavelengths. 


(Ad =2 2). (5.5) 
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Fig. 18. Locus of Z1=Zi¢+Z4; at variable frequency. 
As f is changed, the origin O moves along the 
dotted line normal to the z-axis. 


On the other hand d must not be too large, because 
flat maxima and minima are undesirable. Thus a 
compromise must be reached with respect to a suit- 
able choice of d. 

‘A qualitative consideration concerning the error 
in determining c, from the spacing between two 
maxima or minima of | Z,| can be given as follows. 
A change in frequency changes ? according to eq. 
(3.27). From Fig. 10 it is seen that an increase in 
& corresponds to a clockwise movement of r around 
the circle in the r-plane and consequently to a clock- 
wise movement of Z;; in the Z,;-plane. Moreover, an 
increase in frequency changes the position of the 
centre O,’ in Fig. 18 in the direction of the nega- 
tive y-axis. Let O,’ in Fig. 19 correspond to the 
maximum A’ of | Z,| ; by a positive change 4m the 
centre O,’ may be shifted to 0,” corresponding to 
the adjacent maximum A” of | Z,|. As may be seen 


Increasing frequency 


Decreasing frequency 


Fig. 19, Location of adjacent maxima A’ and A” or A* 
- and A** of Z, corresponding to a departure 
Aq@ from the reference frequency (4@=0). 


from Fig. 19, the corresponding change in #, lead- 
ing from A’ in a clockwise direction to A”, is 
less than 2 (since } = 2 would correspond to a 
full turn around the circle, leading back to AY. 
The change Af in frequency, which, according 
to eq. (3.27), is proportional to the change 
in #, is therefore less than it would be if the 
centre O,’ did not move. In other words, the 
maxima A’ and A” are spaced nearer together than 
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they would be if the crystal impedance Z,,=0,/O 
did not change. If 4f is inserted into eq. (4.25), a 
value of ec, smaller than the correct value results, 
because formula (4.25) is based on the assumption 
that the change in 7, corresponding to adjacent 
maxima (or minima) equals 2x exactly. 

As the frequency decreases from the refe- 
rence-frequency (Aw =0), the centre O,* in Fig. 19 
may shift to O,**, the corresponding maxima of 
|Z, | being A* and A**, A decrease in f corresponds 
to a decrease in , and the circle for Z,; is describ- 
ed counter-clockwise from A”* towards A**, 
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Here too the change in # between A* and A** is 
less than 22, and the maxima or minima are spac- 
ed more closely than they would be if # equaled 
2%. 

Thus the value for c,, found from formula (4.25) 
is less than the correct value, when derived from the 
maxima or minima on either side of the reference- 
frequency. The variable-frequency interferometer 
involves in principle an error of this nature 
in the determination of c,; the error might be 
made negligible if the shift of the centre O, could 
be kept sufficiently small. 


APPENDICES 


A. Numerical values for the constants of an 
X-cut quartz crystal 


1. The bar in lengthwise vibration 


The numerical values needed for the bar are sh. 
dj, and ¢,”7, as seen from eqs. (1.1) and (1.2). 
The constants, converted to mks units, are taken 
from W. G. Capy [5] '*. The page references below 
are from this book. 
sk, = s¥ =1.27- 1071! m?/newton (pp. 55 and 

137) , 
dy» = — dy, = — 2.3: 10-1” coulomb/newton 
(p. 219) , 
eT —4.0-107"! farad/m (k; on p. 414) . 


From these quantities we derive 
gq’ =1/s¥,=7.9 - 10! newton/m? , 
H =d,/s¥, = — 0.181 coulomb/m? , 
H2/eT g® =1.04- 107. 


Thus the term H?/e” g® can be regarded as small in 
comparison with unity, as repeatedly assumed. The 
relative difference between the frequencies for anti- 
resonance and series resonance for the free bar in 
the fundamental mode (n=0) is given by 
eq. (2.10): 

Wa — Ws 4o 4 HF 


=" = ~4°1073. 
Ws Ws Te oO” 


2. The plate in thickness vibration 


The constants needed for the plate are oo ae 
H =e,,, and ¢,5, as seen from eqs. (1.25), (1.26), 
(1.27), (1.28): 

gq’ =c¥, =8.7- 10" newton/m? (p. 137) , 

H =e,, =0.173 coulomb/m? (p. 219) , 

e* =3.9-10~" farad/m (ky on p. 414). 


18 Numerical values are also given in [3], pp. 90 
and 91 and [15], pp. 84 and 85. 


From this we obtain 


H?/e g® =0.88- 10-2, 


HP? 
qe =i l at #3| = 8.8- 101° newton/m? . 
q 


Thus H?/e" g? , which occurs in the treatment of the 
plate, is small compared with unity. The values for 
g® and H are somewhat different for the plate and 
the bar, as stated already in Part I, Section 3, but 
still of the same order of magnitude. 

The relative difference in frequency between se- 
ries resonance and antiresonance for the free plate 
in the fundamental mode (n=0) follows from eq. 
(2.19): 

W,—-W, Aw Ay dak 


=- So Omen. 
Wa Wa a 6S GE 


The value of coty/2 at series resonance 
of the free plate, neglecting crystal losses, follows 
from eq (2.15): 

ol 2 HT? 
elo oat aa, Tat ge 
(2n+1) x Aw 
2 Oa 


and for the fundamental mode (n=0) 


A 
cot p/2= ~ —" -5.5-10-8. 
Ais (0 
The slight departure of cotq/2 from zero corres- 
ponds to the fact that the thickness of the free plate 
at series resonance is slightly different from /,/2 in 


the fundamental mode. Indeed, we have 


Mgt _ (2n+1) a ( 1 aa) 


Digs 2 


and hence for n=0 


Wa, 


a 4 (1 =| — 0.9965 


Os 


As 
= 


€ 
s 


Actually, as is well known, the frequency measured 
experimentally for series resonance of a free plate 
is somewhat higher than the theoretical value, owing 
to the fact that the plate is of finite area and not infi- 


— 


Sn 
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nitely extended as assumed by the theory. The slight 
departure of 1 from /,/2, as found above, is was 
fore chiefly of theoretical interest. 


B. Numerical examples of interferometer 
impedances 


As an illustration of the foregoing theoretical con- 
siderations we offer two numerical examples, one 
for water and one for air. In each case the ultra- 
sonic source is an X-cut quartz plate. 


1.Medium: Water, fp =15 Me/s 


The following data are used: 
0, = 2.65 + 10? kg/m? , 
Co =5.8- 108 m/s, 
Om = 1.0° 10% kg/m, 
m =1.5:10% m/s, 
7p = On Cpa] Oo Coe OL AEG (Sel ave 
A=lcom?=10'm?, |I'|=0.9, 
2o,=0.11cem7!=1llm 
From appendix A 2 we have H =0.173 coulomb/m? 
and q? =8.8- 101° newton/m?. 
The thickness of the plate is 1=/,/2=c,/2 fy= 
1.9: 10~4m, assuming the plate to be driven at its 


fundamental antiresonant frequency (eq. (2.18)). 
zy is found from eq. (3.22): 


The wavelength in water is A4n=Cm/fy=10~4m 
=0.0lcm. The quantity a, 7, =0.00055 <1; ac 
cording to eq. (4.1), m, can be neglected in com- 


parison with my and the centre of the circle for Z;; can 
be assumed to be on the real axis. A change in path 
length of 4,,/2 , corresponding to a full turn around 
the Z;-circle causes a contraction of the radius R in 


the r-plane by the factor e ™™—0.99945; evi- 
dently this contraction can be ignored (see Part IV, 
Section 1b). 

We choose several different values for the path- 
length d; the corresponding R is found from eq. 
(3.27). The location of the centre xq and the radius 
o of the (Z;;/zo)-circle follows from eq. (3.34) ; 
yo=0, since m,=O-and cot p/2=0. Maxima and 
minima of |Z,;| are then given by z)(aco+@) or 
by eq. (4.10). Table I shows the results. 

It is seen that the crystal is loaded considerably 
by the water-column even at the minimal values of 
|Z,;|, as has been remarked at the end of Part IV, 
Section 5. The resistance Z;, of the free crystal may 
be assumed to have a value between 100 and 1000 
ohms; it is of the same order of magnitude as 
|Z1i\min and can therefore not be ignored for fre- 
quencies of this order of magnitude. Since the order 
of magnitude of |Z4;| is proportional to mo/f,? 
(eq. (3.22)), |Z4;| extends to much larger values 
when the frequency is lowered, so that the influence 
of Z;, may become negligible at lower frequencies. 


2. Medium: Air at atmospheric pressure and 


0° C, f=1Me/s 


Using 0» = 1.3 kg/m? and cy, = 331 m/s one finds 
my =0.28:10~*. The reflection coefficient | J” | may 
be computed according to a formula given by Herz- 
FELD [10]: 


Table I 
Water, fp =15 Mc/s, area=1 cm”, X-cut quartz plate (J=A,/2). 

| - | | | ae hinase 

d 2 Ua | R XC Zo QZ [74 | paase |) 1748 Fata pera 

. | 71 | min 

x cm ies tes ey ohms ohms ohms ohms 
1 200 | 0.81 22150 21670 43810 480 91.2 
2 400 0.72 14500 | 13750 | 28250 750 37.7 
5 | 1000 0.52 8010 | 6560 14570 1450 10.0 
10 | 2000 | 0.30 5510 3030 8540 2480 3.4 
oe) | © te) 4600 0 4600 | 4600 1.0 
Table II 


Air at atmospheric pressure, fy =1 Mc/s, area=1 cm?, 
X-cut quartz plate (J=A,/2). 
] [|Z 


li | max 
d | 2 d/Am | R | XC 2 | Q2 | Ani Lae | Zyj Pe / [4s is 
| ‘min 
em / ohms | ohms ohms | ohms 
1 41 2 80.8u), » 0:64 717 651 | 1368 66 | 20.7 
2 ) 21g 0.43 437 317 754 + 120 6.3 
5 | 303.0 | 0.13 311 79 390 232 1.7 
10 |} 606.1 | 0.018 | 301 11 312, | 290 1.07 
a oa 0 | 300 | 0 300 | 300 1.00 


ae 
’ 
€ 
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| £|=1-—4.4- 10-5 yf =0.96. 


The coefficient of absorption may be found from 
[3], Table 64, p. 540: j 
Om = 2° f?- 10-8 =0.2 cm! or 20,=0.4cem'. 

The wavelength in air is Ay =cy/f) =3.31:1074m 
=3.31°10~? cm. The quantity a, 4, = 0.0066 <1; 
thus m, as well as the contraction of a circle in the 
r-plane due to a change of /,,/2 in d, can be ignor- 
ed (Part IV, Section 1). Table II gives the numeri- 
cal values, computed in the same way as for water 


in Table I. 


Owing to the considerably smaller value of mg in 
comparison with the foregoing example for water 
the crystal is much less loaded by the air than by 
the water, although the frequency chosen in the case 
of air is 15 times smaller. This latter fact increases 
the ratio | Z;;|/|Z:.| by a factor (15)? = 225. When 
the frequency is further lowered, |Z; | rises to 
higher values by the factor 1//? and the crystal be- 
comes heavily loaded even by a gaseous medium. 


20000 
0 500 


40000 Q 
1000 2 


A=l1cm’, 
500 ohms, (a) water at 15 Me/s, (b) air at 
1 Me/s. 


Fig.20. Impedance diagrams for R.= 


Fig. 20 illustrates the two cases for water and air. 
In (b) the vector construction for Z,=R,+Z,; is 
indicated. A similar construction would be used in 


(a). 
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C. Note on the “free” (non-radiating) face 
of the piezoelectric crystal 


All the results of the interferometer theory, as 
presented here, have been based upon the assump- 
tion that the crystal radiates into the medium within 
the interferometer path only, whereas the opposite 
surface is ,,free“, that is, non-radiating. 

Actually the latter surface is usually bounded by 
air, rather than by vacuum, and the foregoing nume- 
rical example for air indicates that, especially at 
lower frequencies, the crystal may be considerably 
loaded at a surface bounded by a gaseous medium. 

If the calculations in Part! are carried through 
with a modified boundary condition at x=0, cor- 
responding to a load at x=0 represented by an im- 
pedance Z,,' (or 2m’), the following expression for 
the impedance Z,; replaces eq. (3.20) **: 


/ 
YP 2m Sm __ 


] —2 j cot on ates 7 
pes Zo am : 2m (C1) 
—; —jcot — 
2m + 2m : 2 


At resonance (/=/,/2) the cotangent terms vanish 
and it is seen that z,' simply constitutes an addi- 
tion to z,, the acoustic load of the interferometer 
path (cf. eq. (3.20)). Consequently at resonance 
we have to add an impedance Z,;’, due to the gase- 
ous medium at the free surface, to the impedance 
Z;; of the interferometer path considered before. If 
a plane reflected wave on the free side of the crystal 
can be avoided by the use of highly absorbing mate- 
rial, or suitable configuration of the reflecting sur- 
face, or by having the path sufficiently long, Zi 
equals my 2/2 (eq. (3.29)). As seen from the pre- 
vious example for air, Z,; then equals 300 ohms at 
f=1Me/s for an area A of 1cm?. Z;;' increases 
with 1//?, becoming 30000 ohms at 100kce/s. If 
the radiation into the free side strikes some obstacle 
and produces a reflected plane wave, the gaseous 
path behaves like another fixed-path interferometer 
in series with the true interferometer, and Z,;/ may 
induce a considerable complex impedance in addi- 
tion to Z;;. Proper attention therefore should be 
given to the “free” surface of the crystal. The effect 
of Z;; is minimized when a plane reflector parallel 


14 Tt may be noted that eq. (C1) also describes the 
effect of a piezoelectric plate which is inserted nor- 
mally to the direction of propagation of a plane com- 
pressional wave and which is loaded at its input by an 
electric impedance Zap. One only has to solve eq. 
(C1) for zm and to replace Z;; by the expression 
[Zap/(1—j @ C Zap)]—Zic (see Fig. 8). However, the 
formula so obtained applies only to the case in which 
the centre of gravity of the plate is considered to be at 
rest. 


ers 
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and close to the free crystal surface is located at a 
distance of (2n+1) A/4. 


The author is indebted to Dr. W. G. Capy for 


many stimulating discussions. 
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Symbols 


electrode area, 

cot ~/2, 

parallel and series capacitance in the equi- 
valent circuits of Figs. 3 and 5, eqs. (1.23), 
(1.43), (1.45), 

C in Fig. 6, specialized for bar and plate, 
eqs. (2.6), (2.17), 

velocity of sound in crystal and medium, 
electric displacement, 

length of interferometer path (Figs. 1] and 7), 
interferometer path-length corresponding to 
a maximum or minimum of the electric 
impedance Z4;, eq. (4.21), 

electric field-strength, 

dissipation constant, eq. (1.9), 

frequency = /2 x, 

piezoelectric constant, eqs. (1.3), (1.4) and 
Ge eCle28):, 


current at crystal input, 


= L in Fig. 6, specialized for bar and plate, 


eqs. (2.6), (2.17), 

length of crystal bar (Fig. 2) and thickness 
of crystal plate (Fig. 4), 

Ym m/c c= My + jm, eq. (3.17), 

Ch ASRS. 

elastic stiffness, eqs. (1.3), (1.4), (1.27), 
(30), 

geo +j w F =¢® for a bar or qg? for a plate, 
eq. (1.9), 

Qm0+j oF, eq. (3.3), 

R in Fig. 6, specialized for bar and plate, 
eqs. (2.5), (2.16), 

loss-resistance of the crystal due to absorp- 
tion and mounting, 

|r| e ant, eq (3.27), 

Re J®, eq. (3.26), 

strain and stress, 

thickness of the crystal bar, Fig. 2, 
sinusoidal driving voltage at crystal input, 
x+jy=Zii/z,, eq. (3.32), 

E+jY=% Vii, eq. (4.8), 

coordinates of the centre of a circle in the 
w-plane eq. (3.34), 

coordinates of the centre of a circle in the 
w-plane, eq. (4.9), 


= 1/Z electric admittance, 


electric impedance, Figs. 3, 5, 8, 
mechanical impedance of the interferometer 
path, eqs. (1.14) and (3.1), 

electric impedance of the free crystal, Fig. 3. 
electric impedance of the interferometer 
path, Fig. 8, eq. (3.23), 

P—j Q Hupparn’s electric impedance of the 
interferometer path, eq. (4.5), 

see eq. (3.21) and (3.22), 

(w/e de) Zm, eqs. (1.16), (1.36), 
amplitude absorption coefficient in space 
for crystal and medium (in PartII a is 
written in place of ae), 
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Yo = (@/cc) —ja=WVoc/q® for the bar, and 
@ Vo-/qg? for the plate, eq. (1.11), 
Ym = (w/cm) —j am, eq. (3.24), 
if = |I'|ei® reflection coefficient, eq. (3.8), 
Ete permittivity, eqs. (1.2), (1.26), 
o) phase angle of the reflection coefficient, 
eq. (3.8), 
o = 2w(d/em) —O=4x(d/dm) — O, eq. (3.27), _ 
hesAan acoustic wave length in crystal and medium, 
Oe, Om density of crystal and medium, 
On Gas radius of circles in the w- and wz-planes, 
eqs. (3.34), (4.9), 
He /| Zi Amim_ ] /| ¥1i loin anit 
2 | [Za ae | Vii ae ? ( )s 
p = yel=w(l/ec) —j acel=2 n(l[Ac) —jacl, 
Ws angular frequency for series-resonance of 
the free crystal, eqs. (2.4), (2.15), 
Wa angular frequency for antiresonance of the 
free crystal, eqs. (2.9), (2.18), 
a) = 2quf. 
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ZUR SYSTEMATIK DER ELASTISCHEN EIGENSCHWINGUNGEN 
ISOTROPER KREISZYLINDER 


von E. Srrrie 


1. Physikalisches Institut der Technischen Hochschule, Stuttgart 


Zusammenfassung 

Die elastische Grundgleichung des isotropen Festkérpers liefert fiir den Grenzfall des 
unendlich langen Kreiszylinders aufer den einfachen Lésungen, welche die Typen der 
Torsions-, Longitudinal- und Biegeschwingungen beschreiben, auch noch kompliziertere. 
Es wird gezeigt, daB insgesamt eine dreifach abzahlbare Schar von Resonanzen auftritt, 
sich also eine jede Resonanz durch ein Tripel ganzer Zahlen (n, p, g) indizieren ]aBt. Diese 
Indizes lassen sich mit den Schwingungsknoten auf der Zylinderoberflache in Beziehung 
bringen und sind durch diese experimentell zu bestimmen. 

Da sich das Randwertproblem fiir den Zylinder endlicher Lange nicht streng lésen 1aBt, 
sind bei der experimentellen Nachpriifung gewisse Abweichungen zu erwarten. Fiir den 

Spezialfall der Resonanzen mit verschwindender Axialkomponente kann man jedoch durch 

den Vergleich mit den Werten der entsprechenden Resonanzen einer diinnen Kreisplatte 
den Fehler eingrenzen. 


Summary 

The fundamental equations for the isotropic elastic solid in the form of an infinite 
circular cylinder yield, beside the simple modes of vibration (longitudinal, torsional, 
flexural), more complicated modes. Altogether there is a threefold band of resonances 
with each resonance characterised by a triplet of integral indices (n, p, g) which are related 
to the nodes on the cylindrical surface and can be experimentally determined. 

As the boundary problem of the cylinder of finite length cannot be rigorously solved, 
some departure from experiment is to be anticipated. For the special case of resonances 
with vanishing axial components, the discrepancies can be limited by comparison with the 
corresponding resonances of a thin circular disc. 


Sommaire 

L’équation fondamentale de l’élasticité du corps solide isotrope admet, dans le cas 
limite d’un cylindre circulaire de longueur indéfinie, non seulement les solutions simples 
gui correspondent aux types d’oscillations de torsion, longitudinales et de flexion, mais 
encore des solutions plus compliquées. On montre, dans le présent travail, qu’il apparait 
en tout une famille triple de résonances, de sorte que chaque résonance peut étre indicée 
sous la forme d’un triplet de nombres entiers (n, p, g). On peut établir une relation entre 
les indices et les noeuds d’oscillation 4 la surface du cylindre, et ainsi déterminer expéri- 
mentalement ces indices. | 

Comme il n’est pas possible de résoudre de facon rigoureuse le probleme des valeurs 
aux limites dans le cas d’un cylindre de longueur indéfinie, il faut s’attendre a trouver 
certains écarts dans la vérification expérimentale de la théorie. Toutefois, dans le cas 
particulier des résonances a composante axiale nulle, on peut limiter l’erreur par com- 
paraison ayec les valeurs des résonances correspondantes d’une plaque circulaire mince. 


1. Einleitung 


Fir die Vermessung der elastischen Konstanten 
eines isotropen Festkérpers bieten sich seine elasti- 
schen Eigenschwingungen als bequemes und genaues 
MeBmittel an. Deren Frequenzen kénnen bei Anwen- 
dung moderner elektronischer Methoden mindestens 
ebenso genau wie die geometrischen Abmessungen 
des Probekérpers bestimmt werden. Als solcher wird 
der genauen Herstellbarkeit wegen oft ein kreis- 
zylindrischer Stab verwendet, und mit solchen Kreis- 
zylindern sollen sich die folgenden Ausfiihrungen 
ausschlieBlich befassen. 

Die elementare Theorie unterscheidet im wesent- 
lichen drei verschiedene Schwingungstypen: Tor- 


4) 
& 4 
“rs 


sions-, Dehnungs- und Biegeschwingungen. Die be- 
kannten Formeln fiir deren Phasengeschwindigkeiten 
und Eigenfrequenzen finden sich beispielsweise bei 
Rayxercu [1]. Indessen konstatierte dieser schon die 
Komplikationen, die durch die mit einer Langs- 
dilatation verbundene Querkontraktion verursacht 
werden, und er gab Korrektionsformeln hierfiir an. 
Schon die ersten, nach der Entwicklung der Ultra- 
schalltechnik moglich gewordenen Messungen von 
Boye und Sprovte [2], Rouricw [3], Scudneck [4] 
sowie Suear und Fockxe [5] zeigten aber die Un- 
zulanglichkeit dieser Korrektionsformeln fiir groBere 
Frequenzbereiche. Damit gewann eine bereits seit 
langerer Zeit bekannte theoretische Behandlungs- 
weise an Interesse, die von der durch die allgemeine 
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Elastizitatstheorie gelieferten Bewegungsgleichungen 
des elastisch-isotropen Festkérpers ausgeht. Diese, 
erstmals von PocuHammer [6] angegeben, gilt streng 
nur fiir Stabe unendlicher Lange, stellt aber min- 
destens qualitativ die auftretenden Erscheinungen 
befriedigend dar. Sie liefert fiir jeden Schwingungs- 
typ eine Gleichung, die sogenannte Frequenzenglei- 
chung, die die Phasengeschwindigkeit des Schwin- 
gungsvorganges mit seiner Frequenz verknipft. Diese 
Gleichungen haben im allgemeinen eine abzahlbare 
Schar von Losungen, die geometrische und elastische 
Daten des Stabes als Parameter enthalten. 


2. Die Herleitung der allgemeinen Frequenzen- 
gleichung 


Die oben erwahnte Bewegungsgleichung lautet 
fiir den Fall des  isotropen, kraftefreien Festkér- 
pers [7] 


i = (A+2G) grad div§—Grotrot’, (1) 


worin S den Verschiebungsvektor, A die Lamésche 
Elastizitatskonstante des Materials, G seinen Tor- 
sionsmodul und @ seine Dichte bedeutet. Fihrt man 
die Zeitdifferentiation unter der Annahme sinus- 
formig periodischer Zeitabhangigkeit der Kreisfre- 
quenz @ aus, so nimmt Gl. (1) bei Verwendung der 
Abktrzungen 


/ GC [Ae G 
on=| = n= | — (2) 
die Form 

0 = w? 8+ vy? grad divs — vq? rot rots (3) 


an. Setzt man in Zylinderkoordinaten 


S= (u,, Ue, Uz) , 


so lautet Gl. (3) ausgeschrieben 


> Oo divs a 1 3Q, “) 
—vy a = i 


0=2u, pie GEN" wag 
W* U;, + Vy a; = af a; 
nae 2 1 Odivs (Se 32," 
= W* Uy + Vy * - ree oe = 2 Z 
ees Oc T \ Oz Or 
a 50 divs Pr ouriat,) OQ; 
SE ee — oP (Go — ) 
mit 
, 1. O(r-usiieos Ou 
iui a! POR Lael 
lvs z a; oo z= Ze Ae 
oe 1 Ou, _ Que 
: r Op Oz 
Q, Sur Suz 
see Op’ 
o, 1 (leu) Bur : 
zl ie Or Op ; ) 
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wobei zwischen den drei-letzten Formeln von Gl. (5) 
die Identitat 


1 6@-2,) 
r Or 


1 025°, 9, 
r dp Oz 


=O) (6) 


besteht. Durch geeignetes Weiterdifferenzieren und 
Umrechnen erhalt man aus Gl. (4) die drei Glei- 
chungen 


(2.22.8 ee 
Or +r Or. Pr Og? 322 or) divs=0, (7) 
o? | ers Panes ees o? a? 
(3, Pore Bp toe yg) 20 @ 
2 3090 1 & o2 a? =o] 
(s+ or or Og? Oz? tee =) 2, 
2 02 
Basoces 4 


welche nach Einfiihren des Lésungsansatzes 


u,=U(r) cosnge 17% 


up=V(r) sinn ge i , (10) 
u,=W (r) cosnge iv , 
U0) alice ameter 
in die Besselschen Differentialgleichungen 
2 1 9 Pa sarees 
(sa+ 5 th" Javs <0, (11) 
2. VON Oa ae ears 
(ss++ +18 7)0,-0, (12) 
2 3090 v=] w 2 
Pail, Sey aD Cees Z ses 
& . ro sed bias Ja, ere 
mit h? = w? (2 3 
UL v~ 
1 ] (14) 
[2 eee 
“le as) 
tibergehen. 


Bezeichnen nun A, B,C beliebige Konstanten, so 
lauten die mit der Endlichkeitsforderung fiir r—0 
vertraglichen Losungen 


divs =AJ, (hr), (15) 

2,=Cla(kr) ’ (16) 
B _o| Jn (kr) 

0, =— ],(kr) —i>- Ca 


Diese sind in die Gl. (5) einzufiihren. Man erhilt 
dann vier Gleichungen fiir die drei Komponenten 
U, V, W, aber von diesen sind wegen Gl. (6) nur 
drei unabhingig. Bildet man also etwa aus der 1., 
2. und 4. Formel von Gl. (5) ein unabhangiges 
System, so ist der Lésungsansatz hierfiir 
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A BIn(hr) |B BIn(kr) | C mInlker) 


w/a? Or? /g2 Or ke r ; 
E A’ nJn(hr) Z B Jn (kr) & OJn (kr) 18) 
~~ @?/yy? r wo? vy? r ke cms 2) 
Mead al PG) re Be ke 
Per oto? % In(hr) +i AABF OSB Jn(kr) . 


Die hierin enthaltenen Konstanten A’, B’, C’ sind durch das Randwertproblem zu bestimmen. 
Dieses Randwertproblem wird durch die Forderung des Verschwindens samtlicher Spannungen auf der 
Zylinderoberflache formuliert. Ist a der Radius des Zylinders und L seine Lange, so lautet es 


Pry = AdivS+26 ot =0 fir r=a, yep ha Ze = 0 fir z=0,L, 

1 Ou, 6) & s Que 1 Ou, 5 
Pr = 6|1 setts, (}]= 0 fiir r=a, Peo = (5 as 3 = O00 -tur 2050, (19) 
P, = o( St +S) =Q ftir r=a; P,, = Adivs426 2! == Oe shins z= OF 5 

Oz Or Oz : 


Zunichst seien die durch die ersten drei Formeln von Gl. (19) gestellten Bedingungen behandelt: Setzt 
man die durch Gl. (10) mit Gl. (18) gegebenen Losungen darin ein, so erhalt man ein homogenes lineares 
Gleichungssystem fiir die Konstanten A, B’, C’, das sich nur erfiillen 1a8t, wenn die Determinante der 
Koeffizienten verschwindet. Mit den Bezeichnungen 


2 
ee ra ge On 20) 
nimmt diese die Form 
Reed rig—t) n?—1—27(22—1) 2(n*— 1) [O, (ka) —n] —1?(24—1) 
O,(ha)—n—1 ©,(ka)—n-1 One 2) OF (ba) — n|—72(2¢—1) =0 (21) 
0, (ha) —n — (z—1)[O, (ka) —-n] n é 


an. Sie ist identisch mit der von Hupson [8] angegebenen Frequenzengleichung. Ist Gl. (21) erfiillt, so 
kann man aus zweien der drei ersten Formeln von Gl. (19) die Konstantenverhaltnisse A4’/C’ und B’/C’ 
ausrechnen und erhalt endlich 


u=c"|-4” OIn(hr) py OJn (kr) +n) 
ar 8r r 
Vee [4’n Sn(hr) _ pn deke) _ Shek) ‘ (22) 
T tT; T 
WC" [ A yh) 4B" = In(k | me 
) oolv *: 
mit |2(n?—1)[O, (ka) —n] —2(2x—1) n?—-1—72(2x—-1) | : 
git 12m =21On(ka)—nJ—-P(22—-1)  On(ka)—n—-1 | : 
~ 2n n®?—1—7?(z—1)  n?-1—72 (22-1) | . { 
|On(ha) —n—1 0, (ka) —n—1 
| (23) , 
a L's (2 — 1) 2(n?—1)[O,(ka) —n] —?(2¢—1) | ‘ 
pr} |On(ha)—n-1  2n*-2[0,(ka) —n]—72(22-1) | ‘" 


n?—1—r(z—1) n?-1—71?2(22—1) | ‘ | Be 
O,(ha)—n—-1 @O,(ka)—n—-1 | 7 


ya 


Geht man mit den nunmehr bestimmten Ausdriicken fiir u,, u,, uw, in die letzten drei Formeln von 
Gl. (19) ein, so erweist sich das vorliegende Randwertproblem als tiberbestimmt, indem das gleichzeitige 
Verschwinden von cos @(L/v) und sin w(L/v) gefordert wird. Hierdurch wird der Naherungscharakter ¢. 
der so gewonnenen Lésungsdarstellung bedingt. 


é, 
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3. Diskussion der Frequenzengleichung 


Jeder ganzzahlige Wert n in Gl. (21) liefert eine 
Darstellung des Typs 


paces & =) oder a (eh) » (24) 
vr UL VT 


da ja vy/vy, eine Funktion der Querkontraktions- 
zahl o , namlich 


ae eee. 0<o<0,5 (25) 


ist. Wegen der Quasiperiodizitat der Besselfunk- 
tionen wird diese Darstellung (24) mehrdeutig, sie 
ist also fiir jeden festen Wert von o durch eine 
Schar abzahlbar vieler Kurvenaste wiederzugeben. 
Jeder dieser Aste kann durch eine ganze Zahl p in- 
diziert werden und liefert eine Serie von Resonan- 
zen, deren diskrete Frequenzwerte durch die Periodi- 
zitat der z-Abhangigkeit in den Randbedingungen 
fiir die Endflachen des Zylinders bestimmt sind. 
Diese wiederum konnen durch eine Laufzahl gq in- 
diziert werden. Eine jede Resonanz ist demnach 
durch drei Indizes bestimmt: 

1. n als Index fiir den Schwingungstyp, 

2. p als Index fiir die Schwingungsserie des Typs, 

3. q als Ordnungszahl innerhalb einer Serie. 

Fur n=O spaltet Gl. (21) auf in das Produkt der 
Gleichungen 

20 (ha) =r (2e—1)., (26) 


tO (ha) Oy(ka) +7? (2-1)? Oy (ka) + 
L7?2(2¢—1))0) (he) =05" (27) 


welche zweckmafigerweise mittels der Rekursions- 
formel 


ya So wi Ra 
O,(y) =2n rt (28) 
in @,(ka) =2 (29) 


und 


(x— 1)? 0, (ha) = (2 oe z—1) (0, (ke) =z], (30) 


ubergeftihrt werden. Gl. (30) ist die PocuHammer- 
sche Frequenzengleichung. Einsetzen in Gl. (22) 
zeigt, daB Gl. (29) den Torsionstyp 
See fir k=0 


Cr e(kt) EA ee 


Gi =AUe =| 


beschreibt, hingegen Gl. (30) den Longitudinaltyp 
U=C" [A J, (hr) —B’ J, (kr)), 
a (32) 


tA df ke? 
wc" | 4 i Jy (hr) +B” i ——Jg (kr) 
v olv 


Beim Torsionstyp liegen einfache Verhaltnisse vor: 
Gl. (29) hat abzahlbar viele und nur reelle Losun- 
gen y, =ka, deren jede entsprechend der Definition 


(14) einen Kurvenast 


2 
2¢="—3 +1 (33) 


liefert. 
Die numerische Auswertung von Gl. (29) ergibt 


Tabelle I 
¥1 = 0 
Upc) SRISts 
Ys = 8,417 
Tiss WECPAL, 


y,=90 bedeutet k=0, also v=vy. Dies ist die be- 
kannte, auch elementar zu berechnende Torsions- 
serie. Hiermit ist vy eine anschauliche Bedeutung 
unterlegt. 

Fur die hoheren y, entstehen p zur Achse kon- 
zentrische Knotenzylinder, wenn man die Achse 
selbst als solchen mitzahlt. 

Das Randwertproblem auf den Endflachen ist hier 
streng losbar und fihrt auf die Eigenwerte 


Asie Al) 
Wes ms ie qa 2 ae (34) 
Beim Longitudinaltyp, der durch Gl. (30) und 

(32) beschrieben wird, treten insofern Komplikatio- 

nen auf, als keine explizite Darstellung der einzelnen 

Serien mehr moglich ist, und diese tiberdies von der 

Querkontraktionszah] o abhangig werden. Eine Ta- 

belle des untersten Astes wurde von Bancrorr [9] 

mitgeteilt, eine solche der beiden nachsthéheren von 

Davies [10]. Hirer [11] und Hotpen [12] geben 

graphische Ubersichtsverfahren an. 

Folgende Spezialfalle seien noch erwahnt: 

a) Fir a0 wird v=vp=VE/o. 

Hierin ist F = 2 G(1 +0) der tibliche Elastizitats- 
modul. 

b) Fur v=vy reduziert sich Gl.(30) auf J, (ka) =0, 
was allgemein durch a— ~ erfiillt wird. Damit 
erhalt auch vj, seine anschauliche Bedeutung. 

c) Fur x=1 reduziert sich Gl. (30) auf 


8,(6)=1 mit =r. (35) 
Die numerische Auswertung ergibt 
Tabelle IT 
é, — 1,842 
Eo 5.086 
> 8.533 
z, =11,785 


Diese Werte liefern in jeder Serie einen von o 
unabhangigen Kontrollpunkt. 

Das Randwertproblem auf den Endflachen erweist 
sich hier als nicht mehr streng lésbar, wegen der 
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Phasenverschiebung der Radial- gegen die Axial- 
komponente. Wahlt man die Wellenlinge 4 wie in 
Gl. (34), so wird Pzz=0, und da wegen P-,=P,, 
und wegen des Bestehens von Gl. (30) jedenfalls 
P,,=0 fir r=a gilt, wird fiir kleine t auch P,, iiber 
die ganzen Endflachen hin klein bleiben. 

Der Typ der Biegeschwingungen wird durch n= 1 
charakterisiert. Von Hupson [8] liegt eine Tabelle 
der untersten Serie vor. Hohere Serien sind nach 
Wissen des Verfassers noch nicht behandelt. 

Das Randwertproblem auf den Endflachen ist 
hier nicht streng losbar. 

Typen mit n>1 sind in der Literatur noch nicht 
erwahnt, aber experimentell nachweisbar [13]. 


4. Die Radial-Azimutal-Schwingungen 


_ Rechnet man Gl. (21) auf den frequenzproportio- 
nalen Parameter 


(36) 


um, so 1aBt sich der Grenziibergang x—> © vollzie- 
hen und man erhilt, von einem Faktor [9,(t)) —7] 
abgesehen, 


T=ao/vp, tr=2e 


To" 
cpr =0. 
O,, (4%) —n—1 2n?—2[q(t9) —n] — 142 | (37) 


Einsetzen in Gl. (22) zeigt, da fiir die hierdurch 


n*— )— 


2(n?—1) [On (%) —n| —T "| 
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Gl. (39) eine rein radiale. Gl. (39) ist identisch mit 
der Arrryschen Frequenzengleichung [14]. 


Fiir den azimutalen Typ als einzigen lassen sich 
die Randbedingungen auf den Endflachen streng er- 
fiillen. 


Interessant ist nun der Vergleich mit den entspre- 
chenden Schwingungen einer diinnen Kreisplatte. 
Deren Frequenzengleichung [15] lautet 


2 | 
wa1 XG) [On(t%) —n]—t?| _ 9 
|On(B%) —n—1 2n®—2[O,(t) —n] — 1? | (40) 
lo 
: wae 
mit = oe (41) 


Im Falle des Zylinders verschwinden bei die- 
sen Schwingungen die Spannungen P,, und Pz,. 
Hingegen wird P,,~J,(at)). Mit wachsendem n 
und wachsender Laufzahl p, die hier die Rolle der 
Ordnungszahl der Resonanzen iibernimmt, werden 
nun die Endflachen von einem immer dichteren Netz 
von Knotenradien und zur Achse konzentrischen 
Knotenkreisen fiir P,, tiberzogen, zwischen denen 
aus Stetigkeitsgriinden die Amplitude keinen grofen 
Wert mehr annimmt. Man darf daher eine rasche 
Konvergenz der Werte fiir den Zylinder und die 
Scheibe erwarten, was die Tabelle III auch deutlich 
zeigt. 


Tabelle III 
T -Werte fiir den Zylinder 


n | o6=0 | o= 0,1 Fe 0,28 bes 018 ba Oslo = 0,5 
| | 

O(T) | 5,136 | 5,136 | 5,136 | 5,136 5,136 | 5,136 

O(L) | 2,604 | 2,885 | 3,204 B99 Tan uOsO1L9 yan 100 

1 | 2,474 2,583 | 2,698 2817 | 2,937 3,054 

2 | 2.336 2,341 2.344 2,348 2,351 2.354 


T -Werte fiir die Scheibe 


| 
n | = 0) | o=0,1 


T T ) 
| 


| | 


iota 0:20 por 0.3.4 ch O4aN) oOo 
O(T) | 5,136 5,136 | 5,136 5,136 5,136 | 5,136 
O(L) | 2,604 2,855 3,138 3,463 3,851 4,332 
1 2,473 2,573 2,659 2,734 2,800 2,857 
2 2,336 2,340 2,343 2,345 2,347 2,349 


bestimmten Frequenzen die Axialkomponenten ver- 
schwinden. 

Fir n=0 zerfallt Gl. (37) 
deren einer fiir sich gleich Null gesetzt, folgende 
Gleichung liefert: 


in zwei Faktoren, 


Oo(t%|) =—%?/2 oder O,;(t)) =2, (38) 
der andere 
O,(at%}) = —t)2/2 oder O,(at%)=2a7. (39) 


Gl. (38) beschreibt eine rein azimutale Schwingung, 


Diese Resonanzen lassen sich definitionsgema in 
das allgemeine Indizierungsschema einordnen, in- 
dem man ihnen die Laufzahl g=0 zuordnet. Dann 
kann allgemein n als die Zahl der die Achse ent- 
haltenden Knotenebenen der stationaren Amplituden- 
verteilung einer bestimmten Schwingungskomponente 
gemessen werden, q als die Zahl der Knotenebenen 
senkrecht zur Achse, p schlieBlich kann in bestimm- 
ten Spezialfallen als Zahl der zur Achse konzentri- 
schen Knotenzylinder bestimmt werden. 
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Dem Verfasser sei es erlaubt, Herrn Prof. Dr. 
H. O. Kneser fiir seine Unterstiitzung und Forde- 
rung dieser Arbeit an dieser Stelle zu danken. © 


(Eingegangen am 13. Juli 1956.) 
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MEASUREMENTS OF THERMAL AND ORDINARY DIFFUSION 


IN GAS MIXTURES (H.-He, H2-N, and H:-O,) 
BY USING ULTRASONICS 


by A. van Irrerpeex and J. Nrnoun 


Instituut voor Lage Temperaturen en Technische Physica, Leuven, Belgium 


Summary 

Measurements are carried out on thermal diffusion of H.He, H,—-N. and H,—O, mix- 
tures whereby the change of the concentration is determined from sound velocity measure- 
ments. Using the method formulated by Jonzs and Furry, measurements are made on the 
change of sound velocity caused by the change of concentration as a function of time and 
from this the mutual diffusion coefficient D,. can be determined for the different mixtures 
as a function of temperature. A good agreement is obtained with the values calculated 
from the equations of Hirscure.per. 


Sommaire 

On a mesuré la diffusion thermique de mélanges H,-He, H,—N. et H,—O,, en déter- 
minant la variation de concentration par des mesures de la vitesse du son. On a mesuré, 
en fonction du temps, par la méthode de Jones et Furry, la variation de vitesse du son 
due a un changement de concentration, ce qui a permis de déterminer le coefficient de 
diffusion mutuelle D,. pour les différents mélanges, en fonction de la température. On a 
obtenu un bon accord avec les valeurs calculées a partir des équations de Hirscuretper. 


Zusammenfassung 

Bei Untersuchungen iiber die Thermodiffusion in Gasgemischen (H,—He , H.-N» , H.—Os) 
wurde die Konzentrationsinderung durch Messung der Schallgeschwindigkeit bestimmt. 
Weiter wurde nach einer von Jones und Furry angegebenen Methode der gegenseitige Dif- 
fusionskoeffizient D;, der verschiedenen Gasgemische als Funktion der Temperatur aus der 
durch Konzentrationsanderungen verursachten zeitlichen Anderung der Schallgeschwindig- 
keit ermittelt. Man erhialt eine gute Ubereinstimmung mit den aus Hirscurenpers Glei- 
chungen berechneten Werten. 


1. Introduction 


A year ago we published [1] a few experimental 
data on thermal diffusion measurements in H,—He 
gas mixtures using ultrasonics. We described at that 
time the experimental technique, which consists in 
measuring the velocity of sound with an acoustical 
interferometer. These velocity measurements are 
used to determine the concentration of the compo- 
nents of the gas mixtures and also the change of 
concentration produced under the influence of ther- 


mal diffusion. Other investigators like Isps and Grew 
[2] make use of heat conductivity measurements to 
determine the change of concentration. 

In our previous paper we showed the possibility 
of separating with our experimental method the 
thermal diffusion coefficient from the ordinary dif- 
fusion coefficient. Independently Nerriey [3] pub- 
lished just before ours a paper on the same subject. 
The experimental technique consists in measuring 
the change of the velocity as a function of time, and 
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to calculate from this the change of the concentra- 
tion Ay as a function of time ¢. Corresponding with 
Ciark Jones and Furry’s [4] theory we have 


(Ay) = (47) 4-0(1—e7""). (1) 


The relaxation time t=C/D,. ; C being a constant 
(see further), which can be calculated from the di- 
mensions of the apparatus used and from the known 
experimental conditions. We have carried out meas- 
urements on the mixtures H,-N,, H,—-O,, H,—He, 
from which we calculated the ordinary diffusion co- 
efficient D,. for the mixtures in question and com- 
pared the experimental results with the values 
calculated from the equations given by Hirscuret- 


per [5]. 


2. Experimental apparatus and method 


In order to check the theory of Jones and Furry 
a new apparatus was built. The apparatus is con- 
ceived in such a manner that the conditions mention- 
ed in the theory are approached as closely as possible. 

In our previous apparatus the connecting tube 
between the cold volume (temperature 7,; boiling 
point of liquid oxygen or hydrogen) and the volume 


Fig. 1. Diagram of apparatus. 


at room temperature 7’, was not well defined. This 
was also the case for the ratio between the two vo- 
lumes. In our previous experiments the velocity 
of sound was also measured in the cold volume ly. 
However to determine the exact concentration we 
had some difficulty, because the velocity becomes a 
function of pressure at low temperatures. On the 
present apparatus the velocity is measured at room 
temperature. The apparatus (Fig. 1) is constructed 
of brass. It consists of two parts: V, and V,, V, 
contains a quartz crystal with a large diameter 5 cm, 
to have a great sensitivity. 

The frequency of the quartz crystal (+600 ke/s) 
is determined at each velocity measurement by means 
of a precision frequency meter. a is a tube in ger- 
man silver which connects the two volumes V, and 
V, (total length: 305 mm, inner diameter 3.9 mm). 
The other dimensions are indicated on the figure. 
The ratio V,/V calculated from the dimensions is 
equal to 0.62. 

After the apparatus has been evacuated and filled 
with hydrogen gas so that the different parts reach 
the desired temperature, the apparatus is evacuated 
again and filled with the mixture to be studied. 


The technique consists further in measuring the 
change of the velocity of sound in the mixture as a 
function of time at intervals of 10 minutes. From 
the curves obtained we calculate the instantaneous 
change of the concentration Ay in the volume /,. 

We used the equation 


= ; ———— (2 
Cys bey (1-71) My 71+ M2 (1—74) el 


panee. Var ey (UY) s 


in which the symbols have their ordinary signi- 
ficance. The total change of the concentration is 
obtained from 


Ay = Ay, + Aye. 


Ay, being the change of the concentration in the 
cold volume; Ay, is calculated from Ay, by means 
of the equation 


VTS ar: 
Yo = Ay ="=Ay, -0.62: — 
A) 2 Ay, V, ii V4 ( 62 (ip 
or Ay=kAy,. 


In our experiments the values of / used are: 
T,=290°K, T,=90°K, k=1.16; 
T,=290° K, T,=77°K, k=1.14; 
T,=290°K, T,=20°K, k=1.04. 


So finally graphs are drawn for Ay as a function 
of time. To these graphs we applied the equation (1) 
and calculated the relaxation time t. According to 
Jones and Furry’s theory the relaxation time 7 can 
be written 


at See, 
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N, No L Ts T, 

ES Ome er eee 3 
“aWeaING “a Dia Gots aie (3) 
N, and N, being the number of molecules effectively 
in the volumes V, and V,, n the molecular density 
in the volume V,. Equation (3) can also be written 


VV. (¥ VN Ait Ts 
= ee |i cat a ae 4 
TOWKT +5] De SUE = (4) 


Tv 3 
r, 
AT =T,—T,; L is the length of the connection 
tube, S its cross-section. 
From the values obtained for t and using equa- 
tion (3), the coefficient D,. corresponding to the 


average temperature can be calculated 
: Tors fs 
i Sees 
ToT, L; 
The measurements were carried out with two dif- 
ferent cross-sections of connecting tube, but there 
was a complete agreement. 


3. Experimental results 
3.1. H.He mixtures 


The measurements are carried out for mixtures 
with concentrations of approximately 50%, for the 
reason that for those mixtures the thermal diffusion 
effect is largest. 

In TableI are given the experimental results 
together with the change of the concentration cal- 
culated from the measurements. - 

Measurements were also carried out for the same 
mixture but with 7,=17.1°K for which we ob- 
tained for t (1 atm) =9.52 «10's. 

In Fig. 2 are represented as an example the meas- 
urements corresponding to the fourth series of 
measurements. 

In Fig. 2 three curves are brought together: 
curve 1 represents the change of W as a function 
of t, curve 2 represents y or Ay as a function of ¢, 
and finally curve 3 gives —In[1— (Ay/Ay..)] as a 
function of ¢. From Fig. 2 we see that curve 3 is a 
straight line, which agrees with Jones and Furry’s 
theory. 


3.2. H,-O, and H.-N, mixtures 


The same kind of measurements were also carried 
out on mixtures H,—O, and H,—N,.. The results of 
those measurements are given in Table II. 

In Fig. 3 we have graphically drawn the values 
of —In[1— (Ay/Ay..)] for the H,—N, mixtures as 
a function of t. We see that again straight lines are 
obtained. The same is the case for the H,—-O, mix- 
tures. 7 
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Table I 
Thermal diffusion in H,—He mixtures 
_(a) yo =51.0% He; T,=291° K, T,=—90° K; 
T=153° K; p=32.2 cm; frequency 605.00 ke/s 


i emcee sire =mn(1- | | z-10%s | r- 104s 
es ae ae | % | Ay ~/ | (32.2 em) | (1 atm) 

0 | 1097.5 | 20 1.62! 0.428 

6 | 1098.7 | 30. 1.85 0.509 

13. 1099.5 | 40 2.32 0.691 

20 | 1101.0 | 50 2.75| 0.875 

27 1101.2 | 60 3.07 1.087 

38 1102.4 | 703.40 1.21 

a as © | 3.54 | 8.36 
134 1107.7 | 
193 1108.3 
226 1108.4 ) 


(b) Yo = 54.5% He; T,=290° K, T.=77.3° K; 
T=137° K; p=43.7 cm; frequency 605.00 ke/s 


g tv) | «4088 t- 10s 
tr | 30 em) (Latm) 
0 1085.5 | 20 0.99 0.230 | 
6 1086.6 | 40 1.72 0.456 | 
14 1086.9 | 60 2.30) 0.672 | 
27 | 1088.5 | 80 | 2.85 0.932 | 
41 1090.0 }100 3.28, 1.197 | 
61 | 1092.2 | co | 4.7. | | 5.08 | 8.84 
99 1093.5 ) 
230 1096.4 
253 1096.6 | 
300 1096.8 | 


(c) ¥p9 =37.6% He; T,=291° K, T,=77.3° K; 
T =137° K; p=52.0 cm; frequency 604.72 ke/s 


: “A t-10%8s | t+ 107s 


(1 atm) 


min m/s 


6 | 1137.9 | 20 | 0.78 0.200 . 
22 | 1139.2] 40 | 1.42 0.401 . 
40 1141.4] 60 | 1.92 0.614 
62 1143.0 | 80 | 2.36 0.796 
97 | 1144.6 |100 | 2.71] 0.995 
138 | 1147.1 150 | 3.42} 1.586 
280 | 1149.5 | co | 4.3 6.00 | 8,77 
| 


(d) 7=55% He; 7,=289°K, T,=20.3°K; 
T=58° K; p=42.7 cm; frequency 604.80 ke/s 


t+ 1048 


= _ Ay t-10%s 
in(1 ) | (42.7 em) 


0 | 1082.6 | 20 | 2.29 0.263 
6 1084.9 | 40 3.80 0.484 
17: 1086.4 | 60 . 4.97) 0.697 
31 | 1090.2 | 80 | 5.94} 0.916 
51 1093.8 }100 6.72 1.136 
78 | 1097.7 1150 1.688 
111 | 1100.7 | @o | 9.9 _ 9.48 
128 | 1101.9 
271 | 1107.8 
350 1108.1 
402 1108.6 


g : \ : - ° . 
a 
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0 200 400 min 


i= 


Fig. 2. Thermal diffusion in H.—He mixtures 
(T,=20° K, 79 =55% He). 


Table II 
Thermal diffusion in H,-O, and H,—N, mixtures 
(a) Hy—O,; 9 =53.2% Os; T; = 291° K, T,=90.1° ke 
T=153°K; p=19.2 cm; frequency 604.80 ke/s 


| 
| 
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(c) Hy—-Ns° yp =51.6% Ny; T,=291.5°K, T,=90° K; 
T=153° K; p=24.4cm; frequency 608.55 ke/s 


ieee a 


| | 
| ae | —In(1— v +107 
| (24.4 cm) 


Zl OS 
(1 atm) 


(d) Hy—-No; 79 =51.0% N,; T,=288.2° K, T,=90° K; 
T=153° K; p=53.3 cm; frequency 608.59 ke/s 


t-10%s 
(1 atm) 


| t-10%s 
(53.3 em) 


27 | 
62 


90 

116 : 11.82 | 16.85 
570 | 492.2 | | 

486 | 493.2 | | | 

600 | 493.2 | 

(e) Hs-Ns; yp =53.5% No; T,=292°K, T,=90°K; 


T=153° K; p=28.0 cm; frequency 604.80 kc/s 


| z-10-%s | t- 10-8 
| (28 cm) 


; a atm) 


t-10%s | t-10%s 
(19.2 em) | (1 atm) 


(b) Hy—Og; Yo =62.1% Oz; T,=290° Ks T,=80.1° K; 
T=142° K; p=15.5 cm; frequency 604.77 ke/s 


| | | 
—in(1- =~") | r-10%s | r- 10s 
| AV en | (15.5 em) |) (1 atm) 


lense 
7 407.6 


22 408.5 

34 | 411.0 

57 | 412.1 

78 | 413.7 | o [5.1 3.72 | 18.20 
110 416.4 | 

195 | 417.5 

210 418.4 


158 | 6.24 | 16.94 

4 = 

3 

1 

0 100 : 200 300 min 700 

Fig. 3. H.-N, mixtures; yy=+50%; 

curve 1: p=24.4cem, T,=90.1° K; 
curve 2: p=23.3 cm, T,=77.3° K; 
curve 3: p=53.3 cm, Tz=90.1°K. 
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Using values of the relaxation time t of Tables I 
and II and equation (2) we calculated the diffusion 
coefficient D,,. On the other hand we also calculated 
this coefficient by means of the equations and the 
empirical data given in Hirscuretper’s book [5]. 
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As up to now no direct measurements on dif- 
fusion coefficients at low temperatures are known, 
we can conclude from Table III that this indirect 
method works fairly well. At the highest tempera- 
tures the agreement between the theoretical values 


Table II 
; 7, Room temp. | ah Tem Oe Dyp exp. Dy» theor. 
Mixtures re Tae = | 5 
oK K s em?/s em?/s 
| 

H,—He (50%) eM weed eat Gir Mail eames Rien elt bau os 0.518 

77.3 137} 8.84 | 0.343 0.429 

20.3 58 9.48 __| 0.093 - 0.098 

digi! 52 9.52 0.076 0.078 

H,—He (35%) 77.3 137 8.77 0.346 0.429 

-H,_N, (60%) 90.1 153 17.0 0.207 0.237 

90.1 153 16.8 0.209 0.237 

90.1 153 16.9 0.208 0.237 

Tiles 137 19.6 0.173 0.194 

H,—0O, (50%) 90:19 >") eels 16.4 | 0.215 0.250 

H,—O, (60%) | 80.1 | 142 18.2 0.174 0.218 


D2 =0.0026280 


D,» diffusion coefficient in cm?/s , 

p pressure in atmospheres , 

7 temperature in °K, 

Typ = kT/e 9, 

M,,M, molecular weights of species 1 and 2, 

O19 €49/k molecular potential energy parameters 
characteristic of the 1—2 interaction in A 
and °K respectively, 

Q°-) the well-known integrals for calculating the 
transport coefficients for the Sutherland poten- 
tial. 

The results obtained are given in Table III. 


and the experimental values is not so good as at 
liquid hydrogen temperatures where the change of 
the concentration is about twice the change at liquid 
oxygen temperatures. 


We take the opportunity of expressing our thanks 
to the Belgian Ministry of Education for financial 


aid during these experiments. 
(Received May 3rd, 1956.) 
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MUSIKINSTRUMENTENSTIMMUNGEN UND TONSYSTEME 


Beispiele: Klavier, Fliigel, Akkordeon, Harmonium, Streichinstrumente 


von H. Meren 


Forschungsinstitut fiir Musikinstrumentenbau des Deutschen Amtes fiir Material- und 
Warenpriifung, Zwota iiber Klingenthal/Sachsen 


Zusammenfassung 


Die Stimmung yon Klavieren und Fliigeln einerseits, Akkordeons und Harmoniums 
anderseits, entspricht verschiedenen temperierten Tonsystemen. Fiir die letzteren gilt das 
bisher als alleinherrschend angesehene System der Einteilung der Oktave mit dem Fre- 


LZ 
quenzverhaltnis 1 : 2 in 12 Halbténe der Gréfe 2. Fiir gut gestimmte Klaviere und Fli- 
gel wird jedoch in der Praxis eine Temperierung mit um etwa 3 cent pro Oktave erweiter- 


ten Intervallen beniitzt, HalbtongroBe V 2,0017. Gutgestimmte Klaviere und Fliigel weisen 
deshalb nicht temperierte, sondern annahernd reine (natiirlich-harmonische) Quinten auf. 
Da auch Streichinstrumente in reinen Quinten gestimmt werden, ergeben sich fiir Klaviere 
und Fliigel einerseits, Streichinstrumente (Geigen, Bratschen, Celli) anderseits gleichlau- 
fende Stimmungskurven. Ursachen und Bedeutung dieser Ubereinstimmung und technische 
Anwendungsmoglichkeiten der Stimmungskurven werden besprochen. 


Summary 


The use of differently tempered scales for pianos and accordions or harmoniums must 
be accepted. In the latter the octave is divided into 12 semitones of value 2'% but in the 
former the practice uses octaves 3 cents wider, making the semitone 2.00177" so that 
the fifths are almost pure untempered. As bowed-string instruments also have pure fifths 
they and the piano have similar tuning curves. The causes and results of this agreement 
are discussed as well as its technical possibilities. 


Sommaire 


L’emploi de différents modes tempérés pour les pianos droits et 4 queue d’une part, les 
accordéons et les harmoniums d’autre part, doit étre considéré comme bien établi. Dans le 
cas de ces derniers instruments, on utilise le systeme considéré jusqu’a présent comme 
unique, dans lequel l’octave est partagée dans le rapport de fréquences 1 : 2 en 12 demi- 


12 ,— 
tons de valeur V2. Dans le cas des pianos droits et 4 queue bien accordés, cependant, on 
utilise en pratique un tempérament avec intervalles élargis d’environ 3 cents par octave: 


valeur du demi-ton /2,0017. Les pianos droits et A queue bien accordés présentent donc 
des quintes non tempérées mais a peu pres pures (harmoniquement naturelles). Comme les 
instruments a cordes sont aussi accordés par quintes pures, les pianos droits et 4 queue 
d’une part, les instruments a cordes d’autre part (violons, altos, violoncelles) ont des cour- 
bes d’accord analogues. On indique les causes et signification de cette concordance, ainsi 
que les possibilités techniques d’utilisation de ces courbes d’accord. 


1. Klavierstimmung und Tonsysteme 


Seitdem es méglich ist, Frequenzen relativ schnell 
und genau zu messen, wurde aus verschiedenen 
wichtigen Griinden versucht, tiber Klavierstimmun- 
gen, insbesondere tber die beste Klavierstimmung, 
Klarheit zu gewinnen. Gearbeitet haben dariiber 
Griirzmacuer und Lorrermoser [1], O. L. Raits- 
pack [2], R. W. Youne [3], D. W. Sraurrer [4], 
H. Meret [5]. Um gute Vergleichsmoglichkeiten 
auch mit unseren Ergebnissen zu haben, wurden 
Kurven aus den Arbeiten der obengenannten Auto- 
ren in einheitlichem Mafstab in Abb. 1 wiedergege- 
ben und unserer Darstellungsweise angeglichen. Aus 
der Arbeit von Griirzmacuer und Lorrermoser 
wurde dabei, unserer Zielsetzung entsprechend, nur 


die beste der dort gezeigten drei Stimmungen aus- 
gewahlt, die jedoch noch nicht wirklich gut ist. 
Youne weist selbst darauf hin, daB die von ihm ver- 
offentlichte Kurve nur die Steigung der eigentlichen 
Stimmungskurve darstellt, nicht die Stimmung selbst. 
Sraurrer bringt eine Kurve, die von zwei Klavieren 
stammt, ohne dafi nahere Angaben tber die Ge- 
winnung gemacht werden. Sie weisen je Taste nur 
einen Mefiwert auf und lassen die Gite der Stim- 
mung nicht geniigend beurteilen. Sie muf erkenn- 
bar sein, will man sich ein Bild iiber die besten 
Klavierstimmungen machen. So sind wir zu Ver- 
gleichszwecken im wesentlichen auf die Kurve von 
Ramspack angewiesen, werden aber die anderen 
nicht vernachlassigen. Wir selbst haben aus dem 
Priifbetrieb der letzten Jahre umfangreiches Mate- 
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a) Gritzmacuer und Lorrermoser, b) RatisBack, 


c) Youne, d) Sraurrer. 


rial (etwa 150 Stimmungskurven) vorliegen, das 
sich auf zahlreiche und auch auf Klavierbaufirmen 
von Weltruf erstreckt. 

Die Abb. 1 und 2 bringen Beispiele. Waagerecht 
sind die Tonbezeichnungen der Klaviertasten, senk- 
recht die Abweichungen von der Werckmeister- 
Temperierung aufgetragen. Die von Werckmeister 
(1691) eingefiihrte, schon vorher bekannte Tempe- 
rierung — diejenige, die wir heute schlechtweg als 
temperierte Stimmung bezeichnen: Einteilung der 
Oktave mit dem Frequenzverhaltnis 1:2 in 12 


Halbton-Intervalle der GroBe V2 — hbedeutet also 
in unserer Darstellung eine waagerechte Gerade, die 
durch Null geht. Andere Temperierungen lassen die 
Gerade von links nach rechts fallen bzw. steigen, 


wenn die Halbténe <2 bzw. > V2 sind. Mehrere 
ibereinanderliegende MeBwerte, vgl. Abb. 1b und 
2a—k, gehoren zu den mehr oder weniger gleichhoch 
gestimmten Saiten einer Taste, zum gleichen ,,Chor“, 
Gemessen wurde mit dem von uns _ entwickelten 
elektro-optischen Stimmgerat [5]. 


Die Abb. 2 enthalt Stimmungen folgender Instru- 
mente: 
a) Gute Stimmung eines im Rundfunkbetrieb beniitz- 
ten Fliigels (nicht vom Herstellerbetrieb gestimmt). 


b) Beispiel der guten Stimmung eines Fliigels unmit- 
telbar vor und nach einem Klavierkonzert yon Frau 
Prof. Etry Ney. Die Kurven geben gleichzeitig ein 
Bild der sehr geringfiigigen Anderungen, die ein 
Fliigel durch den Konzertgebrauch erleidet, wenn 
er haltbar gestimmt ist, und der guten Reproduzier- 
barkeit der MeBergebnisse. 

c) Beispiel einer etwas weniger guten Klavierstim- 
mung. 

d) Sehr gute Klavierstimmung. 

e) Sehr gute, von einem Blinden ausgefiihrte Klavier- 
stimmung. : 

f) Etwas zu stark ansteigende, aber sonst als gut an- 
zusehende Stimmungskurve eines Klaviers, vgl. da- 
zu Abb. le und 1d. 

g) Noch, gute Fliigelstimmung, aber mit (ungewollten 
und unerwiinschten) GesetzmaSigkeiten auSerhalb 
einer Temperierung: Héherstimmung yon dis, e. 

h) Gute Fliigelstimmung. 

i) Gute Fligelstimmung. 

k) Gute Fliigelstimmung. 

1) Streichinstrumentenstimmung. 

Der Ordinatenwert O bedeutet gleichzeitig fiir das 
a’ die Normstimmtonfrequenz 440 Hz. Von uns werden 
aber nicht Abweichungen des at von 440 Hz als Ab- 
weichungen vom Normstimmton angesehen; denn sie 
kénnten zufillig sein und daher kein wahres Bild der 
Hohe der gesamten Stimmung des Instrumentes geben, 
auf die es doch ankommt. Wir betrachten yielmehr die 
Abweichung der Stimmungsgeraden — oder, wenn, be- 
sonders bei Blasinstrumenten, keine solche vorliegt: 
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Abb. 2. Gute Klavier- und Fliigelstimmungen eigener Messung. Instrumente der Firmen C. Bechstein, Berlin; 


Julius Bliithner, Leipzig; VEB Pianofariefabrik Kisenberg/Thiir.; August Férster, Lobau; Grotian-Stein- 


weg, Braunschweig; Sichs. Pianofortefabrik VEB Seifhennersdorf/L.; Steinway & sons, Hamburg. 
Die Reihenfolge der Kurven stimmt nicht mit der Reihenfolge der Firmen. iiberein. Weitere Stimmngs 
kuryen, auch yon schlecht gestimmten Instrumenten, unter [5]. 


der Stimmungsmittellinie — bei a! von 440 Hz als Ab- 
weichung vom Normstimmton (vgl. Abb. 2a). So ent- 
sprechen z. B. die Stimmungen in Abb. 2b, 2e, 2f, 2g 
der Normung a! =440 Hz. 


Die Beurteilung der Giite der von uns gemessenen 
Stimmungen fuBte urspriinglich auf der musikalisch- 
kiinstlerischen Beurteilung seitens unserer Gutach- 
terausschiisse. Spater konnte sie sich auf die physi- 
kalischen MeSergebnisse allein stiitzen. Es hatte sich 
z. B. oft gezeigt: Gehdren zu einer Taste mehrere 


Saiten und fallen die davon herriihrenden Frequenz- 
mefwerte nicht zusammen (Toleranz grofer als 
etwa 2 cent), so wird der Klang dieser Taste, durch 
Schwebungen bedingt, nicht als ruhig empfunden, 
er ,,quirlt“, ,,flackert*, ,,.wimmert“ oder wie man es 
sonst nennen mag, geniigte jedenfalls nicht den 
kinstlerischen Anspriichen unserer Gutachteraus- 
schiisse. Die musikalische Bedeutung der physikali- 
schen MeBergebnisse der Frequenzunterschiede zwi- 
schen den verschiedenen Saiten einer Taste, die 
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Giite der Stimmung in dieser Richtung, entspricht 
etwa der Grofe dieser Unterschiede. Am besten ist 
nach den uns vorliegenden Urteilen eine vollige 
Ubereinstimmung der Saitenfrequenzen. Fir die 
musikalische Bedeutung der Unterschiede liegt je- 
doch eine Art Schwellwert von etwa 1 cent vor. Noch 
geringere Unterschiede zwischen den Frequenzen der 
Saiten einer Taste bringen keinen bemerkenswerten 
musikalischen Gewinn mehr. Werden die Unter- 
schiede groer als etwa 2 cent, beginnen sie unan- 
genehm zu wirken. 


Anfangs waren wir wegen verschiedener Hinweise 
in der Literatur und auf Grund von Behauptungen 
mancher Stimmer unsicher in der Beantwortung der 
Frage, ob die beste Klavierstimmung tatsachlich 
moglichst genau temperiert sein soll, oder ob nicht 
doch gewisse Abweichungen von der temperierten 
Stimmung zweckmaBig sind, etwa in Richtung einer 
mitteltonigen Stimmung, d.h. der Bevorzugung ge- 
wisser Ofters beniitzter Tonarten, vgl. [1]. Wir ver- 
wendeten eine beachtliche Zeit darauf, in unseren 
Klavier- und Fligel-Stimmungskurven, die manch- 
mal erhebliche Abweichungen von einer temperier- 
ten Stimmung zeigen, verntnftige Gesetzmafigkei- 
ten und schlieBlich solche irgendwelcher Art auBer- 
halb einer Temperierung zu entdecken (vgl. auch 
[1]. Die Priifung unseres gesamten Materials, ins- 
besondere auch der jeweils von einem Stimmer vor- 
liegenden verschiedenen Stimmungen und unsere mit 
zahlreichen guten Stimmern durchgefihrten Ver- 
suche, zeigte jedoch in keinem Fall gewollte und 
daher wenigstens bei demselben Stimmer immer 
wieder vorkommende Abweichungen  gegeniiber 
einer temperierten in Richtung einer mitteltonigen 
Stimmung. Fast zwangslaufig weisen tbrigens die 
gut temperierten Stimmungen auch geringste Unter- 
schiede innerhalb eines Chores auf. 


Betrachten wir die Stimmungskurven in ihrem 
Verlauf von Klang zu Klang, so sehen wir vor allem 
einen Anstieg von links nach rechts, vgl. Abb. 1 und 
Abb. 2, wobei die Kurven fur die kiinstlerisch am 
besten beurteilten Stimmungen weitgehend eine Ge- 
rade darstellen. Die Halbtone sind also etwas gro- 


12 
er, als einem Faktor y2 entspricht und die Gerade 
selbst zeigt, da} eine Temperierung vorliegt. Die 
bei der Temperierung von unseren besten Stimmern 


sezeigte Leistung (Abweichungen im _alleemei- 
g g g g g 


nen St ; cent) ist tberraschend hoch [5]. Manch- 
mal nahert sich die Stimmung in den mittleren Ok- 
taven einer waagerechten Geraden (Abb. la bis 
d, Abb. 2h bis k), manchmal finden wir, dai bei 
den tiefen Klangen die Stimmung weniger absinkt, 
bei hohen weniger steigt als der schragen Geraden 
entspricht (Abb. 2b, c), manchmal ist es umge- 
kehrt (Abb. le, d, Abb. 2f). Aber das sind 
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Effekte (zweiter Ordnung), die bei unseren Betrach- 
tungen keine Rolle spielen, noch dazu, da wir sie 
bei unseren zahlreichen Messungen nicht als typisch 
feststellen konnten. Vom Stimmer werden sie beim 
Spielen von besonders ftir die Temperierung geeig- 
neten Akkorden gehoért, kaum beim Spiel von Kla- 
viermusik. Typisch ist dagegen das allgemein fest- 
zustellende Steigen der Stimmungskurve '. 

Es bedeutet, da8 neben der bisher allgemein als 
alleinherrschend angesehenen Werckmeister-Tempe- 
rierung bei einem unserer bedeutsamsten Instru- 
mente, dem Klavier oder Fligel, eine andere Tem- 
perierung besteht. Diese offenbar auf praktische 
Bediirfnisse zuriickgehende neue Temperierung ver- 
groBert sogar das Frequenzverhaltnis der Oktave, 
des konservativsten Intervalls, das sowohl im pytha- 
goreischen und natiirlich-harmonischen als auch im 
bisherigen temperierten Tonsystem 1:2 betragt. 
Zwangslaufig erhebt sich hier die Frage, ob auch 
andere Instrumente diese besondere Art der Tempe- 
rierung aufweisen. 


2. Akkordeon- und Harmoniumstimmung und 
Tonsystem 


Wir haben in den letzten Jahren an hunderten, 
zum wesentlichen Teil ganz hervorragend gestimm- 
ten Akkordeons, mehreren Harmoniums und eini- 
gen Orgeln, darunter beriihmten und von Organi- 
sten mit Weltruf bentitzten, die Stimmung gemessen. 
In keinem einzigen Falle haben wir bei diesen In- 
strumenten eine ansteigende (oder fallende) Stim- 
mungsgerade gefunden, sondern immer nur die 
mehr oder minder gegliickte Angleichung an die 
waagerechte Gerade unserer Darstellung, an die 
Werckmeister-Temperierung. Die Abb. 3 zeigt ein 
Beispiel einer guten Akkordeon- und Harmonium- 
stimmung (Akkordeonstimmungen  verschiedener 
Giite unter [5]). An unseren Stimmungskurven von 
Klavieren finden wir dagegen nur ausnahmsweise, 
d. h. in etwa 5 bis 10% aller gemessenen Stimmun- 
gen, eine waagerecht liegende Kurve. Bei derart um- 
fassendem Material und solchen Ergebnissen darf 
man die Beniitzung verschiedener temperierter Ton- 
systeme bei Klavieren und Fliigeln einerseits, Ak- 
kordeons und Harmoniums anderseits als festste- 
hend ansehen. 


1 Die dieser Arbeit zugrunde liegenden Klavierstim- 
mungen gehen auf verschiedene Methoden des Klavier- 
stimmens zuriick, und darauf sind im wesentlichen die 
oben erwahnten Unterschiede in der Stimmung, als Ef- 
fekte zweiter Ordnung bezeichnet, zuriickzufiihren. Uber 
die Richtigkeit dieser oder jener Stimmethode gibt es 
Meinungsyerschiedenheiten. Auf Grund der yorliegen- 
den Ergebnisse diirften sie sich, von der Praxis aus ge- 
sehen, als unbegriindet erweisen. 
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Abb. 3. Beispiel guter Akkordeon- und Harmoniumstimmung. 


3. Ursachen der Beniitzung verschiedener 
temperierter Tonsysteme 


Unter den geschilderten Umstaénden muf eine be- 
stimmte Ursache fiir die verschiedenen Stimmungs- 
arten der erwahnten Musikinstrumente vorliegen. 

Wir priiften die Meinung vieler Klavierstimmer 
und mancher Pianisten, nach der die erweiterten 
Intervalle dem Klang eine grofSere Brillanz verlei- 
hen, selbst nach, und kamen zu gleicher Auffassung. 
Eine solche Wirkung liegt auBerdem im Zuge der 
Entwicklung des Klavierklanges in den _ letzten 
200 Jahren. Etwa seit der Zeit Mozarts wurde die 
Erhohung der Kraft und Brillanz des Klavierklan- 
ges u. a. durch Erhohung der Saitenspannung tech- 
nisch in starkem Mafe verwirklicht und ist durch 
einen Vergleich von originalgetreu besaiteten Kla- 
vieren aus der Zeit Mozarts mit heutigen deutlich 
zu héren, z. B. im Deutschen Museum in Miinchen. 
Der Wunsch nach gro8erer Brillanz diirfte aber nicht 
als alleinige Ursache fiir die erweiterten Intervalle 
anzusehen sein, wahrscheinlich nicht einmal als die 
wichtigste; denn die mit der angegebenen Erweite- 
rung der Intervalle verbundene Erhéhung der Bril- 
Janz ist subjektiv nicht sehr bedeutsam. Dem heu- 
tigen Klavierbau stiinden wirksamere Mittel zur 
Verfiigung. 


Um weitere Ursachen zu finden, liefen wir unter 
gleichklingenden Akkordeons einige mit erweiterten 
Intervallen stimmen und die musikalische Wirkung 
durch Sachyerstandigengremien beurteilen. Die Er- 
héhung der Brillanz wurde auch-beim Akkordeon 
festgestellt und zwar auch hier als nicht sehr beacht- 
lich, wenn auch horbar. Die Sachverstindigen kamen 
jedoch nicht zu einer allgemein positiven Einstel- 
lung gegeniiber der Anwendung der besten Klavier- 
stimmung fiir das Akkordeon, sondern mehrere ga- 
ben der Werckmeister-Temperierung den Vorzug. 


Das mag darin begriindet liegen — was fiir Har- 
monium und Orgel wenigstens teilweise ebenfalls 
gelten konnte —, daf diese Instrumente Zungen- 
stimmen mit ihrem starken Gehalt an hohen Teil- 
tonen besitzen und eine Erhohung der Brillanz da- 
her gar nicht brauchen. AuBerdem scheint es bei 
langer andauernden Klangen und Akkorden — die 
auf Akkordeon, Harmonium und Orgel im Gegen- 
satz zum Klavier moglich sind — so zu sein, da 
sich bei ansteigender Stimmungsgeraden Schwebun- 
gen und Rauhigkeiten des Klanges unangenehm be- 
unerkbar machen. Sie lassen sich durch das Auftre- 
ten zusatzlicher Interferenzen zwischen Teiltonen, 
besonders Oktaven, als Folge cer steigenden Stim- 
mungsgeraden erklaren. Damit sind die auffalligen 
Unterschiede in der Lage der Stimmungsgeraden 
zwischen Klavieren einerseits, Akkordeons, Harmo- 
niums und Orgeln anderseits in ihren Ursachen eini- 
germafen verstandlich. Eine quantitative Betrach- 
tung der Steigung der Stimmungsgeraden fihrt wei- 
ter. 

Die Steigung der Stimmung der besten Klaviere 
betragt ungefahr 3 cents pro Oktave (siehe dazu die 
vergleichsweise in den Abb. 1 und 2 eingezeichne- 
ten schragen Geraden, die alle eine Steigung von 
3cents pro Oktave aufweisen). Diese Steigung be- 
deutet, da die Oktaven um 3 cents grofer sind als 
1:2, also =1: 2.001735, daB die Halbtonerweite- 


rung j5 cent betrigt und die Erweiterung einer 


Quinte (7 Halbténe) re cents, also fast 2 cents aus- 
macht. Da die reine (natiirlich-harmonische) Quinte 
mit dem Frequenzverhaltnis 2:3 auch um 2 cents 
groBer ist als die nach Werckmeister temperierte, 
bedeutet das, daB bei der hervorragenden Klavier- 
stimmung die temperierten Intervalle so erweitert 
werden, da8 aus den temperierten Quinten praktisch 
reine Quinten werden. 
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4. Streichinstrumentenstimmung und Tonsystem 


Die Streicher bemihen sich nach ihrer vielfach 
geauBerten Meinung und unseren Messungen fast 
ausnahmslos, ihre Instrumente moglichst rein zu 
stimmen, d.h. ihren Instrumenten reine Quinten, 
Frequenzverhaltnis 2:3, zu geben”. Damit besit- 
zen die Stimmungskurven von Klavieren bzw. Flii- 
geln die gleiche Steigung von etwa 3 cents pro Ok- 
tave wie die von Streichinstrumenten *. Das ist be- 
achtenswert, da Streichinstrumente praktisch die ein- 
zigen Instrumente sind, mit denen Klaviere zusam- 
men gespielt werden. 

Weitere Ursachen der Ubereinstimmung liegen 
darin begriindet, daf} das pythagoreische System — 
eine Folge reiner Quinten — melodisch bedingt ist 
und nach bisherigen Erfahrungen [6] Musiker 
(Streicher) bei Melodien eine der pythagoreischen 
angenaherte Stimmung bevorzugen. Auch bei 
ersten Versuchen mit geblasener Melodie waren die 
Intervalle gegeniiber der itblichen temperierten 
Stimmung erweitert, so da sich auch hier Beziehun- 
gen zur Klavierstimmung andeuten. Von Sangern 
horten wir ahnliches. Ein strenger Nachweis in die- 
ser Richtung nach bereits durchgeftihrten Messungen 
mit unserem Stimmgerat erfolgt demnachst von an- 
derer Seite aus. So diirfte die Ubereinstimmung zwi- 
schen den Stimmungskurven fiir Klaviere und 
Streichinstrumente zwar auch auf Stimmungserfah- 
rungen beim Zusammenspiel beruhen, kénnte aber 
sehr wohl allgemeinere Bedeutung haben, haupt- 
sachlich in einem gemeinsamen kunstlerischen Emp- 


2 Nur wenn man so stimmt, hat man den Eindruck, 
da8B die Geige wirklich stimmt, jede klangliche Unruhe 
innerhalb der Quinten deutlich wegfallt. Die Stimmung 
der Saiten in reinen Quinten bedeutet tibrigens, dal 
nicht etwa die Geige ,,rein“ gestimmt ist, das wiirde 
nur fiir d-dur gelten, sondern daB sie wegen ihrer auf- 
einanderfolgenden Quinten pythagoreisch gestimmt ist. 

3 Wenn das Cello-a, wie tiblich, mittels einer reinen 
Oktave auf das Normstimmton-a! eingestimmt wird, 
liegt seine Stimmung gerade um 3cents gegeniiber 
Violine und Viola héher (vgl. Abb. 21), um 50% mehr, 
als der Unterschied zwischen reiner und temperierter 
Quinte betragt. Das (und ahnliches bei anderen Musik- 
instrumenten) sollte bei Normvorschriften zum Ein- 
stimmen der Musikinstrumente beachtet werden. Auch 
fiir die Stimmung im Streichquartett ergeben sich Dis- 
krepanzen, lésbar durch Annahme _pythagoreischen 
Spiels. Doch kann auf dies alles hier nicht eingegangen 
werden. Hingewiesen sei lediglich noch darauf, da 
die ,,leeren“, d. h. unverkiirzt schwingenden Saiten der 
Streichinstrumente festliegende Orientierungstonhéhen 
sind, die der Streicher nie wesentlich verlaBt (und die 
auch die Frequenzen der anderen Klange mit festigen). 
Spielt er temperiert, spielt er daher zwangslaufig In- 
tervalle, die um etwa 3 cents pro Oktave erweitert sind, 
spielt also entsprechend der guten Klavierstimmung. 
Selbst wenn der Streicher nicht temperiert spielt, mu8 
die Beriicksichtigung der Stimmung der leeren Saiten 
seiner Intonation eine ahnliche Tendenz aufdriicken. 
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finden, das auch ohne Zusammenspiel zu den glei- 
chen Stimmungsergebnissen fiir die Instrumente ge- 
fiihrt hatte, begriindet liegen. Moglicherweise besit 
zen wir in der auf kunstlerische Bediirfnisse gegriin- 
deten guten Klavierstimmung eine Temperierung. 
den Ausdruck eines musikalischen Gefiihls, von gro. 
Berer praktisch-musikalischer Bedeutung, als der auf 
theoretischer Grundlage fuBenden Temperierung 
nach Werckmeister. 

Die instrumentenbautechnische Bedeutung det 
Aufnahme solcher Schwingungskurven sei kurz zu. 
sammengefaft:-Sie.geben ein genaues, allein auf 
meftechnischer Grundlage fuBendes Bild tiber die 
Kinhaltung des Normstimmtones und die dazu wich- 
tige allgemeine Stimmungshohe, iiber die Giite der 
Stimmung innerhalb eines Chores und von Taste zu 
Taste, iber die Haltbarkeit der Stimmung gegeniiber 
Transport, Temperatur und Spieleinfliissen. Sie zei. 
gen auch Mangel am Instrument selbst; bei Klavie- 
ren und Fligeln: zu schwache Rast, ungeniigend 
stabiler Eisenrahmen, zu locker sitzende Wirbel und 
anderes. Bei Akkordeons und Harmoniums geben 
sie u. a. Aufschlu8 iiber die Giite der Nietung der 
Stimmzungen, tiber Alterungserscheinungen und an.- 
dere Einfliisse auf die Stimmung. Nicht ganz un. 
erwahnt mag bleiben, da manche Stimmer iber 
ihre Leistungen eigene, nicht immer zutreffende. 
Meinungen haben. Auch solche Falle lassen sich klar 
und zugunsten der wirklich tiichtigen Stimmer ent- 
scheiden. Fiir Klavierstimmzwecke allein wird am 
besten ein Stimmgerat beniitzt, das die Erweiterung 
der temperierten Intervalle um 3 cents pro Oktave 
bereits beriicksichtigt. Ein solches Gerat wurde von 
uns nach Art des bereits beschriebenen [5] entwik- 
kelt. 

Fir die subjektive Beurteilung der Stimmungen 
danke ich Frau Prof. E:ty Ney, den Herren Ginter 
Kootrz, Wotreane Scuerericu, Prof. Frrrz Werrz- 
MANN, den Mitgliedern unseres Priifungsausschusses 
Tasteninstrumente: den Herren GrerHarp FO6rsTER. 
Fritz Grisster, Max Knopre, Wii Kronsere, Pau 
NeutscHMann, WALTER STILLER; fiir die sorgfaltige 
Durchfiihrung der Messungen meinen Mitarbeitern 
Frl. Erna Bacumann, Herrn Kurr Heroxp, Herrn 


Fritz KaIsER. (Eingegangen am 1. Juni 1956.) 
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Die Spiralaustastmethode zur Messung stationarer Schalldruckfelder 


von G. Hiner, Berlin 


Abgesehen von wenigen Standardbeispielen ist die 
Verteilung des Schalldruckes in einem akustischen Feld 
mit tragbarem Aufwand oft nicht zu berechnen. Aber 
auch die Méglichkeit, ein solches Schallfeld meBtech- 
nisch zu erfassen, ist heute noch recht unbefriedigend. 
Man erhalt auch hier nur durch eine grofe Zahl von 
Messungen und tiber langwierige Auswertungen den 
xewiinschten Uberblick. Das im folgenden angegebene 
MefSverfahren erlaubt eine schnelle, liickenlose Aus- 
messung groferer Schallfelder mit gleichzeitiger Re- 
sistrierung der Mef{werte. 


1. Verfahren 


Ein Mikrophon wird im Schallfeld auf einer Spiral- 
bahn so bewegt, da der radiale Vorschub pro Um- 
jrehung konstant und gleich der effektiven Mikrophon- 
breite ist. Dadurch wird eine Kreisflache liickenlos in 
einem Zuge ausgetastet, das Schallfeld in einer Schnitt- 
sbene untersucht. Die dem Schalldruck am Ort des 
Mikrophons proportionale Mikrophonklemmenspannung 
wird nach Spannung und Leistung verstaérkt und steuert 
wahrend der Austastung die Helligkeit einer Glimm- 
lampe. Die Glimmlampe als Registrierorgan wird 
schlieBlich synchron zum Mikrophon und im gleichen 
Sinne wie dieses auf einer verkleinerten Spiralbahn 
iber die photographische Schicht gefiihrt und belichtet 
diese dabei. Durch die Anordnung erhalt man ein Bild 


Abb. I. Das Austastgerat; a Mikrophon, b Verstarker. 


des Schallfeldes, wenn sich die belichteten Spiralen 
einander anschlieSen. Durch Einhaltung eines festen 
Abstandes zwischen Film und Registrierlampe und bei 
definierten Entwicklungsbedingungen ist das Verfahren 
quantitativ, man kann durch Photometrieren des Fil- 
mes nach einer Eichung den Schalldruck fiir jeden 
MeBort des Austastkreises in ubar angeben. 


2. Konstruktion, Daten, Schaltelemente 


Die Registriereinrichtung und der Endverstarker wur- 
den von dem im Feld befindlichen Austastgerat ge- 
trennt. Die Austastung geschieht mit konstanter Winkel- 
geschwindigkeit. Die hierdurch bedingte, unerwiinschte 
1/r proportionale, zur Bildmitte hin zunehmende Schwar- 
zung des Bildes kann elektrisch oder optisch bei der 
Aufnahme korrigiert werden. 


Winkelgeschwindigkeit m des Austastfliigels: 


Wahlweise 2,62s 1 entsprechend 25 U/min 
oder 0,873 s-! entsprechend 8,33 U/min; 
Durchmesser des ausgetasteten Kreises 2,61 m_ ent- 


sprechend 5,4 m?, 
Vorschub des Mikrophons pro Umdrehung: 2,6 cm, 
Dauer fiir die liickenlose Austastung einer Flache von 
emcees 2 2s) bzws O1's, 
Verkleinerung von Feld zu Originalbild: 1 : 13, 
Dynamik bezogen auf die registrierte Schwarzung 


(bei w=0,873 s-1): 65 dB. 
Kleinster mebbarer Schalldruck 


Der von der Apparatur noch mefibare kleinste Schall- 
druck ist von der Winkelgeschwindigkeit des Austast- 
fliigels abhangig und wurde hier einzig durch die bei 
der jeweiligen Winkelgeschwindigkeit auftretenden 
Storerschiitterungen des Mikrophons festgelegt. Bei der 
Ausmessung yon Feldern mit Stellen kleineren Schall- 
drucks mufs die Apparatur mit der geringeren Ge- 
schwindigkeit betrieben werden. Die weitere Bedin- 
gung, da die Winkelgeschwindigkeit nur so grof sein 
darf, da8 der vom kleinsten Schalldruck bewirkte 
kleinste Lichtstrom des Registrierorgans noch eine aus- 
reichende Schwarzung (0,1 iiber dem Schleier) hervor- 
ruft, ist allgemein erfiillt. Bei einer Messung von Fel- 
dern mit leisesten Stellen von 50 dB wird die Appara- 
tur der Stérerschiitterungen wegen mit der Winkel- 
geschwindigkeit #=0,873s 1 betrieben, kénnte der 
Schwirzung wegen * aber 10- bis 15mal so schnell ge- 
dreht werden. 


* Die auf Grund der Schwarzungsbedingungen hochstzu- 
lassige Winkelgeschwindigkeit bzw. héchstzulassige Drehzahl 
(Ymax) berechnet sich zu 

Ymax=43,3 Emin: 10%—D/10 d/R U/min; 


Ein kleinste Beleuchtungsstarke der Glimmlampe in lux 
(bei der Ziindung), 
n/10 DIN-Empfindlichkeit des Films, 


d Durchmesser des Lichtflecks, 
R Radius der gréBten Registrierspirale. 
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(a) (b) (c) 


Wandlungskanal 


Austastgerat: Kristallmikrophon mit der Rohre AC 701 
(in der ersten Verstarkerstufe) ; Vorverstarker mit 
den Rohren ECC 83 und EF 40; Schleifringe. 

Registriergerat: Endstufe; Schleifringe; Signalglimm- 


lampe der Firma OSRAM; Rontgen SS-Film. 
3. Niveaulinienmethode 


Die Nachteile der Registrierung bei dem angegebe- 
nen Verfahren im Hinblick auf eine quantitative 
Auswertung sind die genau einzuhaltenden Entwick- 
lungsbedingungen des Films und die Zeitdauer der 
noch notigen photometrischen Auswertung. Die Ermitt- 
lung der Linien konstanten Schalldrucks, also einer 
Karte des Feldes, erscheint eine anzustrebende Losung 
des Problems, ein stationaéres Schalldruckfeld quantita- 
tiv zu erfassen. Die gesamte Feldkarte lift sich in einer 
Aufnahme ermitteln, wenn man die Spiralaustastung 
mit einer solchen Spiralregistrierung kombiniert, die 
auf der photographischen Schicht nur dann Schwarzungs- 
punkte erzeugt, wenn bestimmte Schalldruckwerte — 
z. B. von 5 zu 5dB — vom Mikrophon aufgenommen 
werden. Dadurch werden die Linien konstanten Schall- 
druckes nach einer vollsténdigen Austastung punkt- 
weise zusammengesetzt. 


4. Schallfeldmessungen 


Ein Beispiel fiir die Spiralaustastung mit Intensitits- 
registrierung ist die Abb. 2 und ein Beispiel fiir die 
Spiralaustastung mit Niveaulinienschreibung die Abb. 3. 
In Abb. 2 wurde das Schallfeld hinter einer akustischen 
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(d) (e) 

Abb. 2. Parallelschnitte durch das Feld einer Zonenplatte. Abstand der MeBebene von der Zonenplatte bei der Aufnahme: 
(a) 35cm, (b) 46cm, (c) 80cm, (d) 155 cm, (e) 250 cm, (f) 350 cm. Wellenlinge 4 = 5,3 cm, Radius der ersten 
Fresnelzone 25 cm, entsprechend einer Hauptbrennweite von 255 cm. 
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Abb. 3. Schnitt durch das Feld einer im Freien stehender 
Tonsaule. Lange der Lautsprecherzeile 140 cm, Fre 
quenz 400 Hz, entsprechend 4 = 85,5 cm. Zwei be 
nachbarte Linien gleichen Schalldrucks unterscheider 
sich um 2,5 dB. 


Zonenplatte in Parallelebenen mit verschiedenen Ab 
standen ausgemessen und gleichzeitig registriert. Die 
Abb. 3 ist ein Versuch fiir eine Spiralaustastung mi 
Niveaulinienregistrierung, die noch nicht voll automa 


tisch geschrieben wurde. 
(Eingegangen am 21. Januar 1957.) 
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